Sapienza 

Università di Roma 


NetGrowth: simulating thè growth of biological 
neurons in spatial confinement 


Facoltà di Scienze Matematiche, Fisiche e Naturali 
Corso di Laurea Magistrale in Fisica 


Candidate 

Alessio Quaresima 
ID number 1471802 


Thesis Advisors 

Prof. Samuel Bottani 
Dr. Paolo Del Giudice 


Academic Year 2017/2018 


Thesis defended on 17 Aprii 2018 

in front of a Board of Examiners composed by: 

Prof. Nome Cognome (chairman) 

Prof. Nome Cognome 
Dr. Nome Cognome 


NetGrowth: simulating thè growth of biological neurons in spatial confinement 

Master thesis. Sapienza - University of Rome 

© 2018 Alessio Quaresima. All rights reserved 

This thesis has been typeset by UTpX and thè Sapthesis class. 


Version: Aprii 3, 2018 

Author’s email: q.ale@protonmail.ch 



Aux souris de Paris. 



Ili 


Abstract 

Advances in biological comprehension and technology allows for more complex in 
vitro experiments. C. Villard’s group at IPGG in Paris [Renault, 2015| as well as 
other authors [Smirnov, 2014] performed studies on thè environnrental interactions 
between growth cones and different substrates and mechanical confinements. The 
results of such studies are as relevant for medicai applications as for thè fundamental 
understanding of growth cone inner mechanisms. Furthermore, thè investigation of 
thè neuronal code brings researchers to study spontaneous activity in vitro, shedding 
tight on thè bursting phenomenon and it’s spatial properties [Orlandi et al., 2013a|. 
To reproduce such dynamics, a useful tool as NEST can be used - a powerful open- 
source simulator reproducing neuronal activity with eligible complexity. However, 
no tool is currently available to reconstruct thè network morphology or topology 
with a biological hint. 

The Neurophysics group of Université Paris VII-Diderot collected theories and 
results from biophysics and biology in order to combine them into an efficient 
neuronal growth simulator: NetGrowth. This simulator is built up with thè idea of 
modularity, allowing thè researcher to choose thè complexity level of his simulation. 
Since thè exact dynamics of actin, tubulin, and other constituents is stili obscure, 
thè simulator can act as a valid instrument to perforai in silico experiments. 

The application of such a tool is broad, ranging from thè microbiological lab- 
oratory, looking for a benchmark for in vitro experiment, to theoretical physics, 
requiring more realistic topology and morphology to study in silico dynamics. 

At this stage we tested some models for environmental interaction and elongation 
mechanisms, starting from biophysical theories and biological evidences. The navi- 
gation model results from an advanced study of randorn walk and run-and-tumble 
algorithms. For thè elongation we implemented and studied thè dynamics of a doublé 
SDE System, which reproduces thè perishable resources thè growth cone requires to 
stretch its skeleton and to elongate itself. The branching mechanism is simulated 
as well. We have implemented state of thè art models in an efficient manner. The 
relation among diameter and branching angles is reproduced as observed in thè 
laboratory. The environment sensing model - thè most relevant novelty of thè 
project - is implemented as a functional convolution of probability. The interaction 
between neurons which leads to fasciculation is stili a work in progress. Eventually, 
thè NetGrowth simulator is able to produce network topologies. The actual rnethod 
for synapse establishment is a naive estimatimi of neurite density, which was already 
tested for other environment-agnostic simulators. 

The manuscript describes in details thè state of thè art, thè study behind thè 
realization of thè simulator and thè obtained results. 
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Some notes on my internship at 
MSC Lab 


Before get into thè work itself, I want to submit some considerations: When I started 
this project one year ago I was an Erasmus student, with a curriculum in theoretical 
physics. My experience wasn’t net neither, indeed my Master Degree has been much 
more a survey over whole Physics than a proper specialization course. I get involved 
in this project at NeuroPhysics group without a clear comprehension of its purpose, 
I liked thè people, and that was enough. When I arrived in Paris, I needed a full 
rnonth to understand what I was expected to do, and when I realized I was scared 
by thè huge amount of work. I was used to program in C++, not to build user-proof 
interfaces. The aspect that baffled me most was thè biology of growing neurons. 
All thè knowledge in biology I had dated back to high school and thè course of 
Theoretical Neuroscience by prof. DelGiudice. I remember those days, exactly one 
year ago. When I expressed my doubts to Samuel, my mentor, he listened to my 
preoccupations and encouraged me, he said that it was not as hard as it looked, 
that, step by step, I would build thè simulator and be surprised by thè work and 
difficulties I had overcome. And so it was. It was hard, sure, I had to study an open 
problem like neuronal growth, study it through an incredibly vast literature and a 
wide set of point of views, discovering reviews and articles day by day, dropped off 
thè clouds when a new report contradicted thè one I was trusting. Yeah, it was hard. 
But I never felt thè fever of Science as in those days spent with Samuel and Tanguy 
discussing on thè future implementation of thè simulator. I have never learned so 
much about Science than during thè days reading neuroscience review. I am thankful 
to Samuel, Pascal, Tanguy and Stephane for thè conhdence they gave me. 
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Chapter 1 

Introduction 


Here we begin this journey in thè maze of neuronal computation. What is compu¬ 
tation! Is it some special feature of Intel produced chip-set or a basic property of 
matter? What happens in our brain when we sum up thè costs of a bill? Is it related 
to thè recognition process of a friend’s face? Where does thè computation happen 
in our brain? With this project we wanted to focus on a very sirnple problemi is it 
possible to reproduce such complex operations inside a 2D culture? If so, which are 
thè required ingredients? Is it thè network topology, its geometry, or thè strength of 
thè synaptic connections? And how does varying these elements affect thè properties 
of thè System? 

These questions are on thè edge among epistemology, artificial intelligence, and neu- 
roscience. The quest is on thè neuronal code i.e. thè actual method used by neurons 
to communicate with each other, which is a top-shelf question in neuroscience. Let’s 
start frorn thè basic ingredients to pavé thè way for thè NetGrowth simulator. 


Manuscript outline The thesis deals with thè design and realization of Net¬ 
Growth, a neuronal growth simulator based on biological evidence. The project is 
maintained by thè NeuroPhysics group, at thè MSC Lab, Paris Diderot University. 
The work is divided in three parts: The first part is an introduction to thè project 
itself, thè state of thè art is described and thè research project is motivated. The 
first chapter is an overview of biologically based computation, it is a survey over 
thè research project of thè NeuroPhysics group. The second chapter focuses on in 
vitro neuronal cultures, thè results obtained with different experimental techniques, 
some attempts of reproducing computational units with culture are illustrated. The 
third chapter presents some pre-existent simulators and list their pros and cons, thè 
relevance of NetGrowth in this panorama is discussed. 

The second part presents thè work we have done in thè last year, during thè intern- 
ship and in subsequent rnonths, to develop an efhcient simulator. The first chapter 
gets into details of neuronal microbiology, illustrates elements and processes we 
wanted to reproduce with NetGrowth. The second chapter reports some biophysical 
rnodels we have based our work on. Eventually thè third chapter is thè richest and 
thè most demanding of thè whole work, this chapter illustrates and evaluate thè 
model used for thè simulator, thè most of originai work I have done during thè 
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internship is collected in this chapter. In thè fourth chapter NetGrowth is evaluated 
and discussed, some complex scenarios are presented and future improvements are 
outlined. In thè conclusion I evaluate thè relevance of thè project and possible 
outcomes. 


1.1 Computation as an emergent property of networks 

1.1.1 Neural or Neuronal Networks 

The very fìrst concept to deal with is information. When Shannon introduced thè 
information theory in 1950 [?] he couldn’t imagine it would have affected so rnuch 
our contemporary culture. Indeed information is nowadays one of thè most prohtable 
and abused concept. Information is thè arnount of knowledge we have on a certain 
System, thè density measure of information, as introduced by Shannon, is its entropy. 
The theory was formulated in thè field of Communication Theory and aimed to 
formalize thè transmission of discrete signals; in this particular case, thè entropy 
describe thè uncertainty about thè next character of thè transmitted string. The 
concept of information was soon related to computation: thè ability of process such 
information, to extract relevant parts frorn thè whole, detect similarities and make 
evaluations. Living beings exert these actions all along their life. 

The generai quest about how thè living matters process thè information is on 
going debate, thè discussion is even pushed further to comprehend thè nature of 
consciousness. Anyway a meeting point between Neuroscience and Informatics are 
neural networks, an oriented graph of processing units. Neural Networks are driving 
thè Informatics research on Machine Learning and obtaining thè most striking results 
in reproducing intelligent behaviouiQ attracting thè attention of researchers and 
stakeholders. 


The common points between a formai System running on top of Silicon hardware 
and thè soft matter in our brains are thè network structure and thè existence of an 
activation threshold. A simple sketch of this parallel is offered in Fig 1.1 One of 


thè most important aspects of this analogy is that Computation appears not as an 
intrinsic property but emerges from thè web of connections among single, simple, 
units. 

Computation is not a privilege of biological existence although, computers are 
composed by transistors, and in thè very end, by a well-structured arnount of Silicon. 
The way these systems process thè information is well known, studied and improved 
day by day into researches lab. Processore, complex structures of hardly more than 
nanometres components, are among thè most advanced product of human technology, 
they are composed basically by an ALU (Arithmetic Logic Unit) and a CU (Control 
Unit) thè former operates those actions we usually regard as computation while thè 
latter control thè fìrst, shift information and rnanage thè whole process. 

From a mathematical point of view, thè work of Alan Turing (1912- 1954) laid thè 
foundations of computer Science, and fixed thè constraints to what machines can do 


1 The most famous was AlphaGo winning thè world Champion of Go, in 2016. 
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neuron celi body 



Figure 1.1. Neural or Neuronal Networks The parallel between Informatics formai 
structures or biological agents can be sustained on two main evidences. First both thè 
networks are built upon and with simple elements whose computational properties are 
minimum, they can usually sum a set of input and confront with a threshold. The 
inputs can be both excitatory or inhibitory. Second, ac ommon feature is thè non linear 
response. The all or none behaviour is defìned by an activation function cf>(x), which 
is a sigmoid for thè artifìcial networks and variate by neuron to neuron in biological 
tissue. On thè other end a relevant difference is thè time domain thè agents act on. For 
neuronal networks thè information appear to be coded in time also, while thè ANN are 
generally use discrete time and each evaluation takes one step. 
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or no10 in 1980’s thè introduction of quantum computer changed profoundly thè 
panorama, and new categories and tasks were defìned, maybe a better comprehension 
of biologically based information process will help to solve some of thè rnost puzzling 
questions of our century. 

1.1.2 Graphs in Neuroscience 

The relevance of graph structures appear in various beld of Neuroscience, I will 
shortly present some examples: 


Convergent or Divergent Connections The very interesting work of Martin 
Lindauer |Lindauer, 1996] starts from networks and synapses to reveal thè complex 
simplicity of animals communication. In thè account of sensory systems Lindauer 
depicts thè difference between divergent or convergent networks on thè sensor nerve. 
When thè signal from a celi is shared with upper standing cells thè information is 
reinforced and thè possibility it arrives to higher areas of thè brain increases. It can 
also be processed and refined, as in thè interneurons of eye’s nervous System which 
convolve many signals from thè retina to obtain thè contrasted image. The convergent 
configuration, instead, allows for more complex information: thè perception is 
augmented mixing different signals, e.g. thè bees antenna mixes olfactory and 
tactile signals for a plastic perception of fìavours. Fig 1 1.2 1 present a schema of thè 
connections types. 

Connectomics This field focuses on thè production and study of thè Connectomes, 
comprehensive rnaps of connections within an organism’s nervous System, typically 
its brain or eye. The human connectome can be viewed as a graph. Therefore, we use 
graph theory by its rich tools, dehnitions, and algorithms. By comparing diseased 
connectomes and healthy ones, we should gain insight into certain psychopathologies 
- such as neuropathic pain - and potential therapies. Generally, thè field of neuro¬ 
science would benefit from standardization and raw data. For example, connectome 
maps can be used to inform computational rnodels of whole-brain dynamics. Indeed, 
current neural networks rnostly rely on probabilistic representations of connectivity 
patterns. Connectograms - circular diagrams of connectomics - ha ve been used in 
case of traumatic brain injury to document thè extent of damage to neural network^ 

Integrated Information Theory When talking about information processing 
and biological intrinsic properties, UT plays thè most distinguished role in thè 
attempt to formulate a comprehensive theory of human consciousness. UT lies on 
graph theory and precisely on thè convergerne between graph theory and Information 

2 Touring introduced thè notion of Turing complete machines, and showed thè same, codihed, 
computation can be performed on different architec.tures c.hanging thè partic.ular set of procedures, 
thè algorithms. Together with Godei they set thè limit of such machines by thè very famous Halting 
Problem: for a undefined set of inputs a Touring machine can not say whether it’s computation is 
going thè end or not, these problems are said undecidable. 

2 http://www.scholarpedia.org/article/Connectome 
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Figure 1.2. Scheme for nervous System connections. The pictures reproduce thè 
scheme Lindauer observed in nervous Systems of insects and mammalians, thè sensor celi 
has two type of connections, usually both present in same sensory unit. The divergent 
connection ensure redundancy and signal transmission, while thè convergent connections 
gives rise to complex stimulus when merged into an interneuron. In thè right picture 
there are some examples of converging connections, each one with a specifìc function. 
This simple schema stresses thè relevance of graplr structure even in peripheral nervous 
System. |Lindauer, 1996| 
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theory. The topology of thè network is thè most important aspect. Roughly speaking, 
thè possible partitions represent thè possible states of thè experience [^j 


1.1.3 Morphology and computing dendrites 

The importance of graph topology Comes in action when moving to thè single unit. 
We have briefly mentioned thè simplicity of thè fundamental component of neuronal 
networks i.e. thè nervous celi. However, thè main attention of thè work in NetGrowth 
regards thè complex topology and geometry of dendritic trees. The reader who is 


unfamiliar with neuron biology will hnd a short review in Fig. 1.3 


A hint of thè complex process of arborization in offered in thè comprehensive 
review in [Cuntz, 2014j. The dendrite is investigated as an active actor of thè 
neuron-neuron signaling process, and its properties are outlined. This work was 
a good basis over which to build our simulator, and it is eventually to test thè 
result we are obtaining. The chapters we were most interested are those regarding 
thè optimization principles thè neurite satisfy during thè growth. The physical 
quantities which constraint thè neurite outgrowth are linked with its biological 
limits and scope: thè maximization of neuron connectivity given its limited set 
of time and nutrients [Cunt z et al., 201 0 , ||Tsigankov and Koulakov, 2009| , e.g. thè 
maximization of synaptic repertoire given thè dendritic arbor total length. The 
property of thè dendritic tissue are relevant for signal transmission and processing 
also; thè computation capabilities of individuai neurons have been suspected for a 
long time to be significantly richer in light of thè very typical dendritic architectures 
of some neuronal types. The presence of voltage-gated ion channels and action- 
potentials inside thè dendrites obliterates thè idea that they merely funnel incoming 
signals passively to thè soma. 

Interestingly, it was shown that - under particular conditions - dendritic 
spikes can compensate for thè distance-dependent attenuation observed in pas¬ 
sive dendrites, so that thè efficacy of synapses is nearly location-independent 
[V elte an d Masland.^1999 . The function of active dendrites could therefore be 
a simple way for neurons to free themselves from thè strong metric constraints 
inherent to physical networks. 

But a neuron endowed with excitable dendrites could do much more and perforili a 
complex non-linear summation of its synaptic inputs, depending on their spatial and 
temporal activation pattern. Although thè importance of active dendritic trees has 
been predicted for a long time, their role in computation has been experimentally 
observed quite recently. In thè retina, for instance, direction selectivity was shown 
to depend on active dendritic processing. Theoretical investigations showed that - 
under certain circumstances - even neurons with passive, non excitable dendrites 
can compute non trivial functions (like thè XOR logie gate), so that dendritic 
(pre) processing might be a generai feature throughout thè nervous System rather 
than an exception. 

In practice it might not be evident to exploit such computing power in artifìcial 
devices because it would presumably require precise positioning of synapses on thè 
dendritic tree of each neuron, which is fìrstly technologically challenging and secondly 


4 Nulla di più grande, Massimini - Tononi, Baldini + Castoldi 2013 
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intrinsically regulated through non yet elucidated mechanisms. The description of 
this problem lies outside thè scope of thè manuscript but underlies thè importance 
of a plausible reproduction of neurite morphology. A short résumé on well-known 
processing abilities is reported in Fig. |1.4[ 


1.2 What are neuronal cultures and devices 

Leaving well-behaved models and semi-conductors for thè realm of biology is easier 
said than done. Using living materials imposes a lot of constraints to thè experimenter 
who wishes to build something out of them. As everything that is alive can die, 
thè first challenge is to keep said living materials healthy as long as necessary. 
Providing a physiological environment with thè appropriate nutrients is of course thè 
main determinant for survival. However, like they do in vivo, cultured neurons are 
subjected to aging, giving as a consequence an additional layer of temporal constraints. 
The development and maturation of cultured neuronal networks actually involves a 
large panel of processes through which neurons first grow and establish synapses 
with other neurons, then express different sets of genes as tirne goes, causing some 
major shifts in their electrophysiological properties. Provided that thè basic needs 
of neurons are satisfied, and that they are in thè right developmental stage, there is 
another hurdle : homeostasis. Cells that we grow in vitro have originally evolved 
different compensatory mechanisms to maintain an internai state that guarantees 
proper operation across a wide range of environmental conditions, for thè sake of 
their common good, that is thè survival of thè animals they normally constitute. 
The cells will continue to fiercely regulate themselves when grown in a Petri dish, 
whether or not thè functional implications are relevant to thè experimenter. The 
sad corollary is that properties of interest that are incidentally connected to these 
homeostatic mechanisms can hardly be controlled experimentally, at least under 
physiological conditions. The microbiology labs are nowadays well-equipped with 
facilities for growing and studies neuronal culture on Petri dishes,A closer look to 
this experimental setup is offered in thè next section. 


1.2.1 Basic description of cultures: fabrication, properties, appli- 
cations 

A cultured neuronal network is a celi culture of neurons. It is used as a model to 
study thè centrai nervous System or to investigate fundamentals properties of thè 


living matter. A scheme is depicted in Fig. 1.5 The fabrication stages are generally 
based on soft litography, which enables for resolution at tens on nanometers scale, 
thè most common material is thè PDMS (Polydimethylsiloxane), an elastomeric 
organosilicon with notable visco-elastic properties. Soft litography has a wide 
set of advantages over thè traditional litography techniques and is widely used in 
biotechnology. The properties of PDMS include thè possibility to build sharp edge 
structures where neuron can be confined. Neuronal networks are typically cultured 
frorn dissociated rat neurons, because of their wide availability. Studies commonly 
employ rat cortical, hippocampal, and spinai neurons although laboratory mouse 
neurons have also been used. One of thè most formidable problems associated with 
cultured neuronal networks is their lack of longevity. Like most celi cultures, neuron 
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cultures are highly susceptible to infection. The long timelines associated with 
studying neuronal plasticity - usually on thè scale of rnonths - rnake thè extension 
of thè lifespan of neurons in vitro paramount. The applications of such devices 
range from medicai to fundamental research. Starting from thè laboratories which I 
personally visited or acknowledged, they include basic studies on axon guidance and 
interaction with Chemical and mechanical environment [^J graphene based neuronal 
tissue reconstruction or 3D growth from PDMS micro-pored structures |^J and full 
network studies as those performed on bursting spatial patterns Q The flaws of 
neuronal cultures are numerous, not even counting thè difference in terms of lifespan 
timescales. Cultured neuronal networks are by dehnition disembodied. Therefore, 
thè neurons are influenced in ways that are not biologically normal when outside 
their naturai environment. Among these abnormalities is thè foremost fact that thè 
neurons are usually harvested as neural stem cells from a foetus and are therefore 
disrupted at a criticai stage in network development. Lacking a body, neurons 
are also devoid of sensory inputs as well of thè ability to express behaviour - a 
cruciai characteristic in learning and memory experiments. It is believed that such 
sensory deprivation has adverse effects on thè development of these cultures and 
may result in abnormal patterns of behavior throughout thè network. Traditional 
cultured networks, 3D structure is a promising exception, are fiat, single-layer sheets 
of cells with connectivity only two dimensions. Most in vivo neuronal systems, to thè 
contrary, are large three-dimensional structures with much greater interconnectivity. 
This remains one of thè most striking differences between thè model and thè reality, 
and this fact probably plays a large role in skewing some of thè conclusions derived 
from experiments based on this model. 

The spectrum of accessible experiments and application has increased with thè 
introduction of more sophisticated techniques of production. The Microfìuidics has 
emerged as an important technology in academic research and even more remarkably 
in thè industry, to control thè architecture of neural networks growth in vitro. Using 
microfiuidic devices offers indeed many advantages over sirnple patterning techniques, 
including high manufacturability, robustness, controllable environment, etc... 

In thè next section a basic idea of tools and methods used to study neuronal culture 
is presented, then some higlight on thè activity of thè neuronal culture are reported 
and eventually this section is concluded with some notes on a particular type of 
culture: thè one produce at IPGG and whom I will refer often throughout thè 
manuscript. 

1.2.2 Instruments to observe and study neuronal cultures 

The investigation of neuronal cultures is performed with several techniques, whose 
perspectives and invasiveness change slightly. I had thè chance to study and, 
in some case, perform some experiments during thè Neuron Summer School in 
Trieste, two weeks of seminars and hands-on sessions focused on thè experimental 
techniques of neurobiology. A fully analogical investigation can be conducted in thè 
held of electrophysiology: measuring thè electrical potential at various resolution 

5 Villard, Ipgg lab, Paris R£nault^et^P_2016| _ 

6 Ballerini, Sissa, Trieste BosiujtjaPJlM et al., 2012| 

7 Soriano, UBV, Barcelona |Orlandi et al., 201 3a| 
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uncovers thè main communication pathway for neuronal networks. Electrophysiology 
can access thè voltage profile of a single neuron, or larger areas field potentials. 
These experimental methods are thè rnost long-lived in experimental neurobiology 
and where fundamental for thè formulation of neuronal dynamics models. The 
intrusiveness of such experiments depends on thè setup, but usually thè single celi 
recordings lead thè celi to death in thè range of a few hours. A wider perspective can 
be obtained with thè Multi Electrode Arrays (MEA), at thè cost of less resolution. 
MEA are grid structures of in/out electrodes. They allow for largely attenuated and 
temporally filtered signals as well as simultaneous recordings of large populations 
without celi labelling. The electrode diameter is 100 /un and a standard component is 
Silver Chloride, which converts ionie current into electricity. Reversely, it is possible 
to electrically stimulate neurons through thè electrodes. MEA is not recording thè 
internai properties but thè locai field potential (LFP). The LFP signal is blind to 
thè subtreshold area of thè neuron activity, and EPSP or IPSP-released present thè 
sanie LFP. Recent approaches consist in changing thè shape of thè array because thè 
piate shape has a limited coupling efficiency. When thè morphology, thè topology, 
or thè overall network dynamics is thè investigated property thè electrophysiology 
gives way to imaging techniques. The morphology, for instance, can be investigated 
by fluorescent microscopy, e.g. Audric ?? used SiR fluorescent dymers to evaluate 
tubulin amounts in outgrowing neurons. An exceptional investigation tool has 
recently been introduced with Calcium Imaging [Grienberger and Konnerth, 20121]. 
The Ca++ intracellular current occurring in spiking neurons proxies for a spatial 
localized pictures of activity. The Calcium Imaging is widely used to test for 
pharmacological response, for bursting patterns, and for numerous other experiments 
that associate single neuron labelling (i.e. by its position in thè culture) to culture 
dynamics on very large scales. Calcium Imaging can be combined with MEA for 
high quality spatial profiling of culture electrical response. There are severe flaws for 
this method too. First, thè Ca++ marking dimer decays over tirne and poisons thè 
celi. Second, thè measured signal is coarse in tirne and is attenuated on scales much 
longer than thè neurons’ spikes. The Calcium signal cannot eventually distinguish a 
inhibitory signal frorn an excitatory one. 

1.2.3 Activity in neuronal cultures 

When thè culture is constructed and thè neurons displaced on thè Petri dish, thè 
complex process of growth takes place. This process itself is detailed in thè thesis and 
thè scope of Net Growth is to reproduce it. The importance of electrical activity in 
thè shaping of thè network topology is well assessed in thè neurobiologist community, 
e.g. neurons die precociously when fading to connect, but this topic belongs to thè 
future perspective of this research project and will not be tackle. Here, we want to 
discuss thè electrical activity that animates thè mature culture, once thè growth is 
Corning to an end and before their feeble life ceases. The literature on this topic is 
wide, for simplicity I will concentrate on thè experiments thè group of J. Soriano is 
pursuing in his lab in Barcelona. 

Once thè most of synapses are established, and thè topology becomes more Constant 
in tirne, thè culture enters in thè bursting phase. This phase is essentially dominated 
by thè spreading of activation waves throughout thè culture celi; these waves can be 
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observed, e.g, with Calcium Imaging. The generation of these waves is investigated 
both theoretically |Fardet et al., 2018| and experimentally |Orlandi et al., 2013a 
and reveal many similarities with thè physical phenomenon of coupled oscillators: it 
is proved a sufficient condition for network bursting is thè presence of an excitatory 
population of oscillatory neurons which displays spike-driven adaptation, this result 
was proven on fully connected graph. The experiment conducted at Soriano Lab, 
UVB, Barcelona reveals that thè waves nucleate at specihc sites of thè culture and 
propagate very rapidly (with an increasing speed during maturation), and then are 
followed by long periods of silence or very low activity. More importantly, these 
sites are specihc (frorn culture to culture): there are preferred zones in thè System 
where waves nucleate, and these are not correlated with any locai properties of thè 
network (hring rate, density, connectivity, clustering, ... )(^J 

To solve thè puzzling phenomenon of bursting is an essential step towards thè 
realization of neuronal device and, more important, thè understanding of neuronal 
encoding of information. In this perspective thè use of activity simulators, e.g. 
NESTfJ assumes a fundamental role: it can help to test hypothesis and make 
predictions. Whether thè dynamics of fully connected graph can be easily performed 
thè reproduction of a realistic network has not be taken for granted, because thè 
actual adjacency matrix of cultured network remains mostly unknown. In thè next 
chapter we will see thè important role NetGrowth can fulhll in this context. 

1.2.4 Neural circuitry in Lab-On-A-Chip technologies 

In this section we will report thè work done by R.Renaud during his PhD in thè 
MSC Lab, focusing on thè studies and experiments he accomplished towards thè 
realization of an in vitro functional celi, also said a Lab-On-a-Chip. 

R. has focused on shaping thè networks at larger scales, organizing thè connectivity 
between groups of neurons instead of guiding single neurons. The organization 
of neurons on this larger scale can be controlled by direct physical confinement; 
this is generally achieved using micro-fluidics chips. They consist in rnolded Silicon 
polymer (polydimethylsiloxane or PDMS) covalently attached to glass coverslips 
through plasma- activated bonding. He proposed a new paradigm for neuronal 
devices: Primarly, microwells were used to compartmentalize cultures into neuronal 
populations, individually acting as super-neurons. These super-neurons are both 
robust and customizable as their properties emerge from thè mean composition and 
density of thè seeded celi suspension, which can be controlled accurately. They 
are also easy to read, as neuronal populations can display all-or-none bursts of 
activity which are easily detectable. Even if thè amplitude of such bursts can vary 
and therefore carry information relevant to neuronal computation [58], they can 
also be seen as stereotypical events corresponding to super-neuron spikes ; thè 
investigation of temporal codes is therefore straightforward in such scheme. Secondly, 
thè populations are connected through arrays of axon-selective micro-channels which 
control topographically thè initial connectivity between them. Since thè directed 
connections in thè network are spatially separated, each can be reinforced individually 
without affecting thè others. 


s Javier Orland website 
9 http://www. nest-simulator.org 
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This type of architecture can bypass thè limitations of previous approaches where 
learning is hindered by highly reciprocai and overlapping connections in thè culture. 
Additionally, network topologies developed in thè held of ANNs can be easily 
implemented using this approach (Fig. |ì~8j ) . The fact that thè network topology 
inside each population is not constrained in such devices means that a lot of 
computational power is potentially lost. In thè worst case scenario, and provided 
that populations can actually act as super-neurons, thè number of useful computing 
units would be thè number of populations, which is far less than thè actual number 
of neurons. However, random sub-circuits inside each population could be exploited 
to multiplex different functions on each unit, as it is done on unconstrained cultures 
growing onto MEA. The computational capabilities of population units are therefore 
difficult to predict, but they would at least support thè properties of thè units inside 
a perceptron. 

The learning protocols, Renaud carried through, failed for most of thè structures 
where neurons were cultured, but thè idea stays and always more groups and projects 
are studying thè possibility of complex and confined structures in cultures. The costs 
of creating a functional celi, from thè realization of thè PDMS structure up to thè 
learning protocol, unveil thè remarkable necessity of creating a cheap test bench. The 
simulator is, in this case, becomes a suitable instrument to explore thè fog and test 
culture shapes and learning protocol. Furthermore in thè context of axon guidance 
thè availability of a modular tool, to implement and verify chemio-mechanical 
interactions model, is very handy as well. In thè next section thè simulation of 
culture is detailed, and some other cases which can benefit of NetGrowth simulator 
are indicated. 

1.3 In silico neuronal cultures, applications and limits 

Whether thè complexity of life beings, thè multitudes of Chemical reactions and 
micro-physics process go far beyond thè computational possibilities of modern 
computer^] thè computers have revealed powerful tools at thè Service of neuroscien- 
tists |Amit, 1998| . Nowadays each research lab has an IT infrastructure, many have 
computing servers, and very often researchers are skilled in one or more programming 
languages. The possibilities opened by thè spread of IT technologies are numerous: 

• Non analytical theoretical models can be implemented and checked, e.g. simu¬ 
late thè propagation of an electrical signal on a complex dendritic tree. 

• The research can be performed on synthetic data, avoiding thè complex and 
expensive infrastructure needed to grow thè neurons, nor thè expertise to 
perform thè experiments. 

• The results become very reproducible and can be tested by thè scientihc 
community. Other researchers can recreate thè experiment and adjoin new 
research pathways. 

10 It actually depends by thè scale of thè process we want to reproduce, indeed thè folding of 
proteins can be simulated, either, but once at time and for just few seconds 
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In this section we will zoom on thè problem of reproducing neurons and neuronal 
cultures with computers. In thè first section I will have a look to thè perspectives and 
applications of neuronal simulations, stressing their relation with in vitro experiments. 
Then a brief review of existent simulators will be presented, focusing on those most 
related to thè NetGrowth simulator. 


1.3.1 Why simulate neuronal growth 


There are many reasons why creating effective and realistic neuronal morphologies 
can be interesting; they come frorn both neurobiology and theoretical neuroscience. 
This doublé importance is thè consequence of thè many scales observed in neu¬ 
roscience. Let us introduce both these quests and look at how thè simulation of 
neuronal morphology can pavé thè way to a deeper understanding of both neurobio- 
logical and cognitive processes. 

Regarding thè former, thè growth of neurons, and thè various morphologies they 
assume, is a puzzling question for biologists and biophysicists. The mechanisms 
that regulate thè growth, both Chemical and mechanical, are a subject of signih- 
cant interest and experimentation in modern research, as I had thè opportunity 
to understand during thè Neuron Technology summer school. The study of thè 
mechanical interactions between neurons and substrate opens thè doors to relevant 
medicai applications, e.g. in thè held of rehabilitation. The resulting morphology 
sheds light on thè Chemical process going on below thè celi surface and accounts for 
thè role of specifìc neurons in thè population: its ability to integrate signals and 
process Information depends mainly from thè dendritic tree. All this variable are 
kept into account in biophysical models as those presented in |Recho et al., 20161 , 
or |Hely, 2001 . Neuronal growth simulators offer thè opportunity to test these mod¬ 
els and observe thè consequences over thè whole morphology, opening thè door to 
effective comparison between theoretical model and experimental evidences, i.e. thè 
morphology observed and catalogued in morphology databases as Neuromor phcp^ 
When we are looking at thè neuron from thè top of cognitive basic functions, like 
spatial or working memory, thè morphology appear much less important. Its relevance 
is covered by thè multitude of neurons we are looking at, but it is thè fundamental 
process that provide for thè establishing of a neuronal network. Whenever thè 
researchers are interested only to thè adjacency matrix of thè network, neglecting thè 
synapses localization and thè neurite extension, thè simulation of realistic morphology 
is necessary. The function of thè graph, and thè macroscale properties of thè network 
depend by thè number of connections a neuron successfully establishes. A more 
concrete example is showed from three relevant publications on thè shelf of culture 
activity studies: all thè articles use a growth model to establish connections and 
generate thè adjacency matrix of thè culture. In case of |Orlandi et al., 2013b thè 
neuron morphology was almost absent and thè connections between neurons were 
established with a very abstract method, this is detailed in Fig. ??, thè generated 
network was used to simulate thè culture and thè Calcium Imaging measures 
compared with those in silico. The second example is more accurate: starting from 
NetMorpfp^l simulator, thè authors have studied thè variation of synapses density, 


11 http://neuromorpho.org/ 
12 described in next section 












1.4 NetGrowth, a brand new neuronal simulator 


14 


during thè outgrowth of environment independent neurons. The effectiveness of 
this simulations cannot be proven directly, but thè idea is more empirical: since it 
produces plausible independent morphologies thè juxtaposition of many of thern will 
produce something acceptable. 

Eventually thè project in | A Gritsun et al., 2012 starts from thè NetMorph simulator 
to build a more consistent network: here thè growth cones are driven by chemotactic 
gradients; details in ?? 

In thè next section some of thè available and published simulators will be described, 
with a short description of pros and cons. 

1.3.2 State of thè art in neuronal growth simulation 

The history of neuronal growth simulator has evolved in thè past twenty years, first 
examples were related to thè publication itself and thè implementation was very 
basic, usually thè model implemented was detailed in some aspects and neglect 
some others. After more sophisticated model were proposed thè simulator started to 
grow in complexity and took into account thè branching functions, thè environment 
and so on. A complete review of thè active simulator could appear as a boring 
enumeration, I have preferred to select thè relevant aspects of each one and stress 
thè new possibility thè simulator was offering. 

Here I will elencate thè simulators: +SelfWiringModel +Luczsak +Netmorph 
+NeuroMac 

1.3.3 NEST: thè activity simulator 

The overall flaws of thè simulator is thè impossibility to reproduce thè exact dynamics 
of thè living rnatter. 

1.4 NetGrowth, a brand new neuronal simulator 

In this section I will present thè NetGrowth simulator, I will describe thè collabora- 
tions from which it is promoted and point out thè new elements it introduces in thè 
panorama of neuronal growth simulators, an example is presented in Fig. |1.9| 
NetGrowth is a parallel simulator of cultured (2D) neurons, it simulates thè growth 
process, from thè first phase of neurite differentiation up to mature phase, a period 
of about three weeks. The neurons can internet with thè environment through their 
growth cones, thè neurite extruding tips. The user can set models and parameters of 
thè simulated cells; thè parameters can be changed during thè growth process itself, 
and, these features permit to reproduce empirical observation that are not comprised 
in any model. All thè models are correlated with biological evidences and biophysical 
models. The environment interaction objects to reproduce thè fasciculation between 
neurites and thè formatimi of synapses, based on thè dendritic spines evidence. The 
informatics infrastructure of NetGrowth reproduces closely thè one from NEST, 
thè aim of thè project is, indeed, to complement thè NEST simulator, towards an 
integrated simulation of activity and growth. 
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The NetGrowth project is thè result of thè collaboration among thè Neurophysics 
group at MSC Lab, thè group of Catherine Villard at IPGG, Paris, thè group of 
J.Soriano at UVB, Barcelona and E.Moses, Weizmann Institute, Rehovot. The 
research groups are involved in different stages of neuronal cultures studies, from 
thè axon guidance up to thè implementation of learning protocols for thè neurons’ 
population. I had thè chance to discuss about NetGrowth project with thè group 
of Villard and Soriano. They represent two of thè future, possible, users of thè 
simulator: thè former is looking for in silico studies to compare thè confinement 
strategies for neurons. The latter needs realistic network topology and nrorphology 
to base thè NEST simulations. I had thè opportunity to visit Catherine Lab and 
I pian to visit Soriano Lab in May, and, for thè occasion, help them to realize thè 
simulation of a new culture PDMS structure. The interaction with these groups is a 
splendid example of coordinated research: different experiences and competences are 
joined together at distinct level of thè investigation chain, in order to obtain more 
reproducible and consistent results. In this fraine we fulfil thè role of theoretical 
physics and bio informatics, translating experimental evidences into feasible models, 
that can be used for larger scale investigation. On thè other hand thè group of C. 
Villard, in thè person of Ian Audric, has furnished relevant and irreplaceable data 
whorn thè simulator is based on. 

With thè NetGrowth simulator, and thè NNGT module developed by Tanguy 
Fardet, thè scientific community of neuroscience is enriched of a sophisticated and 
performant tool. Starting from NetMorph, which is thè state of thè art in 2D and 3D 
growth simulation, our project adds consistent improvement. From thè scientific point 
of view thè insertimi of mechanical interactions rules, thè possibility of confinement 
and, in soon, of fasciculation and chemiotaxic guidance is a signihcant step forward. 
This improvement responds to thè technologically advancement in thè realization of 
micro-fluidic chambers. On thè hand of technical features it presents a change of 
paradigm: NetGrowth has a parallel infrastructure that responds adequately to IT 
resources and trends of today and tomorrow. Since T.Fardet, who I collaborated 
with in thè developing of thè project, has a role in thè Nest developers and users 
community, thè similarity to thè de facto simulator of neuronal activity are numerous. 
We hope to integrate thè projects in thè next future, opening thè door to very new 
perspectives: an efficient simulator of activity and growth, where thè electrical 
stimulation influence directly thè outgrowing properties and vice versa. Eventually 
thè simulator is realized in favour of a broader community of users in respect to 
its predecessors: while thè efficiency of thè compilation is granted from thè C++ 
language, with which thè models where implemented. The simplicity of use is offered 
by thè user interface written in Python, thè rnost used language in research labs. 
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Figure 1.3. Scheme of neuronal structure The picture presents a schematic neuron 
whit all its fundamental components. The neurons communicate with each other 
throughout a electrochemical signal which is originated in thè soma and travels across 
thè axon up to thè neurite of thè post-synaptic neuron. The dendrites transport thè 
signal towards thè soma and play a signifìcant role in amplifying or attenuating thè 
signal. About a million of synaptic connections can arrange over a mature dendritic tree. 



1.4 NetGrowth, a brand new neuronal simulator 


17 



Figure 1.4. The dendritic computational toolkit A schematic figure highligkting four 
key dendritic mechanisms, mapped onto a layer 5 pyramidal neuron morplrology, which 
can allow dendrites to act as computational elements. These mechanisms can coexist in 
thè same neuron and be active in parallel or in a hierarchical manner. 

Bottoni left: Passive dendrites act as passive filters. A high-frequency fluctuating current 
injected in thè dendritic pipette will evoke high-frequency and large-amplitude locai 
voltage responses, but thè response recorded by thè somatic pipette will be attenuated 
and smoothed (low pass filtered). 

Top left: Nonlinear interaction between excitation and shunting inhibition on small 
dendritic branches can implement logicai operations. The branch marked by an arrow 
sums up thè current from thè two subtrees, such that its output would be a logicai 
OR on their output. Each of thè subtrees will in turn inject current if and only if thè 
excitation AND-NOT thè inhibition will be active. 

Bottoni right: Dendrites can help reduce or amplify thè mutuai interaction between 
excitatory inputs. Excitatory inputs to thè same branch tend to sum sublinearly, whereas 
inputs on different branches sum linearly. Thus mapping input to different branches can 
reduce this effect. In neurons with active dendrites, however, clusters of inputs active 
synchronously on thè same branch can evoke a locai dendritic spike, which leads to 
significant amplilication of thè input. Synapses onto a different branch (open circles) 
are only slightly influenced by this spike. 

Top right: In layer 5 cortical pyramidal neurons, as depicted here, coincidence detection 
between thè apical and basai dendritic compartments is achieved by active dendritic 
mechanisms. A backpropagating action potential, which coincides with a distai synaptic 
input, will trigger a dendritic Ca 2+ spike, which depolarizes thè whole apical dendrite 
and drives a burst of spikes in thè axon. London and Hàusser, 2005] 
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Figure 1.5. Neuronal cultures a , Schematic representation of thè culturing prò cess by 
using pierced PDMS molds (blue) attached to glass coverslips (yellow). The preparatimi 
process (top) included piercing, autoclaving and poly-lysine coating; thè culturing process 
(center) yielded to thè formation of mini-cultures in thè pierced areas; thè final process 
(bottoni) consisted in thè removal of thè PDMS mold and thè preparation of thè culture 
for thè measurements. b , Actual image of combined glass-PDMS structure 11 days 
after plating. The rectangular area depicts thè maxi- mum held-of-view (FOV) of thè 
camera ( 8.2 x 6.1 nini), c , Bright fìeld image of a small region of a culture. Round 
objects are neurons’ celi bodies. [Orlandi et al., 2013a 



Figure 1.6. MEA culture dish and close-up look on thè electrodes. A weekold 
culture of circa 50000 neurons and glia from embryonic rat cortex, growing in a MEA 
and forming a dense network 1-2 mm across. Fifty-nine 30 /ira electrodes spaced at 200 
firn intervals connect a few hundred of thè network’s neurons to thè outside world, by 
allowing their activity to be extracellularly recorded or evoked by electrical stimulation. 
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Figure 1.7. High-contrast bright-fìeld image of a neuronal culture in CaH—H 

Grown on glass and confined within a circular diameter. The culture contains about 
3000 neurons. b Bright-field image showing a detail of thè culture and thè distribution of 
neurons. The circle identifles a single neuron. c Corresponding fluorescence image during 
a spontaneous activity event. Bright spots are fìring neurons. The resolution of thè 
image is thè same as thè actual measurements. d Fluorescence signal from a 30 minutes 
recording of thè spontaneous activity in thè culture shown in a, averaged over thè 500 
brightest neurons. The top plot corresponds to measurements with both excitation 
and inhibition active (E+I network); thè bottoni one corresponds to excitation-only 
measurements (E network), with inhibitory synapses blocked with 40 firn Bicuculline. 
Fluorescence peaks identify network bursts. The symbols below each burst identify its 
initiation in a specifìc area of thè culture. 
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Figure 1.8. Minimal neuronal devices Top : Two populations of circa 700 fluorescent 
neurons are connected through micro-tunnels, thè geometry of which can be controlled 
to yielcl unidirectional connectivity. Only a fraction of neurons are fluorescent here. 
Bottoni : Tliree populations device (two inputs and one output) showing neurons stained 
for thè cytoskeletal protein MAP2. Two and three populations devices constitute thè 
building blocks of more elaborate devices. 
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D Chandelier-like interneuron generateci with NetGrowth 
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Figure 1.9. A sample culture generated with NetGrowth NetGrowth simulator can 
simulate biologically relevant neurons in confìnement. The example presents thè neurons 
generated in a Petri dish like culture A. Neurons interact with thè environment and 
create complex shapes, merging thè forces imposed by thè mechanical traction with 
thè internai attitude to grow out from thè soma B. It is possible to evaluate whether 
two neurons intersect or not and hence create thè adjacency matrix for thè simulated 
network C. When thè neurons are not constrained they assume thè lrabitual morphology, 
e.g. thè neuron in D presents a Chandelier like morphology. 
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This part of thè thesis will deal with thè NetGrowth simulator. In thè first 
chapter thè biological process of neuronal growth will be described. Then we will 
recali thè biophysical models in scientific literature and eventually discuss thè models 
implemented in NetGrowth. 
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In this chapter thè preparatory work done for thè Net Growth simulator is pre- 
sented. The chapter is divided in two sections. In thè first we focus on thè growth 
cone behaviour, looking at thè interactions with thè environment, thè steering mod- 
els, and thè elongation rates. The second section gives a hint of thè complex world 
of neurite’s branching. 

But first, a short description of neuronal growth and its several stages is offered 
before entering in thè details of growth cone biology. 


Neuronal growth stages Let’s briefly describe thè stages of neuronal outgrowth 
to familiarize ourselves with this very complex biological process. 

During culture lifetime, neurons will extend their neurite and create synapses with 
other neurons. This process takes about 21 days. At this stage of thè project we are 
considering thè neuron fixed on thè substrate by thè use of some adhesive proteins. 
Neuronal growth occurs in several stages: 

first, a neurite ìnitiaiion process occurs in which hand-like extensions of thè celi 
membrane - thè lamellipodia - develop into short neurites with a rich actin frontier. 
This process is enhanced by thè actin-waves. 

After a period of outgrowth and retraction of these newborn neurites, one of them 
develops thè most and elongates rapidly over 2 days while inhibiting thè outgrowth 
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Stage 1: 
inltiation 


Stage 2: 
differentiation 


Stage 3: 




Figure 2.1. Stages of neuronal arborization The growth process is composed of three 
main stages of neurite development: (1) initiation, (2) differentiation into an axon and 
dendrites, (3) arborization, including elongation and branching. The whole process lasts 
3 weeks. Afterwards, thè growth slows down and a reshaping phase initiates during 
which most of thè dendrite is pruned. The topograplrical plasticity slowly disappears in 
favour of synaptic plasticity. 

of thè remaining neurites. This predominant neurite becomes thè axon and assumes 
distinct biological characteristics: thè neurites now have differentiated themselves 
into axon and dendrites. 

The dendritic arbor cornes out after a long period of elongation and branching 
of each neurite. At first thè arbor grows rapidly with respect to thè total neuritic 
length, and thè structure changes repeatedly. Eventually, thè elongation rate slows 
down and thè global shape of thè arbor becomes stable. 

The presence of obstacles, mechanical cues, or Chemical ones drastically changes thè 
shape and thè timing of thè arborization. The Growth Cones (GC) either show a 
repulsive or attractive behaviour with respect to different external stimuli. Such 
properties can be used to drive and control neuronal network formation. 

2.1 Growth Cone: searching for biological targets 

The Growth Cone is one of thè most amazing actors in thè outgrowth process of 
neurons. It is thè tip of thè neurite. It senses thè environment and assembles tubulin 
in micro-tubules in order to elongate towards its biological target. A sirnple sketch 
is offered in Figure [2~2j 

To describe all aspects of Growth Cone elongation and retraction would require a 
long discussion. Moreover, no description would be completely exhaustive, since 
there are rnany diverging opinions. The interested reader will hnd a detailed review 
in |Franze and Guck, 2010| and in |Lowery and Vactor, 2009] . 

Here, we will give an overview of thè principal forces acting on thè stage. The goal 
is to develop a sirnple model for thè Netgrowth simulator. Following thè review in 
[Suter and Miller, 201 1| , we’re going to present thè recent studies on thè elongation 
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Biological structure of a Growth Cone 



Figure 2.2. The structure of thè growth cone is fundamental to its function. The leading 
edge consists of dynamic, fìnger-like fìlopodia that explore thè road ahead, separated by 
Lamellipodia-like veils - sheets of membrane between thè fìlopodia (see thè figure). The 
cytoskeletal elements within thè growth cone underlie its shape. The growth cone can 
be separated into three domains based on cytoskeletal distributionl4. The periplieral 
(P)-domain contains long, bundled actin filaments (F-actin bundles), which forni thè 
fìlopodia. It also contains mesh-like branched F-actin networks, which give structure to 
lamellipodia-like veils. Additionally, individuai, dynamic, and ’pioneering’ microtubules 
(MTs) explore this region, usually along F-actin bundles. The centrai (C)-domain 
encloses stable, bundled MTs that enter thè growth cone from thè axon shaft, in addition 
to numerous organelles, vesicles and centrai actin bundles. Finally, thè transition 
(T)-zone (also called T-domain) sits at thè interface between thè P- and C-domains, 
where actomyosin contractile structures - called actin arcs - lie perpendicular to F-actin 
bundles, forming a hemicircumferential ring within thè T-zone33. The dynamics of tliese 
cytoskeletal players determine thè growth cone shape and movement during its journey. 
Images from [Lowery and Vactor~2009| 
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of axons, their response to traction, and thè influence of thè substrate on rate and 
direction of thè elongation. The approach will be bottom-up, starting from thè vary 
basic components of thè growth cones up to thè physical and computational models 
hopefully reproducing thè elongation and thè directionality of thè growth cone. 

For thè following sections I thank Ian Audric, who has written a precise review on 
growth cone internai biology, for his work at IPGG in thè group of Catherine Villard 
□ 

2.1.1 A pointed Little Head: biology of a growth cone 

The growth cones are composed mainly by 3 regions, thè Central, Transition and 
Peripheral domains. The first domain is dominated by thè presence of microtubules 
while thè rest by lamellipodia and filopodia, two actin filament structures. The 
macroscopic region containing thè P ans T domains is called kinetoplasm , or Growth 
Cone, while thè axoplasm connects thè GC to thè soma. 

The equilibrium between these regions, thè propensity to branch, thè microtubule 
assembling rate and thè direction are regulated throw hundreds of intracellular 
protein and extracellular signals, sensed and mediated throw thè filopodia, as shown 
in Fig. ?? Let’s get into details describing thè main actors in thè outgrowth process. 

The cytoskeleton is a very dynamic network of protein filament. It is re¬ 
sponsive for sustaining thè structure of thè celi. The growth cone is a particular 
arrangement of cytoskeletal filament, protruding from thè main cellular body. 

The neurite presents three main polymers: thè F-actin, thè micro-tubules, and thè 
neurofilaments [Coles and Bradke, 2015). The first two dominate thè short time 
scale which we are interested in when looking at thè dynamics of thè growth cone. 

The actin filament (F-actin) is a semi flexible polymer organized in a doublé 
helix structure. As explained in Figure ??, thè balance of polymerization and 
depolymerization is oriented. The actin is widely present in many cellular processes 
and it has thè ability of building chains of forces. It is fundamental in processes 
such as cellular motility or mitosis. It can be organized in a dense network or in 
(anti) parallel bundles. See Fig 
In thè axoplasm F-actin is mostly organised into a cortical mesh around thè micro¬ 
tubule’s inner core. The presence of myosin II motors in this cortex leads to thè 
presence of contractile stresses | Jiilic her et al., 2007|. In thè GC, F-actin is organised 
in a lamellipodium similar to thè ones found in cells specialised in crawling. 

The microtubules (MTs) are rigid polymers of a diameter of around 25 nm 
composed of tubulin,! - a proteic dirner. The tubulin is arranged laterally to create 
cylindric and hollow MTs. The MT has a persistence length of thè order of millimetre, 
as defined and measured in [Martin, 2013| . Its stiffness gives a robust skeleton to thè 
whole celi. It oscillates continuously between a phase of elongation and retraction, 
also known as catastrophe. The active equilibrium is finely regulated by a protein 
called MAP (Microtubule Associated Protein). 

They can be cross-linked by molecular motors (e.g. myosin II for F-actin, Dynein or 
Kinesin for microtubules) able to exert active stresses inside thè meshwork. Most of 
microtubules lie in thè axoplasm and are connected by passive cross-linkers (MAP) 
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5 min 


Figure 2.3. A biomechanic experiments to measure traction forces The growth 
cone produces different amounts of traction force during adhesion-mediated advance. 

(A) DIC images of a growth cone interacting with a 0.0025 N/m Con A-coated mi- 
croneedle at thè beginning of thè experiment, at thè end of latency and traction phases. 

(B) The same growth cone shown in (A) interacting with 0.0140 N/m Con A-coated 
microneedle. Needle stiffness (k), needle deflection, and traction force are indicated. 

(C) Kymographs showing thè position of thè microneedle throughout thè time course 
of thè experiments shown in (A) and (B). Arrows indicate thè end of latency phase. 
Images from Athamneh et al., 2015a] 


that generate a network with a quasi-lattice structure. These microtubules are highly 
stable with turnover duration of hours, possibly thanks to thè presence of MAPs. 

The research is active on this topic and technological developments allow biolo- 
gists to go further in thè understanding of thè complex behaviour of cytoskeleton 
polymerization and thè roles of thè many proteins active in thè neurite filaments 
Von Der Ecken et al., 2015 . 


Filopodia and Lamellipodia The growth cone contains an actin cytoskele¬ 
ton which adds mechanical strength to thè cone keeping its shape and as well 
helps driving and guiding thè cone’s movement. Filopodia are small actin fil- 
ament bundles extruding from thè tip of thè growth cone. They grow and re- 
tract at a high rate constantly probing thè environment and picking up guidance 
cue^J Filopodia can adhere to thè substrate, producing tension within thè growth 
cone |Mogilner and Rubinstein, 2005 Betz et al., 201 1[ |Craig et al., 2012] . 

In [Atham neh et al., 20l5b thè author offers a review over thè many roles forces 
play during neurite outgrowth. How filopodia sense thè environment is a relevant open 


Many videos are available on youtube . com 
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question in this research fìeld. Despite thè lack of understanding of thè mechanisms 
underlying mechanotransduction, experimentalists are using recent techniques - 
such as optical tweezers - to characterize thè response of fìlopodia to external 
mechanical stimuli. Fig. 2.3 shows an experiment with a coated microneedle, in 
which experimenters measure thè force exerted by thè Growth Cone on a microneedle. 


2.2 Elongation, steering and external stimuli 

Elongation is a key feature for neurites’ outgrowth. The Growth Cone senses 
thè environment during this process. Its ability to act depending on thè locai 
properties and thè external stimuli, is both thè most relevant ability of this biological 
micrometrie walker and thè main limit for thè neuron’s growth simulators. 

The direction of elongation and its intensity allow or not thè GC reaches thè target 
and fulfil its biological aims. This ability of navigating onto thè substrates and reach 
its target is widely studied in biology labs [^j 

The desire of controlling elongation and enhancing thè elongation process in vivo 
is very ambitious and offers a wide spectrum of applications, mostly in medicine and 
rehabilitation. Stili, thè mechanism of finding thè axon path is really complex, with 
hundreds of chemical |Wen and Zheng, 2006| and mechanical |Chua et al., 2014] 
cues interacting with a biological and active walker! Fortunately, thè development 
of bio-mechanical technologies, as thè AFM (atomic force microscope) or thè optical 
tweezers, are offering a wider spectrum of measures. They are cruciai to understand¬ 
ing this process and investigate thè mechanism underlying thè complex behaviour of 
axon and dendrites elongation. 

We will show how thè forces - intended as mechanical pressures and pulling - are get- 
ting a relevant place in biology and neuroscience and how thè possibility of quantify 
them is changing models and approaches [Athamneh et al., 2015b . Oversimplifying 
thè process, we can account for thè elongation with thè synthesis and fixation of 
tubulin in thè micro-tubule formation. The steering can be assumed as thè direction 
where such microtubules grow longer, as in Figure |2.4| Most diagrams of axonal 
elongation are, in essence, models of microtubule and actin dynamics occurring in 
thè growth cone as it grows in response to extracellular cues. The debate is stili open 
on which between thè pulling of actin cytoskeleton and thè pushing frorn thè back 
plays thè most significant role |Rauch et al., 20131. However, this simple scheme 
hides most of thè open questions in modern research on growth cone. In a recent 
article |Kahn and Baas, 2016 thè author wondered whether thè microtubules play 
thè main role in thè invasion of thè peripheral domain and in thè consequent turning, 
or if thè actin cytoskeleton is alone responsible of Growth Cone steering. This leaves 
us with open questions as: 

“What are thè means by which each relevant molecular motor protein is activated 
or deactivated in a manner that contributes to thè turning of thè growth cone in 
response to environmental cues?” 

3 http://neurofluidics.org/program/http://phdneurobiology.sissa.it/eng/faculty/ 
laura-ballerini.aspx 

http://phdneurobiology.sissa.it/eng/faculty/vincent-torre.aspx 
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(A) Nonpolarized axon growth 


(B) Growth cone turning 



Figure 2.4. Scheme of Growth Cone turning The centrai domain of thè growth cone 
is thè microtubule (MT)-rich region contiguous with thè axon shaft. The periplieral 
domain is thè outer-most part of thè growth cone. The peripheral domain comprises 
a broad fìat lamellar region in which actin fìlaments are arranged as a meshwork, as 
well as elongated thin fìlopodia in which actin fìlaments are arranged as aligned bundles. 
The transition zone is thè region between these two domains. Retrograde flow of thè 
actin cytoskeleton in thè peripheral domain pushes back most microtubules, compacting 
them in thè centrai domain. Individuai microtubules from thè centrai domain are able 
to penetrate thè transition zone to enter thè peripheral domain during growth cone 
advance. (B) During growth cone turning, microtubules extend from thè centrai domain 
through thè transition zone preferentially in one side of thè peripheral domain. Image 
from |Kahn and Baas, 2016|. 


The leading interpretation describes thè process in thè following manner. Poly- 
merization of actin fìlaments at thè periphery of thè growth cone and centrai 
depolymerization effectively result in actin bundles rnoving backwards through thè 
growth cone, giving thè growth cone a caterpillar track style movement. The space 
that is created behind thè growth cone by this pulling action is filled by microtubules 
and elongation occurs ||Kiddie et al., 2004] . 

The complex mechanism by which thè flux of actin determines thè synthesis 
of microtubule is too detailed for thè scope of this work, thè interested reader can 
refer to Bornschlògl et al., 2013 which offers a purely mechanic model, synthesized 
in Fig. 2.5 Following thè excellent work of P. Recho in [Recho et al., 2016], it is 
possible to distinguish whether they imply that thè GC pulls thè trailing axoplasm 
thanks to F-actin polymerisation pushing thè membrane in thè P-domain and myosin 
II contractility pulling from thè T- domain or whether it is rather microtubules which, 
from thè axoplasm, polymerise against thè T-domain and propel thè GC. Even whit 
such basic knowledge on Growth Cone dynamics thè relevance of interaction with 
environment come out. To study thè reaction to mechanic or Chemical stimuli leads 
to different research projects and questions, and thè literature is very broad on thè 
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Mechanical model for Growth Cone elongation 
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Figure 2.5. Mechanical model for fìlopodia. Force production via thè tip arises from 
thè parallel action of membrane forces ( F mem ) and from actin dynamics. Cortical 
retrograde flow ( v rc ) couples with high friction to thè filopodial actin shaft, modeled as a 
Kevin-Voigt body with an elastic modulus kfu 0 and a viscosity Cfn 0 . Polymerization at 
thè tip ( Vp 0 iy ) slower than thè cortical retrograde flow leads to retraction. Cytoskeletal 
force transduction can fail due to weak actin-membrane linkage at thè tip ( bi„.t )• An 
elastic force with stiffness k e ff controlled by thè optical trap is exerted on thè bead. 
Image from |Bornschlógl et al., 2013j. 


topic. 

Let’s now focus on understanding how thè mechanical stimuli offered from a 
wall or from a generic stiff object will infìuence thè behaviour of thè growth cone. 
This subject is commonly designated as mechanical signalling. Even in this case 
thè research is far from being concluded; thè details how thè different events are 
orchestrated to gradually build up traction force, up to hundreds of nanoNewton, and 
guide axonal growth in thè direction of thè adhesion site. In Athamneh et al.. 2015a 


thè authors show that traction force increases gradually over time as thè growth cone 
encounters a new adhesion substrate. The maximum level of thè generated force 
depends on thè stiffness of thè new substrate, implying continuous strengthening of 
thè clutch and/or active recruitment of molecular motors. From previous studies 
it is also known that thè growth cone response to a physically restrained adhesion 
substrate includes adhesion formation, SRC-tyrosine phosphorylation, slowing of 
retrograde F-actin flow, increased actin assembly, advancing of microtubules, and 
leading edge advance. 


2.3 Growth cones in constrained environments 

Here we discuss some recent experiments performed to study thè interactions between 
thè active growth cone and thè physical environment. We have seen that thè variety 
of behaviours and stimuli is vast; therefore we will focus on thè experiences closer to 
our laboratory and to thè airns of NetGrowth. 


Experiment on environment interaction 

Let’s now discuss some recent experiments performed to study thè interactions 
between thè active growth cone and thè physical environment, we have seen thè 
variety of behaviors and stimuli is spread so thè attention will be focused on thè 
experiences closer to our laboratory and to thè aims of NetGrowth. 

When talking of culture over substrates, and thè consequence of, a net example of me- 
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Figure 2.6. In (a) (top row) confocal micrographs show hippocampal cultures grown on 
2D-PDMS (left) and 3D-PDMS (right) immune-stained for b-tubulin III (in red), GFAP 
(green) and DAPI (blue), scale bar is 100 mm. 3D platforms exhibits properties and 
belraviours that differ from their equivalent 2D confìgurations, neverthless neurons that 
grow on supports, such as multi-layered artificial substrates, are only in part exposed 
to thè 3D environment, due to celi adhesion to thè platform(s). It’s stressed that 3D 
networks, besides displaying a similar degree of active cells when compared to 2D ones, 
display a higher rate of bursting, under both spontaneous and disinhibited conditions, 
due to an improved efficiency of neurons and neuronal connections in synchronizing their 
activity. |Bosi et al., 2015| 


chanically enhanced System is offered by thè work of Ballerini’s group at Sissa. They 
study cultures growing on a substrate of carbon nanotubes (CNT) |Gangaraju and Lin, 
Here, thè presence of CNT enhances thè growth of neuronal tissue and accelerate 
thè stabilization of synaptic cross-links up to thè standard connectivity reached in 
control cultures. In this case thè signaling is proved to be not only mechanic but 
electric too. Some applications of this experimental setup were in a rehabilitation 
frarne, helping thè reconnection of bone ganglia. 

The same group worked on 3D cultures built with transparent PDMS and stud- 
ied thè change thè third dimension carries on topology and geometry of thè net¬ 
work [R osi et al., 2015] 


an example of 2D and 3D cultures in Fig. 


?? 


2009], 


Another project we have collaborated with for thè NetGrowth initiative is held 
by C. Villard and collaborators. 

The group, located in thè microfluidic lab of IPGG in Paris, is interested in measuring 
thè modifications thè mechanical environment produces on Growth cone navigatimi. 
A first and seminai experiment was carried by R. Renaud, whom this thesis dues 
very much. He has measured thè affìnity of growth cones with PDMS structure, 
revealing thè formers’ attitude to follow thè structure edges. Such afhnity has been 
evaluated numerically as explained in Fig. ?? The group is also working on thè 
diode project discussed in Fig. EH and in generai on thè modifications that a 
narrow environment acts on thè growth cone behavior. They have implemented 
a selective “return to sender” strategy by exploiting this phenomenon, redirecting 
axons growing in thè unwanted “reverse” direction to their originai compartment via 
U-shaped lateral connections; and thè efhciencies of different variations of this basic 
concept regarding their direction selectivity, and demonstrated that they constitute 
a major improvement over thè tapered channels already used by thè brain-on-a-chip 
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Figure 2.7. Growth Cone affìnity to PDMS structures In this simple set of exper- 
iments R. Renaud studied thè tendency of a Growth Cone to follow thè edge of thè 
PDMS structure when it suddenly turns of a certain angle. The growth cones are 
highlighted with a tubulin bonded green dye. To measure adequately thè affìnity and 
thè percentage of adhesion is very hard with this setup, nevertheless it is possible to 
evidence a qualitative displacement from thè structure when thè angle is greater than 7r 


community. 

In a recent work, Ian Audric, exposes thè GC to a set of different shapes looking 
for common behaviors and aiming to a comprehensive description of constrained 
growth. Unfortunately biology is hard to forecast and thè experiment ends up with 
a too poor evidence to asses any result. The author is working thè project again 
and we hope thè NetGrowth simulator will benefice by his work in next rnonths. A 
picture of thè experimental set up is offered in Fig. 2.8 [p] 


As I have remarked until now, thè bibliography is really huge and I have done 
a big effort to summarize thè relevant aspects, both in this relation than in thè 
research work itself. 

I hope thè reader have a reduced but clear understanding of thè neurite outgrowth and 
growth cone navigatimi to comprehend thè next sections, where biology effectiveness 
is put aside for more affordable and computable rnodels of neuronal outgrowth. 


2.4 From single GC to Neurite arborization, branching 
and competition 

Why do neurons branch? How do thè branches relate to each other? Is there a 
main branch? Is branching similar in 2 and 3dimensions? Does exist, and which 
is thè optimal configuration for thè dendritic tree?All this question are very deep 
into thè comprehension of biological mechanism of neuronal arborization. This field 
has been highly investigated in thè past and is nowadays, thè understanding of 
neurite branching is far from be completed and while some biological mechanism for 
thè generation of new growth cones has been successfully reduced to their minimal 
components thè overall understanding of arborization process rnisses some relevant 
elements. To report thè evidences and offer a complete review of thè rnodels used to 
mimic neuronal branching will require a separate work. Instead we will introduce 
thè problem by thè work of |Gallo, 201 1| and |Cuntz et al., 2010 
. The objective of nervous System is to produce a functional information process 
System. The units of this System are thè neurons and neuronal population, how to 
connect each unit to others in an efficient way is thè challenge thè evolution had 
to solve. The human brain cells solved thè problem with thè branching process: 
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Figure 2.8. Axons grow in tiny PDMS microchannels, from I. Audric experi- 
ment at IPGG. (A) A growth cone advancing into narrow micro-channel. Between 
time zero and 2h thè GC (white arrow) retract itself and pause until a following GC 
joins it (red arrow). They keep on going forward witlrout further pauses. Another cone 
is visible at time 16h (black arrow). (B ; C) Self-rolling phenomena in GC encountering 
an obstacle. 
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Figure 2.9. Experimental results on a new kind of axon diode. A In this design, 
axons are expected to experience thè junction differently depending on their directions. 
The red axon grows along thè edge which undergoes continuous changes of orientation 
and makes a U-turn : this is thè blocking direction. 

The blue axon continues straight as it can not fold around thè corner, whose angle 
(in black) is smaller than thè criticai release angle : this is thè forward direction. B- 
Experimentally, axons growing from thè bottoni to thè top do preferentially go straight 
(left fìeld) or make U-turns (right fìeld) depending on thè relative orientation of thè 
junctions. 



Figure 2.10. A mechanical diode for neuronal cultures Intrinsic direction selectivity 
in microfluidic chamber experiment. The bottom chamber is seeded with neurons, where 
they are left to grow. The affinity between Growth Cones and PDMS edges forces thè 
axon to grow along thè walls of thè chamber, using this attitude, and thè persistent 
direction of Growth Cone navigation it is possible to impose, mechanically, a preferred 
direction of outgrowing. 

The showes thè correlation between thè number of archs in thè diode structure and thè 
probability of reach thè top chamber. Image from [Renault et al., 2016| 
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each neuron makes contacts with multiple targets through branching of its one axon, 
thereby reaching out to targets that are not in its direct trajectory, Fig. |2.11 presents 
an overview of thè wiring problem 


2.4.1 Evidences of branching in vitro 

There are basically two ways of undergoes a branching process: 

• Bifurcation of thè growth cone at thè tip of thè axon during axon 
extension. Growth cone bifurcation generates two axon branches from thè tip 
of thè extending axon. This process contributes to axon guidance and to thè 
development of thè basic organizational scheme of thè nervous System. However, 
branching through growth cone bifurcation is not thè major mechanism that 
contributes to thè sculpting of axons in their target fields. 


Formation of axon branches from thè axon shaft behind thè advancing 
growth cone. Fig. 2.11 B: Axon collateral branching - also referred to as 


interstitial branching - denotes thè de novo formation of an axon branch from 
thè axon shaft, independently from thè growth cone present at thè distalmost 
segment of thè axon. The formation of axon collateral branches is widely 
regarded as thè major mechanism used to establish axon arbors in synaptic 
target fields. 

We are not focusing on thè complex biological machinery to generate a new growth 
cone. The main idea, however, is thè presence of a higher amount of Actin-F. Such ex- 
cess can be caused by thè growth cone stalling or thè arrivai of an actin wave. Once a 
protrusion is formed, thè microtubules can invade thè bulge and forrn thè growth cone 
structure we have discussed before. Fig. |2. 12 offers a schematic picture of thè process. 


To conclude this short overview about thè biology of outgrowing neuron, it is 
necessary to discuss briefly thè problem of pruning. . In order to elongate thè GC has 
to assemble tubulin in microtubule, but thè building blocks it disposes are generally 
limited by thè rate of synthesis in thè soma [Goldberg, 2003| . Actually recent studies 
showed some protein are produced in thè neurite itself, but anyhow thè law of 
cytoplasm conservation by Ramon and Cajal has to be respected [?], establishing 
a competition arnong all extruding tips. Hence thè neuron has to adopt strategies 
to optimize its growth towards effective directions: exuberant axonal projections 
generated during development must be differentially regulated so that beneficiai 
branches are elongated while aberrant branches are eliminated |0§an et al., 201J]. 
This large scale axon degeneration has been documented and studied in a variety 
of developmental systems. Ollustrative examples of both branching and branch 
elimination are shown in Fig. |2.13 


2.4.2 Competition 

We have remarked thè evolution shaped thè genetics of neuron to fulfil their primary 
biological objective, create a network of information processing units. To reach 
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A Hypothetical Solutions to thè problem of connectivity 
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Figure 2.11. Modes of axon branching. (A) A single neuron could make contacts with 
multiple targets by: (i) projecting a single axon, guiding it alternately to each target; 
(ii) project multiple axons, each axon being guided to an individuai target, or possibly 
to multiple targets as in (i); (iii) project a single axon which undergoes branching and 
each branch contacts one or more targets. Naturai selection has continuously solved thè 
problem of neuronal connectivity by adopting thè strategy of axon branching (iii). (B) 
Growth cone bifurcation begins by thè loss of protrusive activity at thè leading edge 
wliile maintaining lateral activity. This results in thè formation of two independently 
active growth cones giving rise to two axon branches from thè tip of thè axon. In 
contrast, axon collateral branches forni de novo from thè axon shaft after thè growth 
cone has advanced past thè site of branching. Collateral branches are initiated by thè 
protrusion of fìlopodia or lamellipodia from thè axon shaft. Image from |Gallo, 20111 
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Figure 2.12. Rise of a collateral axon A basic sequence of thè reorganization of 
thè axonal cytoskeleton during branching. The fìrst step involves thè formation of 
actin-filament-based protrusions from thè axon shaft. The protrusions are subsequently 
invaded by axonal microtubules, giving rise to a nascent branch. Maturation of thè 
branch involves thè entry of additional microtubules and cellular components into thè 
nascent branch, e.g. organelles (not shown). [Gallo, 2011| 



Figure 2.13. Time sequences showing branching and pruning of dissociated 
dorsal root ganglion neurites. (a) Branching (red arrow) and extension (blue 
arrowheads) of primary axons. (b) Extension and retraction (blue arrowheads) of neurite 
tip. (c) Tertiary branching and pruning (encircled). |0§an et al., 20lT] 
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thè post-synaptic neuron is thè primary task of outgrowing axons, we can consider 
pruning a method thè neurons uses to optimize their neurites in this regard. While 
thè navigation should be considered a semi-local behavior -both thè chemotaxic 
gradient and thè mechanical confinement are sensed by thè GC itself, thè neurite 
influences thè Growh Cone with thè microtubule rigidity, and, maybe, realising 
chemorepulsive signaling - branching and elongation are thought to involve thè whole 
neurite or thè neuron even. Indeed thè law of cytoplasm conservation is showed to 
be - almost - enough to shape thè morphology of a dendritic tree [Cun tz et al ., 2010 


This constraint can be partially violated by thè outgrowing dendrite, and these 
violations are thè reasons of thè broad zoology of dendritic morphology. A similar 
idea is sustained in thè brilliant, thought old, review titled How does an axon 
grow ? |Goldberg, 2003| . Here thè relevance of thè question -Does thè axon requires 
building blocks from thè soma to grow or it’s self-sufficient?- is addressed and many 
experiments are reported to sustain thè two hypothesis. Beyond thè complexity of 
cellular biology sirnple questions can be drawn -thè newly created cellular membrane 
is produced? Do thè axon synthesize protein on its own or an diffusion/advective 
flux is received from thè soma?- Thus when simulating thè growth cone shaping 
thè dendrites we need to respect these generai constraints and broaden our look at 
thè whole neuron. In thè respect of these optimization principles we have settled 
from thè point of view of thè growth cones: they compete each other for thè limited 
resource thè neuron disposes. Some experiments were done to support these idea 
|Hj oriti et al., 20l4 |Singh and Miller, 2005j , and to fìnd out which are thè criticai 
nutrients that affect microtubule synthesis, and therefore, elongation. Stili is nec- 
essary to remark that thè growth cone behavior is pretty variable and to make 
assertion on it requires a huge amount of measures to obtain a result statistically 
valid. 


In thè next chapter thè straight observation of neuronal behaviour leaves thè 
way to a biophysical attitude for modelling, towards a simplified behaviour that can 
be reproduced with a computer. 
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Chapter 3 

Biophysical models for Growth 
Cone navigation and Branching 


The evidences discussed in thè previous chapter will now be presented in a biophysical 
fashion. The main difference is thè perspective with which this discussion is addressed: 
thè trade-off between exact description and mathematical feasibility is pushed and 
thè experimental evidences give thè way to reproducible models. In this chapter a 
bunch of already published models for GrowthCone navigation and neurite branching 
will be commented, these models were thè basis for our modelling work. The chapter 
is divided into three sections, first thè steering models will be discussed, then thè 
GC elongation and eventually thè branching process will be take on. 

3.1 Modelling thè growth 

A full model of thè neurite development and outgrowth, starting from thè positioning 
of neurons in a Petri dish, would require a proper model of neurite initiation 
and differentiation. However, we can consider thè axon distributed with uniform 
probability around thè soma because of thè lack of correlations between separated 
neurons and thè isotropy of thè space. 

Once thè first phase is over, thè neurites, leaded by thè growth cone, will be 
able to elongate and retract. The most relevant neurite is out of doubts thè axon, 
whose elongation rate is way larger then thè rest of dendrites, indeed thè most of 
thè articles concentrate on axon outgrowth. The study of dendrites and axon will be 
qualitatively thè same: any evidence of a substantial difference between axon and 
dendrites elongation mechanism is not reported |Suter and Miller, 2011 
In thè absence of a definitive, biophysical, and solved model, thè scientific literature 
offers plenty of particular models - often with a related in vitro experiment that 
tackle thè many faces of thè problem-. The following sections will review these 
modelsm highlighting pros and cons. 

3.2 Steering models 

It was very hard to find reviews or articles about thè steering mechanism. This 
basic process seems to be relevant only for computational neuroscientist willing 
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Figure 3.1. Memelli and Torben-Nielsen steering algorithm The algorithm used 
in |Memelli et al., 2013| is a superposition of directional vectors computed through thè 
parameters of thè different self-referential forces. 

The fìrst parameter, thè inertial force can be computed on locai information only, wlrile 
second and third forces require thè position of thè soma and thè neighbour cones. The 
former can be easily stored in memory but thè latter require full space perception at 
each step of thè simulation. 


to reproduce thè neuronal outgrowth in silico. The work we have kept in highest 
consideration during thè implementation of thè simulator is a previous project 
of Torben-Nielsen and Memelli |Memelli et al., 2013 . The authors suggested a 
environment agnostic model for neurite outgrowth: morphologies are generated by 
self-referential forces only. 

In this model, branching probability follows a Galton-Watson proces^] and thè 
geometry is determined by a drifted random walker, governed by internai forces. 
They modelled these forces, namely an inertial force based on membrane stiffness, a 
soma-oriented tropism, and a force of self-avoidance, as directional biases over a 3D 
random walk. 

These authors didn’t publish a detailed study of such forces, preferring to look at 
thè morphologies and benchmark over them. In Fig. 3.1 thè algorithm is described 
accurately. 


3.3 Elongation models 

In this section we consider three articles that most influenced our project. A preview 
of different models already used to describe thè elongation of growth cones. None of 
thè following articles clairns to provide a full model of neuronal growth. 

standard branching process introduced by Francis Galton’s for a statistical investigation in 
thè extinction of family names. 
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A mechanical GC Let’s start from a completely mechanical model as thè one 
proposed in |Q’Toole et al., 2008 . The axon is considered as a viscoelastic System: 
a series of Burgers elements, springs and dash-pot, as presented in Fig. 3.2 This 


System was used to model towed growth of thè axon, exerted with atomic force 
microscopy; thè assumptions are that thè axon reaction to force exertion enters 
in a steady state after few minutes and thè System saturates its elastic properties. 
By these assumptions thè elongation is due only to thè dash-pot and is completely 
passive, induced by thè pulling force; and it is possible to compute thè elongation 
rate. 

Assuming: 


1. In thè dashpot thè resulting speed depends linearly by thè exerted force, 

v = F/G 

2. The neurite contributes to thè elongation with all thè segments and not only 
by thè GC tip 

3. Some segments of thè neurite are anchored to thè substrate with a dash-pot of 
Constant ij 

<+-(-> then a simple System of differential equations can be solved and thè following 
equations state: 

L(t) = sinh(/3 exp{^t)) (3.1) 

The model forecasts a stable rate of elongation for each element of thè neurite, 
which ends up in a non linear growth rate at thè neurite tip. This model is incon- 
sistent with thè growth cones we want to model, which are not pulled but instead 
grow spontaneously. However thè model is quite simple and solvable analytically. 


MAP2 diffusion model A completely different model is proposed in |Hely, 2001 . 
Here thè author consider a compartmental neurite and thè diffusion of a certain 
protein on it. Such protein is thè Map2 which has a relevant impact on microtubules 
fasciculation and stability. The protein can phosphorylate or dephosphorylate [^] The 
activation (phosphorylation) of MAP2 enzymes reduce thè probability of branching 
and enhance thè elongation rate. In thè hypothesis 

1. thè phosphorylation is regulated by thè amount of Calcium in thè compartment. 

2. thè MAP2 protein is produced only in thè soma. 

3. thè Calcium influx depens only by neurite diameter. 

2 In biology, phosphorylation and its counterpart, dephosphorylation, are criticai for many cellular 
processes. Phosphorylation is especially important for protein function as this modifìcation activates 
(or deactivates) alinosi half of thè enzymes, thereby regulating their function 
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Figure 3.2. Model of a towed neurite as a series of dashpots. A Consider thè 
axon as a series of Burgers elements. Each element consists of two elastic elements 
and a free dashpot (with Constant G), which simulates towed growth. B Diagram of a 
neurite during towing. The distai region of thè neurite is free of thè substrate, whereas 
numerous adhesions in thè proximal region cause thè neurite to remain firmly fìxed. C 
Under Constant tension (Eo), thè behavior of a Burgers element is dominated by its free 
dash-pot. 

The model treats a neurite under Constant tension as a series of dashpots. Attachments to 
thè substrate are represented as friction dashpots (Constant h). Tension is Constant in thè 
distai region but is dissipated by adhesions in thè proximal region. IQ’Toole et al., 2008| 
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Map2 dynamics 



Figure 3.3. The celi dynamics simulated in thè branching model. Diffusible MAP2 U is 
produced in thè soma. First it is converted into thè microtubule-bound MAP2 b complex 
and then phosphorylated into MAP2 p . Calcium enters across thè membrane and diffuses 
throughout thè dendrite. The calcium-dependent functions F and G control thè on/off 
rate of plrosphorylation. The shape of tlrese functions is shown for various values of thè 
steepness parameter k used in thè simulations. 


thè author sketches a schema of 4 linear differential equations, reported in Fig. ??. 
Once thè MAP2 dynamics is integrated in a compartmental neuron, thè elongation 
rate E and branching probability P can be computed by thè following equations: 


E = E°k E MAP2 b MAP2b 

MAP2 p 

(3.2) 

P _ p0 /„ i\j a po MAP2 p 

P br - P br h B MAP2 p MAp2b 

(3.3) 


This model is very interesting for thè simulator project since it is quite omni- 
comprehensive and accounts for thè elongation rate as well for thè branching proba- 
bility. Nevertheless thè solution of thè linear System requires a compartmentalization 
of thè neuron and thè compartments number grows exponentially with thè branches 
in thè neurite. We have maintained thè generai idea of elongation as a byproduct of 
locai factors, thè Calcium influx, and soma generated proteins, using respectively 
neurite diameter and competition as proxies. The link between elongation and 
branching is also maintained in thè model inrplenrented in NetGrowth. 

Growth, collapse and stalling Collapse and stalling are introduced in a fully 
biophysical fashion in (Recho et al., 2016 . The author describes thè neurite as a 
viscoelastic System, imposing and solving set of conservation equations - conservation 
of mass, momentum and rheological assumptions. The neurite undergoes thè forces 
exerted by microtubule invasion of growth cone and F-actin meshwork contractility. 
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This model, which is based on only biological and measurable parameters - actin and 
tubulin meshwork friction parameters, polymerisation rate, etc... - predict success 
thè existence of three stages of neurite outgrowing. These three phases, growth, 
collapse and stalling, are regulated by thè load at thè end of thè neurite, generally a 
pulling force. Whether thè growth cone is free to elongate thè pulling force is thè 
actomyosin molecular rnotors that pulì on thè focal adhesion created by thè growth 
cone on thè substrate. Let’s dehne some variables: a = ^ is thè ratio between 

SM 

thè friction of microtubule with thè neurite cortex and thè friction of cortex with 
thè substrate; hence thè Growth Cone exerted force is Q gc = y/av p , where v p is 
thè polymerisation rate of G-Actin in F-Actin. The derived equations are quite 
complicate but can be summarized in three regimes: 

• Retraction, whether Q gc is less then thè stress exerted by thè microtubule 
turnover Q p 

• Stalling, if Q gc is comparable with thè friction of thè neurite cortex Q s and 
microtubule turnover, and thè total length is predicted in this case. 

• Elongation with velocity 

Y _ Qgc ~ Qs 

2-\/( a ) + B 

where B comprehend parameters we have not introduced, generally linked with 
neurite width and actin viscosity. 

This model is very well done, both for thè deepness of computation, where each 
assumption is justified with experimental evidence, both for thè overall simplicity of 
thè predicted results. To conclude its study thè author has performed pharmacological 
studies, in very good agreement with thè predictions. Fig. |3.4| summarizes thè resultss 
We have replicated thè existence of this three thresholds in NetGrowth, but there 
is not an equivalent model for thè stresses and frictions exerted by thè neurite 
components. We hope to implement this model and verify thè reproducibility of 
Recho’s results in a future model of NetGrowth. 

3.4 A mathematical picture of neurite branching 

Let’s now tackle thè complex problems of GC generation. The initiation of branches 
-and thè consequent arborization- shapes thè neurites, their computational properties 
and thè number of synaptic connections. It is fundamental to reproduce thè neurite 
tree adequately. 

Several modelling approaches have been used in thè in silico synthesis of dendritic 
trees. These can be characterized either as Growth Models or Reconstruction 
Models. Growth Models are based on principles of dendritic development, using 
rules of outgrowth associated with dynamic growth-cone behaviour and microtubule- 
mediated neurite elongation ??. In contrast, Reconstruction Models use an algorithm 
based on a canonical set of elementary properties which are originally derived from 
thè characterization of existing dendritic structures |Polavaram et al., 2014). Note 
that Reconstruction Modeling is a purely descriptive approach which uses minimal 


(3.4) 
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rules to synthesize topologically-realistic neurons. 

In contrast, GrowthModeling adopts an exploratory approach by using biological 
rules of development and observations of thè outgrowth process to explain or predict 
variations in full-grown arbor structures. 

These two modes of branching give rise to different geometrical patterns of axon 
branching and generate different aspects of neuronal circuitry. For thè NetGrowth 
simulator thè choice was obligatory, a full generative model was due to mimic thè 
wall-GC interactions. 


BEST model Altough thè biological mechanisms and thè relevant parameters were 
not fully understood, Dutch researchers Van Pelt and Van Oyen gave a great contri- 
bution to thè understanding of outgrowing neurites. They have published several ar- 
ticles |van Ooyen et al., 2014 Van Pelt and Uylings, 20031 |Van Pelt et al., 1997| to 
face this problem with a simple and effective mathematica! model. In |Van Pelt and Uylings, 2003 


they propose thè BEST model: thè branching nrechanism is GC splitting only and 
thè parameters are completely phenomenological. Given thè number of segments 
n, thè competition parameter among GCs E and a baseline distribution for thè 
branching rate thè growth cone branches with probability: 


Pi(t\n(t)) = Nn~ E 2~ liS D(t) 
with thè normalization factor 

n(t) 

a r = J2 21iS 


(3.5) 

(3.6) 

(3.7) 


It is possible to show, for an exponential decaying baseline, thè number of branching 
events is limited: 

ng(t) = r 4o~ e Ì/T f° r E>0 (3.8) 

Here 7 \ is thè centrifugai ordei)^] and thè parameter S is a symmetry parameter a 
sketch of such process is offered in Fig. |3.5[ 

BEST model is an artificial solution, where parameters have to be fitted with 
experiments and can variate slightly from a neuron type to another. The converging 
integrai dtD(t ) = AT can be associated to thè average number of branching 

events B, while thè time scale is fixed by thè characteristic time T. The parameters 
optimization is performed considering thè mean and thè variance of measured reai 
neurons, as those retrievable in thè NetMorpho database. 


3.4.1 Many GCs compete with each other 

The pruning introduces thè problem of competition for thè growing neurites. Which 
branch is going to collapse? Is thè optimization process of thè tree structure locai or 
global? Is possible to model arborization as a competition between growth cones? 
The microbiology of neurotransmitters that trigger this process is a complex subject, 

3 The branching baseline is necessary to give a temporal trame to thè BEST model: thè fìrst 
implementation [Van Pelt and Uylings, 2003] , used discrete distribution and Constant time binning, 
without any relation to time. Such a model whether reproduced thè right statistical measures were 
wrong in terms of branching time. 

4 Centrifugal order is thè number of nodes separating thè node i from thè root 
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Figure 3.4. Growth, collapse and stalling from biophysical quantities The model 
proposed by Recho, and summarized by Eq. |3.4[ starts from a sirnple depiction of thè 
growth cone internai forces, thè one proposed in A and develop a set of viscoelastic 
equations. After some evidence-based assumptions thè model is simplifìed and a two 
thresholds are derivate: a stalling threshold and a retraction one. The fìrst depend by 
thè microtubule turnover rate, while thè latter by thè contractile load of actin cortical 
network. Picture B show thè growth, stalling and retraction regions varying thè pulling 
load, i.e. growth cone force in our case, and thè parameter e, this is thè ratio of acto- 
myosin molucular motor over thè microtubule network viscosity [Recho et al., 2016] . 
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Figure 3.5. VanPelt BEST model The picture represents thè time evolution of two Van 
Pelt trees. The Constant baseline function - namely D(t) = B/N - is thè average nurnber 
of branch events. The sketch shows two different sequences, where thè sanie stochastic 
model generates two different trees. The time component D(t) is introduced to fìt reai 
neuron measurements. Switching from a discrete-time to a continuous-time probability 
function impaets thè theory and prevents thè nurnber of branching events from divergence 
thanks to thè exponential clecay of thè function D(t) |Van Pelt and Uylings, 1999] 
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Figure 3.6. Tubulin dynamics in thè model. Tubulin molecules (green spheres) are 
produced in thè soma by biological neurons via translation of mRNA on ribosomes 
(thè brown structure). Tubulin is then transported by diffusion and active transport. 
In biological neurons, thè microtubule bundles (thè light green fìbers) act as railway 
tracks on which thè tubulin molecules are bound via motor proteins (thè red molecules). 
Tubulin is transported to thè growth cones at thè tip of thè neurites. At thè growth 
cone, tubulin is either integrated or polymerized into thè microtubule cytoskeleton (long 
polymers of tubulin dimers - green fìbers), which elongates thè neurite. When thè 
microtubule depolymerizes, thè neurite retracts and tubulin becomes free again. 


but we can give a generai idea of it by a simple mathematical model. In fact, thè 
competitive dynamics introduced by |Hjorth et al., 2014] offers a simple perspective 
over thè whole phenomenon. Some elements - required by thè GC to synthesize 
microtubules and then elongate - are built into thè soma, thereby needing to 
be transported to thè very end of neurite’s tree in order to be available for thè 
growing GC. This transportation can be diffusive or active - i.e. done by molecular 
motors - and thè two cases lead to very different results. This simple idea was 
implemented with a compartmental simulator, where each segment of thè neurite 
stores or consumes some of thè limited material. In thè cited work, thè referred 
protein is tubulin. Fig. 3.6 presents a scheme of its production and polymerization. 


During retraction tubulin is depolymerizated. 

The equations of thè model can be condensed in thè following one that represents 
thè evolution of thè amount of tubulin Qi in a certain compartment: 


^ = Df(Q, A, V ) + ug(Q, A,V)-bQ-X^ 


(3.9) 


thè equations govern thè evolution of thè tubulin amount Qi\n a certain compartment, 
thè r.h.s first terni is thè diffusive one and is moderated by thè diffusion coeffìcient 
D, thè second terni represents thè active transport regulated by v, b is thè tubulin 

corruption rate and, eventually, thè last term X — accounts for thè consumption 

at 

of tubulin in thè growth cones. The tubulin diffusion ODE is associated by an 
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elongation ODE for each terminal segment: 

d ^-=pQ-q (3.10) 

The elongation rate is defined as thè result of thè polymerisation rate p times thè 
tubulin Q in thè leave, minus thè depolymerisation rate q. An experiment with this 
model is detailed in Fig. |3.7[ 

The model is quite sirnple but requires thè integration of thè whole ODE step by 
step and is not very suitable for an efficient simulator. The competition model for 
criticai resource we will introduce later is largely inspired by this model. 


3.4.2 Optimal wiring and scaling laws 

Competition among growth cones can be by another perspective, like thè one 
proposed by Cuntz and collaborators|Cun tz et al., 2010[ Cuntz et al., 2012] . Instead 
of generate thè neurite step by step through growth cone migration, thè authors 
of TREESToolbox [^] have proposed an optimization-driven approach for neurite 
reconstruction: thè algorithm builds tree structures which minimize thè total amount 
of wiring and thè path from thè root to all points on thè tree. On thè basis of Ramon 
y Cajal cytoplasm conservation law, they propose three more detailed principles: 


• By adjusting thè balance between thè two wiring costs -total dendrite length 
and distance from thè soma-, a dendrite can efficiently set its electrotonic 
compartmentalization, a quantity attributable to computation. 

• The density profile of thè spanning field in which a dendrite grows determines 
its shape dramatically. 

• A few weaker constraints such as thè suppression of multifurcations, thè 
addition of spatial jitter or thè sequential growth of sub-regions of a dendrite 
are helpful for reproducing thè dendritic branching patterns of particular 
preparations. 


Unfortunately there is not any direct application of these additional constraints in 
a generative simulator. Nonetheless they might shed light on further functional, 
computational, developmental or network determinants for certain dendritic struc¬ 
tures, we have tried to take in to account these constraints in thè competitive model, 
e.g. setting an explicit attenuation of thè criticai resource received when thè -user 
defined- total dendritic length is exceed. Most of all thè sirnple principles presented 
in this study can be used to efficiently edit, visualize, and analyze neuronal trees 
and will be used to validate thè morphology generated. 

On this purpose thè result accomplished in (Cuntz et al., 2012] is even more rele- 
vant: although thè wide diversity of dendritic trees, they have developed a generai 
quantitative theory relating thè total length of dendritic wiring to thè number of 
branch points and synapses. The optimal wiring presented above predicts a 2/3 
power law between these measures and they show thè theory is consistent with data 
from a wide variety of neurons, across many different species. 

’http : //treestoolbox . orgTREESToolbox 
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Figure 3.7. Competition between neurite branches in a complex morphology. 

(A) Example morphology of a reconstructed pyramidal neuron with apical and basai 
dendrites. (B) In thè control case, starting from thè reconstructed morphology, thè 
neuron was allowed to grow out for 10 hours in thè model. The simulation was then 
repeated with thè sanie initial conditions, but with an increased polymerization rate 
for one of thè growth cones. The dendritic morphology obtained in this last simulation 
is represented by a dendrogram, coloured according to thè tubulin concentration in 
thè branches. The gray vertical lines at thè terminal segments indicates thè starting 
morphology, and thè black vertical lines show thè neurite’s length after 10 hours - in 
thè control case. The black dot marks thè growth cone with increased polymerization 
rate. (C) The competition between branches increases with increasing path distance to 
thè soma. The graph shows thè total retraction of all neurites divided by thè growth of 
thè modified growth cone, as a function of thè path length between thè modified growth 
cone and thè soma. |Hjorth et al., 2014 
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3.4.3 Branching angles and forces 

3.5 Gateway to stochastic models 

The discussed models are deterministic as a consequence of solving thè macroscopic 
equations of mechanical or Chemical processes. However thè simulation of such 
complex events cannot be fully deterministic since thè phenomenon undergoes non- 
equilibrium statistical mechanics. A weak argumentation is to consider thè many 
undergoing processes we are neglecting, proteins diffusion, synthesis, and thè chain 
effect on thè growth itself; thè branching itself is depicted by a System with a variable 
number of ODEs (equation are doubled after first brandi). This is enough to expect 
a deeply stochastic process and to introduce stochastic elements into thè simulation. 
The Net growth simulator models thè out growth dynamics of axons and dendrites 
with stochastic rules, at many levels. The steering is regulated by a random walk, 
thè branching is mediated by a stochastic model and thè neuron position is drawn 
by a spatially uniform distribution. 
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Chapter 4 

Implemented models in 
NetGrowth for neuronal 
outgrowth in constrained 
environments 


Most of researchers focus into study a particular aspect of growth cones migration, 
e.g. thè elongation rate neglecting thè steering mechanisms and vice-versa. Although 
this rnethod is scientihcally correct, complex simulators, as NetGrowth, requires to 
mix these models to reproduce a plausible neurite outgrowth. Purthermore some of 
thè algorithms and mechanism discussed are computationally expensive, e.g. thè 
compartmental model which is widely used to describe thè growth cone elongatior0 
Hence thè ideas are kept and thè algorithms rearranged to fit with NetGrowth scopes. 


Outline of thè chapter In this chapter we are going to describe in details thè 
models used for thè NetGrowth simulator, which I personally developed and tested. 
We will tackle each problem separately: thè steering process, thè interaction with thè 
environment and thè elongation rate will be discussed. Some of thè largest digression 
are rnoved in thè Appendixes, to keep thè text linear and coherent. The chapter is 
divided in two sections, in thè former I will look at single growth cone migration: thè 
way they interact with thè environment and move towards their biological targets: 
thè isotropie model is discussed previously and then thè environment-aware one. 

In thè latter section I will introduce thè whole neurite: first a stochastic equation for 
neurite branching is presented, later thè competition among growth cones of same 
neurite is studied with a System of two SDE. 

1 It requires thè integration of an increasing set of ODEs: thè integration time goes at least with 
thè number of segments in thè tree, which is exponentially growing in time 
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4.1 Growth Cone migration: Steering as a random walk 
in 2D 

The easiest way an agent can explore thè space is carrying out a 2D Random Walker, 
a Markovian process largely used in physics. We want to verify if thè RW reproduces 
adequately thè navigation of thè Growth Cone on thè substrate. The hypothesis 
is thè biological procedures to steer and explore thè environment can be mimicked 
with a proper algorithm. Hence thè algorithm is expected to simulate thè migration 
of a Growth Cone in a physical constrained environment. 

In thè hrst section I will introduce thè problem and offer a generai picture of thè GC 
dynamics in thè NetGrowth simulator. In thè second section thè relevant properties 
required to thè algorithm will be presented, some algorithms will be outlined and 
thè results discussed. Eventually some reai neurons will be analysed and compared 
to thè simulator. This short section focused a relevant part of thè project: thè 
navigation is a fundamental aspect of thè whole simulation. Some tools, theoretical 
and computational were used to characterize thè algorithms, they are discussed in 
Appendis[A] 

4.1.1 Modelling thè elongation with a random walk 

The random walk is very easy to simulate and code, to associate it to a Growth 
Cone is straight-forward: thè active walker is thè GC while thè previous location are 
joined up to forni thè neurite. This simple fact implies thè walker cannot cross itself 
constantly and thè locai aspect of thè trajectory is relevant. Indeed in physics or 
fìnance thè random process are often 1D and in generai thè history is not relevant. 
A larger category of walker are thè self-avoiding paths, but they have not been 
used in thè present work for they are computationally expensive. Some classical 
algorithms, whose analytic feasibility make interesting, like thè Ornstein-Uhlenbeck 
process, were set aside for they are definitely not self avoiding and thè trajectory 
cannot resemble anything biologically consistent 

The Growth Cone moves on thè xy piane but it’s evident thè random walk is actually 
on thè angle. For a subset of algorithms thè step can be hxed and thè whole process 
depends by thè distribution of new angles. In this case a simple consideration is that 
thè GC goes as much straight as thè angle changes slowly. Such a naive observation 
hides thè complexity of thè relation between thè mean square displacement in 2D 
and thè one regarding thè angle variable. 

Hence thè coordinates System used is thè polar System with thè origin centred on 
thè soma, where thè RW starts, and thè angle initial condition set to zero: 

r n - r n -1 = x ■ pcos 9 n + y • psin 9 n (4.1) 

The step is generally not hxed, but thè speed is. Generally speaking to compute thè 
evolution only thè stochastic variable 9 is drawn. 
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4.1.2 Geometrical requirements for thè RW algorithm 

Outgrowing GCs are constrained to go straight-forward frorn thè rigidity of micro- 
tubules. This simple argument exclude thè possibility of drawing thè angle from an 
isotropie distributions, as thè uniform distribution over 2tt. The next step direction 
has to be influenced by thè previous in thè angles variable space, thè process, where 
coordinates have a temporal correlation, is also known as Brownian motion. We 
can set whichever functional forni to determine thè walker stochastic evolution 
but we started from thè Gaussian. Such idea relays on some generai results from 
cultured neurons’ experiments [Renault, 2015| ; and a micro-physic explanation of 
such behaviour is based on simple statistical mechanic System fluctuating around 
thè minimum energy. Let’s consider thè micro-tubules stiffness k m t and thè energy 
required to bend thè neurite with a certain angle, looking at thè steering as a 
fluctuation from thè equilibrium: 

rAO 

A E= / k(6 — 6o)dd = kA9 2 (4.2) 

Je o 

A fl 2 

P(A6) ~ ~ exp(-- ~/3 e ff ) (4-3) 

Where /5 e // collects all thè System Chemical and thermal details and is inter- 
preted as thè variance a. Other dynamics, with a more robust theoretical back¬ 
ground, as those leading to Levy-flights or power-law distributions are presented 
in |Metzler and Klafter, 20041. To study thè Gaussian model and check its results 
is an useful work, at least to show its limits. 

Let’s recali briefly thè growth cone we have described so far, thè factors we want to 
take into account are thè following: 

• thè stiffness of micro-tubules Neurite’s microtubule core prevents thè 
System from abrupt turning. [Martin, 20131. 

• thè soma-tropism of neurite: Neuron’s biological objective is to reach 
out other neurons. The neuronal targets are reached by thè growth cones 
which undergo self-referential forces preventing from self-rolling and back- 
turning. [Memelli et al., 2013 

• Space isotropy: The growth cones haven’t any spatial preference in an 
isotropie environment and everything is symmetric respect thè extrusion angle. 

• Ballistic assumption: The measures performed on reai neurons show ballistic 
behaviour, thè GC direets straight-forward away from thè soma. 

Furthermore we have to consider thè requirements set by thè simulator itself: 

• The computational cost must not increase with thè distance from thè soma. 

• The algorithm has to be expanded to an environment interaction model. 
Furthermore we define a physical property, thè persistence length l p . 
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Correlation and Persistence Length The correlation in 2D can be defined 
through a vectorial picture of thè process: Let’s bi thè element of thè trajectory 7 
which thè segment n belong to. Let’s each step has thè same length b. The 

correlation between thè segments is: 

= = { c °s(A0ìj)) (4.4) 

Assuming thè process in a stationary state, thè correlation function will depend only 
by thè distance s = i — j, whether this is not thè case we set thè site i Constant to 
zero and j stays thè only variable. Hence, developing thè function around Ad ~ 0: 

C(») = (cos(A0,)) = £ VA «A».) 2 ") 

0 V z,1 9 


while A 8 S remains dose to zero, thè expansion can be reduced to thè quadratic term. 


C(s) = l- 1 -((A0 S ) 2 ) 


(4.5) 


In this picture of thè process thè correlation is stili an abstract quantity, it does 
depends by thè time interval and misses a characteristic length. It is way more useful 
to define a characteristic length that can be obtained measuring mature neurons, an 
experimental measure of microtubules persistence length is showed in Fig. 4.1 The 


persistence length should be invariant on coarse graining of thè trajectory, while thè 
correlation between two segments distant s time intervals changes when thè time 
step is augmented or reduced. 

Such characteristic length is a generai property of locai processes in physics, while it 
disappears in phase transition and long range correlations: thè functional form of 
correlation depends by thè recursive equation, it is exponential decaying for short 
range phenomena and power law for long range. The algorithms we will use belongs 
to thè first class and thereafter we can define thè persistence length, regarding to 
reai neurons this cannot be said, because some process, i.e. self-avoiding walks, lacks 
a characteristic length and thè correlation decay slower. Whether thè reai neurons 
should belong to thè latter category different algorithms will be implemented. To 
characterize thè neurite properties thè persistence length is a precise tool. It is 
expected to be a macroscopic property of thè System, thè average length of a straight 
segment originating in thè soma. The persistence length depends by thè parameters 
of thè problem and is a property of thè physical System, which has a well defined 
space-scale. The persistence length can be obtained placing thè function C(i) in 
thè reai space. Assuming thè System varies continuously with a function x(t), we 
have some points of thè curve, thè resolution doesn’t affect thè problem since thè 
scale is fixed. The variable s is thè length of curvilinear path, that is thè length of 
thè dendrite itself. The variable s is measured in firn and thè persistence length is 
defined as :l p = min{s : C(s ) < W 1 } 

C(l) = exp(—^) 

ip 


(4.6) 
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Figure 4.1. Persistence length of dissociated microtubules.The persistence lengtli 
can be calculated from mature cultured neurons in vitro. Measures thè persistence length 
of dissociated microtubules can lrelp to forecast thè relation between thè neurite diameter 
and its stiffness. This article is particularly relevant because implement a method to 
measure thè persistence length for hxed microtubules. The pictures report results from 
thè experiment documented in Martin, 2013| . (A) TIRF image of microtubules sparsely 
decorated with fluorescent GTP analogs. Scale bar is 5 mm. (B) Microtubule trajectories 
from (A). (C) Average of cos(A6b.) as a function of path length, s. The data shown are 
an average across all trajectories in (B). Fewer independent data for long trajectory 
lengths lead to statistical fluctuations. The persistence length for these microtubules is 
150 pm 
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4.1.3 Choosing thè appropriate algorithm, pros and cons of differ- 
ent descriptions 

Bearing in mind thè properties outlined in thè previous section we will describe some 
algorithms and analyse their weakness and strengths. A detailed characterization 
for thè algorithm is offered in Appendix [A] here it is limited to maintain thè unity 
of thè manuscript. There are two main classes: Those implementing a stochastic 
equation for thè angle variable; and those presenting an unified dynamics, with thè 
elongation speed infìuencing thè directional dynamics directly. Some algorithms are 
compared in ??. The fìrst algorithm on thè stage is thè Correlated Random Walk, 
CRW. The dynamics is simple and is governed only by thè following equation, r 
represents a correlated Gaussian variable, with unity variance and zero mean: 


r n = P ■ r n ~ i + \Jl - fi 2 • £ n —1 
0 n +1 = d 0 + (r n )y/a 
ro = Co 


(4.7) 

(4.8) 

(4.9) 


Where Cn is a Gaussian distribution with mean zero and variance one. The angle do 
is set and thè neurite produce a river like path around its direction. The persistence 
length of CRW is set by thè parameter /3 for thè angle variable, while thè neurite 
does not experience turns. This algorithm produces a naive behaviour, thè GC 
outgrowth and keeps going straight for all thè growth process, whether there is a 
resemblance for some neurons’ sarnple it is not usable in thè NetGrowth simulator. 
Indeed it depends by thè initial conditions, whether this is fine for free space neurons, 
those constrained by thè environment will present bizarre behaviours. 

A plausible correction for thè CRW is to set dynamically thè angle do, avoiding thè 
cruciai problem described above. The new algorithm presents a short terrn mernory 
that variates thè direction, thè equation for thè Memory CRW follows, while a 
descriptive picture is offered in Fig. |4.2| 


V n rv r ' 

2 a=o a 


n n ~ i 

2_a=o Q 

r n = P ■ r n -1 + \Jl - fi 2 • Cn-i 


dn +1 — d n T 


(4.10) 


Initialized with: 


ro = Co 

9 0 = d 0 


The term 4.1.3 averages over thè trajectory: takes into account previous directions 
and allows thè random walker to variate its originai direction. Memory acts as a 
rotation and doesn’t take part to thè stochastic computation. It is parametrized with 
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a and dominates thè persistence length. Parameters a, f3 and a cannot be chosen in 
thè whole space or pathologic behaviours will appear: i.e. thè correlated Gaussian 
cannot be larger of memory constraint or thè walker will rotate on itself. They are 
not orthogonal, as it will be shown later, and it’s not possible to map one-to-one 
to stiffness, soma-tropism, etc. This model fìts better for thè NetGrowth scopes 
and it has been implemented in thè simulator, whether it cannot be related easily 
to thè biology of thè microtubules its worst drawback is another. The algorithm 
presented leads to very complex discrete differential equations, which cannot be 
solved (thè detailed computation is in thè Appendix|A|. Hence thè persistence length 
cannot be obtained analytically by thè parameters. This won’t be a problem if a 
numerical estimation could be performed, but thè algorithm produces 2D processes 
whose persistence length changes in a super-exponential fashion in respect to model 
parameter; it becomes very hard to fit via thè three different parameters, as illustrated 
in Fig. 4.4 This annoying properties, added to its non linear dependence by thè 


resolution of time step (indeed for higher resolution thè number of turns is reduced 
abruptly) which cannot be correct with multiplicative factors, makes thè MCRW 
very hard to pass thè "reproducibility" test. Whether it can produce interesting 


shapes as showed in Fig. 4.3, it is not suited to predict neuronal network from 
biological data. The impossibility to produce sensible behaviours with such sirnple 
algorithms constrained us to rethink thè whole stochastic process. The first interest 
remains to produce a persistence length l p that can be set by thè user before thè 
simulatimi, maybe from data obtained in neurobiology laboratories. The complexity 
showed by thè MCRW requires to make a step backwards toward a more predictable 
algorithm: The simplest case we can imagine is thè Diffusion RW, as thè name 
suggest thè angle variable changes with a sirnple diffusive stochastic equation, and is 
governed only by thè variance parameter a: 


@ n +1 — $0 + 


(4.11) 


In this case thè angle shows a completely diffusive behaviour and, as it is showed 
in thè appendix, it leads to a doublé scale dynamics for thè GC in thè xy piane. 
For a scale inferior or comparable to thè inverse of thè Diffusion coefficient thè GC 
presents ballistic pace, for larger scales thè walk is diffusive. Assumed thè elongation 
rate Constant, v, a Constant turning rate v can be defined as a funtion of thè time 
resolution Ai: 


v — —~r— (4.12) 

v A t v ' 

v it is thè number of rotation done by thè GC in l^tm. The diffusion process for 

thè angle impose a linear equation for thè Mean Square Displacemennt of thè angle 

variable: 


(A 0 2 {l)) = avi 


(4.13) 
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Figure 4.2. Memory constrained algorithm. The algorithm sunis previous versors 
with a dumping factor a. The growth cone rotates and aligns itself with thè angle 
weighted over thè trajectory. The fìuctuation due to normal distribution are attenuated 
and thè correlation remains high despite thè process has a large variance a 2 . In thè 
figure thè red semicircle and arrow represent thè direction of last segment, thè green 
circle is obtained with thè memory algorithm instead. 



Figure 4.3. Neurons with different persistence length. In thè upside plots a set of 4 
neurons is growth with different parameters for thè random walk, thè branching model is 
Van Pelt and is discussed further in next sections. The persistence length decreases from 
1 to 4. The adimensionai parameters for Gaussian correlation and memory coefficient 
are set equal for simplicity. Plot assumes respectively values 0.7, 0.6, 0.5, 0.4. 


Persistence vs Memory in MCRW 



Figure 4.4. Persistence length is super exponential for MCRW The picture shows 
thè dependence of thè measured persistence length for thè MCRW algorithm, thè presence 
of memory makes thè algorithm much harder to approach, even numerically. In log-scale 
for thè ordinate, thè persistence length appears super-exponential over a. Different 
colors represent different values of <r 
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A Contraction rate 




B xy position MSD: (A.Y 2 ) 




Figure 4.5. Correlation and contraction measures for Diffusive walkers D coe ff = 
0.025, v = 1, v = 1/i/s Some statistical measures applied to thè diffusive process reveal 
thè simplicity of this model. The four plots show thè presence of 2 scales: The fìrst lasts 
about 100 /im while thè second dominate for longer distances. A The contraction rate 
converges in thè diffusive regime B The xy MSD changes from quadratic to linear and 
C thè correlation between segment direction goes to zero when diffusive regime starts, 
thè persistence length can be obtained by these plots. Eventually thè angle D diffuses 
with rate D coe f f. Detailed explanation of methods in thè Appendix |Aj 


From Eq. 4.6 and Eq. 4.5 
length: 


looking for distance l rnuch shorter than thè persistence 


C(l) 
C(l ) 


exp(-Z/Z p ) ~ 1 - - 

ip 

1 - ^(A 0 2 (l)) = l-aul 

21 “2vqc • Af 
•(A d\ì)) = a 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


Eq. |4.15 is nothing else than thè centrai limit theorem. By Eq. |4.16 thè persistence 
length is analytically computable. 


It is possible to verify its consistency measuring thè persistence length with Eq. 4.4 
which confirms thè exponential decay. Fig. 4.5 resumes thè results. The picture 


of thè process in Eq. 4.11 is very clear, but thè drawback of this method appears 
immediately when implemented: in order to have persistence length consistent with 
biological measured ones, thè value of a must be very little, i.e. less then 0.1 rad. 
This becomes a relevant problem when upgrading to environment interaction: thè 
GC interacts with thè walls in a broader angle, around 180°, whether thè walls 
interact as a potential on thè angle distribution it is necessary that thè probability 
of turning in thè wall direction is more than one. This is to say in such model thè 
GC will be blid to thè walls except for a narrow angle. 
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A Measured Persistence Length \/s model parameter 




Figure 4.6. Persistence length as a function of turning rate In this two pictures we 
verify thè linear equation for thè Uniform Turning Rate model |4.18| Focusing on thè 
concept of turning rate shifts thè attention from thè speed and thè time resolution, to 
thè actual spatial realization of thè process: measuring thè turning rate in micrometres 
instead of second is a robust defìnition that is invariant in respect of thè resolution and 
thè time-step; even it is well defìned for neuron images. A shows thè persistence length 
measured for varying parameters cr and A. The persistence length is measured as thè 
ensemble average of ( cos(9 )). The data are produced by thè NetGrowth simulator itself. 
In B thè angular coefficient obtained above is plotted over thè expected value with 
an excellent precision. The plot confìrms that for thè UTR algorithm thè effective 
persistence length can be arbitrarily set by 6 and A. 


Run and tumble To solve this issue we have followed thè method used in Net- 
Morph simulator: thè turns are diluted in time with an exponential distribution. In 
thè fashion of a run an tumble model. 

The growth cone, has a certain probability to turn, depending only by thè length 
of thè step: Ptum(l) = j- This model is called uniform turning rate (UTR) and 
produces straight intervals, with mean A and exponentially distributed. In this case 
thè coefficient v = j depends by thè turning rate strictly. It is necessary to stress 
we are talking of rate in terrns of spatial units. 

Pstraightil ) = j exp(-Z/A) E [l] = X V ar [l] = A 2 

O \ 

l P = ~ (4.18) 

(4.19) 

With thè last equation thè model was largely enhanced: thè persistence length is 
now independent by other quantities, like growth cone speed or resolution time, and 
can be fixed, case by case, by thè parameter A. In Fig. ?? thè linear dependence of 
thè persistence length is shown for variation of both a and A. 
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4.1.4 Methods an Validation 


To study thè algorithm a vast amount of experiments with single growth cones was 
run, they make up a statistical ensemble. Hence thè percentile distributions were 
computed and thè fits too. To study reai neurons we have looked at NeuroMor- 
phcNeuroMorpho database database, we have used rat retinal neurons since any 
other two dimensionai neuron was available and there weren’t SWC files of cultured 
samples. To measure their properties we have decomposed each branch tree into 
single paths and computed ensemble average. 

This approach has some inner lirnits: study stationary properties of non equilibrium 
process is hard, even more if thè dynamics is lost and a 2D static image is thè only 
track of thè process. Furthermore thè axon is just one and was impossible to study 
its statistical properties on thè reduced set we disposed. With thè dendrites we were 
luckier and a complete characterization was performed. 

Unfortunately thè neuron growth in cultures are not on thè database, which concen- 
trates on invivo samples: to study 2D neurons we had to adoperate neuron from thè 
retina of rats. These neurons are slightly different from those present in neurons 
cortical area, used for thè cultures. Anyway a sophisticated tool to reconstruct thè 
neurite from thè SWC file has been done and can be used in future with effective 
data. The principal measure remains thè persistence length that can be measured 
from mature neurons or even in laboratory with microtubule dissociation as evinced 
in Fig. 4.1 Some of thè analysed neurons are reported in Fig. ?? 


4.1.5 Conclusion 

The analysis of thè navigatimi mechanism for thè growth cones, and thè simulation 
of it, has required a deep insight into bi-dimensionai stochastic processes. The 
properties required for thè simulator have been outlined in thè first section and thè 
algorithms have been evaluated over them. Eventually a modified diffusion process 
has resulted enough simple to be analytically solved stili meaning full for thè scopes. 
The benchmark of thè model are cultured neurons, but relevant data were not found 


on thè Web and so a precise verification is postponed. In Fig. 4.9 


some neurons 


built with UTR model varying thè rate parameter are plotted. Further model for 
environment sensing and direction election will be studied in future, a particular 
attention will be reserved to self-avoiding walks. 


4.2 Growth in a heterogeneous environment: sensing 
thè surroundings 

The way thè growth cone sense thè environment biologically is complex and very 
hard to reproduce in a simulation, it’s roughly described in section ?? and more 
accurately in |Roth et al., 2014 . 

We have summarized over thè whole biological process (focal adhesion, pulling, 
etc... ) and considered thè environment interaction as a bimodal process: thè 
obstacles can attract or repel thè GCs. In case of attraction thè probability thè 
active unity will move towards thè obstacle is higher in respect to a growth cone 
indifferent to thè environment, opposite otherway. The growth cone is kept into thè 








4.2 Growth in a heterogeneous environment: sensing thè surroundings 


63 



A Contraction rate 



0 10 20 30 40 

length fim 




D Angle MSD: (A e 2 ) 



Figure 4.7. Characterization of dissociateci neurons A set of neurons from Neuro- 
Morpho database are measured to extract thè characteristic persistence length. The 
images, in SWC format (see Appendix ??) have been processed and reconstructed to 
analyse thè property of thè whole neurite. Ech neurite has been decomposed in shorter 
and longer branches and thè last one have been averaged for thè salient measures. 
The neurons plotted are two of thè analysed ones. Unfortunately NeuroMorpho has 
not 2D reconstructions, hence we had to use retina neurons and not cortical cells or 
interneurons. In thè lower set of pictures we reproduce thè standard measure to calculate 
thè persistence length, this is attested around 50 micrometers. We cannot affimi this 
measure are sound, indeed thè persistence length appears to be much longer and thè 
diffusive regime far. Hence we cannot say whether or not Random Walk regime will 
appear: In picture C thè correlatimi results strong after 40 /ìto and we cannot predicts 
it will go to zero or not. 
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Figure 4.8. Varying parameters for thè UTR model In thè pictures from 1 to 4 a 
neuron is generated with an increasing persistence length. 

The simulator can reproduce a vast amount of neuronal structure, varying thè branching 
rate, thè speed and thè persistence length. The user who wants to simulate a precise 
neuron has to fìt his experiment or fìnd data of thè type of neuron he wants to simulate. 


cultures walls and cannot overpass culture’s shape, thè obstacles are walls in thè 
culture that occlude thè access to growth cones, as black cells in a crossword puzzle. 
The sensing operation is done at each time step, it has to be unexpensive from a 
computational point of view and deactivable when it isn’t required. 

The model we proposed was hrst sketched in R.Renaud [Renault, 2015 , we have 
changed it to fìt with thè computational requirements of NetGrowth and to integrate 
with thè random walk rnodels previously discussed. 

The wall-filopodia interaction is powered by a geometrical engine. The environment 
model presented in this section allows to simulate complex culture and behaviours 
as thè one described in Fig. ??. 


4.2.1 Interaction with environment as probability functional con- 
volution 

The model implemented in NetGrowth to sense thè environment recalls to growth 
cone biology: it checks thè superposition between a hlopodia and thè wall. When 
thè walls are distant enough that thè GC can step towards it within a time step 
thè wall affiniti/ is considered and thè probability of going that way is modihed, 
increased or decreased. In case thè wall is to dose to thè growth cone centrai domain 
thè probability is set to zero and thè GC is forced to go in another direction, see 
Fig. 4.11[ 

The fìrst characteristic this model presents is thè discretization of space, thè direction 
thè growth cone can move to depends on thè angular sensibility of fìlopodias. The 
probability distribution becomes a histogram and environment perception is coarsen. 
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Figure 4.9. Density distribution for different models of GC navigation To choose 
an appropriate model for tire growth cone navigation we have evaluated tire outcome of 
some classimi algorithms and created a new one. In pictures we show thè average path, 
and a sample path in red, of thè outgrowing neurite. Pictures A and B appear similar, 
yet thè former (UTR or run tumble) is resolution invariant and thè latter (MCRW or 
persistent random walk) is not. The MCRW is governed by three parameters, tlrat makes 
very hard to predict its behaviour for a changing time resolution. Picture C reproduces 
thè algorithm proposed by M eme ll i et al., 2013| , although it is demonstrated it produces 
relevant morphologies it has four parameters which are not directly related to a possible 
experimental measure, also thè algorithm was very hard to tackle analytically and to 
predict its time-resolution dependence, furthermore it is not a generalized Random Walk 
and it is not sound to speak about persistence length in this case. Eventually in D a 
simple random walk is illustrated. This process does not produce plausible morphologies: 
it is regulated only by thè diffusive coefficient, if it is very low thè growth cone goes 
straight and ignore thè surroundings, other way it turns and self-crosses repeatedly. 
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Figure 4.10. Complex shapes and environment interaction. The picture portraits a 
culture with 10 neurons simulated with NetGrowth; neurons with positive wall affinity are 
presented in red, while repulsive are green. The green neurons repel thè walls and spend 
most of their lengtlr in away from thè borders, wit is thè opposite for red neurites, as can 
be noticed on thè up left corner of thè picture. The picture shows thè variability that thè 
environment interaction model can produce with a single parameter, thè edge affinity. 
The wall-interaction algorithm is arnong thè completely new features introduced with 
NetGrowth simulator. It is possible to affect thè Growth Cone navigation - changing 
both thè speed, thè persistence and thè elongation threshold - as a function of thè 
substrate it stands on. This feature becomes particularly relevant in respect of thè 
several experiments conducted on mixed substrates. 
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Figure 4.11. Interaction as functional convolution of probability The growth cone 
senses thè environment with its fìlopodia (whose length is tunable) and changes thè 
distribution of thè navigation model: The probability of going in a certain direction, 
which is computed within thè RW or RT model, is convolved with thè potential offered by 
thè environment, be it attractive or repulsive. In thè picture, obtained by |Renault, 2015j 
thè intrinsic probability is represented in blue, while thè environment potential is pink. 
The result convolution is drawn in red. It is a continuous distribution which is reproduced 
in thè simulator. On this purpose we have tested thè model varying thè number of 
discrete angles. The wall on thè upper side acts an attractive effect, stronger than thè 
fasciculation with thè other axon on thè right. 


To simplify this issue thè fìlopodia density is taken uniform and maintained in a 
certain angle ahead thè cone. 

The probability is convolved with thè probability distribution computed in thè 
frame of random walk, in such a way thè stiffness and other properties of growth 
cone spatial dynamics are kept by thè convolving probability distribution, while thè 
environment can pulì thè axons by its side and affect thè overall picture. This model 
has some drawbacks which are introduced and solved in next sections. 

4.2.2 Implementation in NetGrowth 

In order to implement thè model in NetGrowth we needed to have a simple picture 
of thè sensing unit and to fit it inside thè simulator growth cone unit. The easiest 
way is to use a geometry engine and let it compute thè intersection between wall 
and fìlopodias, then create thè distribution histogram, in each fìlopodia direction, 
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and choose stochastically one of them. 

The sensing happens if thè distance between thè GC and thè obstacles is less than thè 
hlopodia length, which is estimated 20 microns |Cheng et al., 2002' , thè potential of 
attraction with thè wall can change depending on thè distance, i.e. be stronger if 
thè wall is distant more than half of hlopodia and even be repulsive in thè proximity 
of thè cone centrai domain. The parameter to define thè number of hlopodia in 
thè GC is angular resolution, their length is identihed with filopodia length and thè 
attraction is regulated by thè edge affiniti /• 

The model interacts with thè random walk model: thè direction is chosen after both 
thè distribution are computed. The limitation introduced is to use a discrete set of 
angles for thè direction model instead of a continuous variable, we have performed 
some tests to verify this drawback is irrelevant, given thè angular discretization is 
enough dense. 


Renaud model for probability election, from continuos to discrete prob- 
ability 


Let’s introduce thè model itself with a step by step description than can be traced, line 
by line, with thè code in |4.1[ The hrst action of thè GC is to sense thè environment 
computing thè intersection by its hlopodias and thè wall. The hlopodia are numbered 
and this number is depends by thè time resolution. In function compute pulì and 
accessibility thè probability for each hlopodia-direction is initialized to one and then 
increased or decreased by thè edge affinity parameter in case of contact, maintained 
otherwise. This step produces an histogram, each direction has a multiplicative 
factor that can be grater (affine) or lower (repulsive) than one. 

In a second step thè histogram is convolved with thè Gaussian distribution of thè 
angle, thè procedure is detailed in Fig. 4.12 At thè end we have a histogram of 
values, doesn’t matter if not normalized, and we have to choose one direction among 
them, with a weighted probability. 


Listing 4.1. Implemented algorithm for environment sensing 

void GrowthCone::compute_pull_and_accessibility( 

std: : vector<double> &directions_weights, mtPtr rnd_engine) 

I 

//set all filopodia to 1 

for (int i = 0; i < filopodia_.size; i++) 

I 

directions_weights.push_back(l.); 

} 

while (all_nan and not stuck_) 

I 

// initialize directions_weights 

// test thè possibility of thè step (set to NaN if position is not 
accessible). 

//thè filopodia are checked within a [-4 sigma, 4 sigma] interval, 
sigma is thè only physical parameter here. 
kernel().space_manager.sense(directions_weights, filopodia_, 
geometry_.position, move_, sigma, move_.module, 1., nan("")); 

//if any direction isn’t available set ‘all_nan‘ variable 
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Figure 4.12. Environment aware directions’ distribution The picture (A) offer a 
schematic but consistent idea of thè NetGrowth model for environment interactions, thè 
GC in green elongates its fìlopodia towards thè exterior, they touch thè surroundings 
and return thè information to thè GC centrai domain. Depending on thè parameters 
chosen thè fìlopodia can be attracted or repelled by thè obstacle. In this case thè three 
centrai fìlopodia (red) are to dose to thè obstacle A and they are repelled, laterals (blue) 
are attracted instead, while thè four fìlopodia on thè right are to dose to thè obstacle B 
so as to set thè probability of going that direction to zero. The environment interaction 
is summarized in picture (B) which multiplied by thè normal histogram (C) produces 
thè convolved histogram (D). The last histogram is thè baseline probability to compute 
thè next step, thè direction is drawn from thè distribution with a reservoir sampling 
method as discussed in its section. 
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all_nan = allnan(directions_weights); 

//this is thè advantage of variable sigma method: whether thè cone 
looks stucked try to enlarge thè sigma and accept broader angles. 
if (all_nan) 

if (abs(move_.sigma_angle - M_PI) < le-6) 

{ 

// completely stuck 
stuck_ = true; 

> 

else 

{ 

// increase sensing angle and reset weights 

move_.sigma_angle = std::min(2 * move_.sigma_angle, M_PI); 

for (int i = 0; i < directions_weights.size(); i++) 

directions_weights[i] =1.; 

> 

> 

> 

// if thè cone isn’t stucked feel thè sourroundings and fili thè 
discrete distribution. 

else 

{ 

// test thè presence of walls to which thè growth cone could 
// be attracted 

kernel().space_manager.sense(directions_weights, 
filopodia_, geometry_.position, move_, 

filopodia_,finger_length, 1., filopodia_.wall_affinity); 

// check stronger interaction if filopodia is wall is dose enough 
if (0.5 * filopodia_.finger_length > move_.module) 

{ 

kernel().space_manager.sense( 

directions_weights, filopodia_, geometry_.position, 
move_, filopodia_.finger_length * 0.5, 1., 
filopodia_,wall_affinity * 2); 

} 

> 

> 

> 

void GrowthCone::compute_intrinsic_direction( 
std:: vector<double> &directions_weights) 

{ 

for (int n = 0; n < filopodia_.size; n++) 

{ 

directions_weights[n] *= filopodia_.normal_weights[n]; 

> 

} 
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Sampling discrete distribution 

Once thè directions’ histogram all we need to do is to choice a direction. To draw 
an element from a weighted list there is a naif and robust method: fìrst we sum 
up thè weight of each element and obtain thè normalization factor, then we draw 
a number from uniform distribution, multiply it for thè norm and, ordering thè 
elements somehow extract thè one whose interval contains thè drawn number; in a 
more figurative way we will assign each element an area proportional to its weight 
and choose a random point on thè piane. With this simple method thè normalization 
isn’t required and this is useful to simplify thè method of variable sigma that we are 
going to describe. 

Geometry with GEOS and variable sigma method 

To compute thè intersection we have used a simple and powerful geometry library 
for C++. GEOf0 allows for python binding, is thread safe in parallelellised codes 
and very simple and efhcient in respect to more famous projects, which often require 
a huge modificatimi to thè code to include them into. The code we have used is very 
easy: we create a bunch of straight lines, as many as thè number of set fìlopodia, 
where one of thè two extremal is thè centre of thè fìlopodia and thè other is thè 
projection of thè angle at thè sensing distance. Once thè fìlopodia are created 
thè engine evaluates whether they intersect or not thè environment, thè latter is a 
geometrical object passed to GEOS at thè simulation kernel initialization. Further 
details will be offered in thè documentation of NetGrowth. The interesting point in 
such a paradigm is thè method we used to test at thè proper angle and convolve thè 
environment distribution with thè intrinsic Gaussian without compute thè Gaussian 
at each step, see Fig. |4.12 


We have introduce a simple and smart method we have called ‘variable sigma 1 . Indeed 
we use this fact: a set of points where thè normal distribution can be evaluated 
X = [xo, ..x n ] are self similar to thè distribution computed over thè sanie set of point 
multiplied for a fixed Constant: X' = [ax o,.. ax n ] that is 


exp(W) = 


exp(X') 


a 


Such a property of thè normal distribution makes very easy and simple thè computa¬ 
tion, indeed thè whole operation can be computed over a distribution with <7 = 1 and 
thè normal weights don’t require to be computed again. The useful of thè method 
comes up when thè GC is stucking in a corner too, after few runs with environment 
distribution set to zero for each fìlopodia thè sigma augments and let thè GC to get 
out of thè corner. 


4.2.3 Turning angles 

A fìrst check of thè quality of thè algorithm can be devised reproducing in silico thè 


experiment discussed in Fig. 2.7 When thew wall affinity is positive thè GC will 


adhere thè wall, it can follow thè edge when it is straight and this will be either in 


2 Geometry Engine Open Source, https://trac.osgeo.org/geos 
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agreement with thè intrinsic to avoid turning, due microtubules’ stiffness. Hence 
when thè wall turns abruptly if it will follow thè edge or keep going straight will 
depend by thè strength of attraction and its counterpart, thè microtubule stiffness, 
which we have parametrized with thè persistence length. In Fig. |4.13 we discuss an 
experiment with NetGrowth with fixed persistence length and variable wall affinity. 
It reveals there is a strong dependence by thè latter factor. Devising an in vitro 
experiment ad hoc will be possible to relate thè parameter to some reai properties 
of thè neurons. 


4.2.4 Mechanical diode 


The second benchmark for thè model is offered by thè mechanical diode devised 
by R.Renaud. Such a device can be very useful in a culture because it selects thè 
direction of axons (and thè connectivity then) and doesn’t require any chemical 
treatment for thè axon guidance, that is, can be used massively for an unsupervised 
culture, thè experiment is discussed in detail si in Fig. |2.9[ The diode is tested 


in Fig. |4.14 and offers quite good results. The study is purely qualitative since 
thè dynamics of thè axons i much more complex than thè one reported. Further 
experiment, indeed, reveal an high correlation between thè narrowness of thè walls 
and thè speed of thè GC. Such evidence lead to a complicate model, where thè 
environment interaction module is interacting with thè elongation module. Further 
studies are required in this direction and for thè moment thè qualitative overview is 
enough. 


4.2.5 Conclusion 

The module for environment interaction has been qualitatevely tested with experi- 
mental data, a correct test would require a dedicated experimental setup. However 
changing thè parameters produce relevant and comprehensive behaviours. The diode 
implemnted at IPGG was tested and offered consistent results. On thè other hand 
experiments at IPGG, and other authors too |Bak and Fraser, 20031, describe thè 
importance of axon fasciculation both in two and three dimensionai cultures. Fascic- 
ulation isn’t yet included in thè simulator since a more sophisticated computation is 
needed for. The group planes to work on a new model to include this feature. 
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Figure 4.13. Growth Cones’ response to turning angles In this picture we present 
some quantitative data that can be matclied with in vitro experiments. The neurons are 
growth with thè simulator in thè confinement illustrated in picture B. Different set of 
parameters have been tested, varying both thè length of thè fìlopodia, i.e. thè distance 
up to thè wall is sensed, and thè affinity with thè wall itself. For each parameter thè 
average angle of thè outgrowing neurite was measured, thè angle is averaged on thè 
spatial area illustrated by thè red stripes in B. It is shown in thè histogram in box A that 
thè affinity affects sensibly thè angle distribution, thè histogram represents thè amount 
of neurites with a certain angle, it is indeed thè distribution of thè angles. The mean 
value is expected to be zero for thè symmetry of thè confinement, hence thè variance was 
calculated and thè values are plotted as function of wall affinity in thè bottoni plot of B. 
It is evident that for an increasing wall affinity thè distribution fìatten, and even thè 
tails are more populated then thè center. Pictures C and D show thè effect of different 
fìlopodia length: it changes thè distribution variance too. When thè fìlopodia is short, 
i.e. 10 /iro, thè wall is not sensed by thè growth cone and it is indifferent to thè wall 
affinity parameter. In thè two picture two different angles of turning are tested. The 
black line shows thè most of growth cones does not follow a sharp turns of 90 °, neither 
for high affinity and longer fìlopodia. 
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Figure 4.14. Mehcanical diode tested with NetGrowth This figure qualitatively 
verifies thè diode implemented at IPGG, reproducing thè result with NetGrowth. In 
thè picture thè orange elements on thè right are diodes up-to-bottom while thè grey 
elements are bottom-to-up. 

In picture A thè effect of a non interacting model (edge affinity is 1.): thè neurons are 
kept inside thè culture but any attraction nor repulsion from thè walls isn’t present. 
We see thè presence of many axons in thè bottoni of thè picture, decorrelated with 
thè presence of thè diodes. In picture B thè GC-wall attraction is switched on (edge 
affinity is 20.): thè correlation with thè diode is very high, most of thè axons in thè 
lower chamber come from thè right (up-to-bottom) diode, as expected. 

The other parameters are identical for both thè pictures: thè persistence length coefficient 
is set to 20. while thè filopodia are long 10 fini, a is 8.2 degrees and thè angular resolution 
is 40. 
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Chapter 5 

From single Growth Cone to 
Neurites 


Neurons connect each other through growth cones, and these connections have to be 
multiple to turn on thè complex behaviour of brain tissue. We have outlined thè 
strategies for this objective are many, but branching is thè rnost widespread for all 
thè neurons. On thè other hand, from thè descriptive point of view branching is 
thè de facto nature of neurites: a biologically relevant simulator must reproduce 
this mechanism as dose to thè reai phenomena as statistic measures can evaluate. 
The branching itself leads to more complex processes that affects thè whole neurite: 
competition, fasciculation, retraction and pruning etc... The wide variety of these 
events and mechanism prevents from an accurate description of thè arborization, but 
it’s possible, as generally is done in physics, to stress thè relevance of some over thè 
others. In a more abstract way we would say that thè best model is thè simplest one 
that is enough complex to describe thè actual phenomena, all further complications 
are avoidable and can cloud thè interesting behaviours. An example of such choices 
is thè following: in thè hrst stage of Netgrowth realization we decided to describe thè 
Actin Waves, an extraordinary phenomenon that shapes thè neurite initiation phase, 
after a while we understood thè minor relevance of this in branching formation, and 
preferred to remove from thè simulation. The behaviour we chose to implement 
are those that can reproduce a complex arborization as thè one showed in Fig. ??. 
First thè Growth Cone is enabled to split, thè split is forced to respect a volume 
conservation law. The neurite can also brandi laterally and originates a new GC. 
Eventually GCs compete each other and while some elongate some other retract. 
When retraction is continuous thè branch is eventually pruned. 

This chapter is divided in three sections: thè branching mechanism, in terms of 
events’ probability is described in thè first section, here we present an adaptation 
of Van Pelt model for arborization, and a uniform branching model for lateral 
branching. The second section introduces thè diameter to compute branching angles 
and branches’ diameters. In thè third section thè branches’ elongation is accounted, 
a model for criticai resource based competition is presented and studied analytically. 
Eventually thè produced neurons are evaluated with standard measures for neuronal 
morphologies and confronted with other academic results. 
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5.1 Modelling branching processes 


Here thè NetGrowth algorithm for branching events are accurately described, fìrst 
thè Growth Cone splitting, and later thè shaft branching. Let’s discuss briefly one of 
thè simulation constraints that infìuenced thè implementation of both thè models. A 
generai request for thè simulator is to reduce thè operation wherever is possible; in 
case of thè branching a naif approach suggests to fix thè branching probability and 
generate a check whether it happens for each step. In such a scenario thè number of 
operation grows with thè number of growth cones and is infìuenced by thè resolution 
time. Another approach consists into switch frorn thè event probability to thè time 
interval distribution. Hereafter to predict thè branching events and apply thè branch 
at thè proper time step. 

The time prediction depends by thè state of thè System, i.e. thè number of growth 
cones, thè age of thè neuron, and so on. With this approach we can reproduce exactly 
thè BEST model |van Pelt et al., 2003a and whichever branching mechanism whose 
interval distribution is known. For thè interstitial branching a similar approach is 
pursued. 


5.1.1 Growth cone splitting, reversing Van Pelt distribution 

To estimate an interval distribution consistent with BEST modeQ introduced in 
Eq. 5.1 we have started frorn thè branching rate: 


Pi{t\n{t)) = Nn~ E 2~ liS D{t) (5.1) 

with thè normalization factor (5-2) 

n(t) 

Af~°° = Y, 27iS (5.3) 


where pi{t\n{t) is thè rate of ith Growth Cone splitting (namely thè transition rate 
frorn a state with nton+1); thè function D{t) is a baseline branching distribution, 
exponential decreasing, determined by experimental data: this distribution states 
for thè diminished probability of observing a branch event in mature cultures, 
frorn [van Pelt et al., 2003a we set D(t) = Aexp(— t/r). The other parameters are 
E which dehne thè competition among growth cones and S which accounts for 
tree asymmetry [^J thè last one is less important for this description because thè 
normalization factor ensure thè branching rate independent by thè parameter S. 

1 One of most cited branching model, it is a fully phenomenological model whose parameters can 
be obtained from statistical measure of mature neurons 

2 For a tree a of degree n: 

n — 1 

A p = ~r Ap ( rj ' ( 5 - 1 2 * 4 ) 
j=i 

is thè mean value of partition asymmetries in thè tree, which, at each of thè n bifurcation points, 
indicate thè relative difference in thè number of terminal segments rj and Sj in thè two sub-trees 
emerging from thè j th bifurcation point. The partition asymmetry A p at a bifurcation is defined as: 

h-s| 

r + s — 2 


A p (r, s) 
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Figure 5.1. Expected number of branching events Varying thè competition parameter 
thè number of growth cone in thè final tree changes, this dependence agrees with thè 
competition model presented lateron. The number of branching events is regulated by 
thè differential equation in Eq. |5.8[ The picture presents thè number of total branches 
expected. 


The average number of tips, i.e. GCs, in thè neurite is governed by thè differential 
equation: 


— = n{t)pi(t\n(t)) 


(5.5) 


We can integrate Eq. 5.1 over all thè possible segments and obtain thè total rate of 
branching p(t\n{t) as: 


n(t) 

P(t\n(t )) = = D(t)n(t)~ E (5.6) 

i 

= ( n(t))p(t\n(t )) = D{t){n) l ~ E {t) (5.7) 

And set thè boundary condition n(0) = 1 which corresponds to start with one 
Growth Cone. The Eq ?? can be solved as a separable first order ODE and we can 
explicit thè (n(t)) as a function of time and branching parameters A,t,E 


(n)(t) = ( 


AE{1 


- p-t /'O 


1 ) 


1 /E 


(5.8) 


( n(t )) is a continuous, reai, variable accounting for thè average number of growth 
cones. Its evolution depends by for varying thè competition parameter is presented 


in Fig. 5.1 


A formai approach we should compute thè probability of having n events, and then 
compute thè average time between two events, given thè parameters (A t)(A,c,E) 
and its distribution. This process it is hard to study correctly, indeed thè probability 
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changes both in time and for thè discrete component n(t) thereafter a Fokker-Planck 
description would require a left-side repulsive barrier evolving in time. Anyway, 
thè interesting distribution is not thè n(t) at thè end of thè growth but thè events 
interval. A simpler but less rigorous approach is to estimate thè amount of time 
required to increase thè number of GCs, given p(t\n(t)), thè branching probabilities 
of thè individuai terminal segments per unit of time. Let’s proceed by step and 
assume D(t ) is Constant in time, that is thè rate of branching is uniform during thè 
growth, this rate, v, depends only by thè number of growth cones and is Constant 
in thè interval between two branching events. In thè limit of a short interval A0 
dehned P(X(t) = n) = P n (t) with X thè number of branching events, thè following 
equations hold: 


P n (t + A t) = Po(t)(l - vXt) (5.9) 

Po(t) = ~vPo{t) 

P 0 (t) = P 0 e~ ut (5.10) 

This distribution is thè widely used in physics, and neuroscience, and it is an 
exponential decaying probability distribution with mean value 4. Up to here thè 
interval distribution can be computed easily and is exact. Let’s see what happens in 
our case study which is explicitly time dependent. The correct equation to solve, for 
thè first event, becomes: 


Pi(t) =-u(t)P 1 (t) (5.11) 

= - Ae~ ct n~ E Pfit) (5.12) 

Pi(t) =Pi(0) • e ~fi v ® dt e - 1 (5.13) 

(5.14) 

which becomes, for P\ (0) = 1: 

Pi(f) = e ^( 1 - e_tc ) (5.15) 


The hrst relevant aspect of thè Eq. 5.15 is thè probability of having only one branch 
does not disappear for large time. This fact can intuitively explained in this way: 
whether thè baseline has a short characteristic time, thè rate decays so fast thè any 
branching event can happen. Therefore thè interval mean time is: 


(A t) = r ±c ) 

J o 


(5.16) 


For second and high orders interval thè computation becomes tricky, indeed thè base¬ 
line is an independent process. The explicit time dependence prevents to compute 
interval distribution with simple methods, a correct computation should comprehend 
thè interval distribution of all previous extracted intervals. In principle, after p n (t ) 
is obtained thè first return theory (^Parisi, 2010) can be used, hence first return 


probability f(t) computed and eventually thè average interval and its distribution. 
In order to make thè computation feasible, and avoid a demanding analysis, thè 


3 A t 'K<1 
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time dependence can coarse grained and substituted with average values: 

The methods can be outlined this way: First, thè rate dehned in Eq. |5.6[ is approxi- 
mated by its value at thè beginning of thè interval, 

D(t) ->• D(t 0 ) = exp(—t 0 /r) 

, this approximation is valid if thè rnean interval obtained this way is much shorter 
then thè characteristic time. Second thè problem is not solved in generai but 
dynamically: we are, in principle, interested to thè interval distribution of thè 
average process, but, to forecast thè next event with thè proper distribution is 
enough for thè simulator. Hence thè number of growth cones is assumed as a 
parameter of thè process and thè interval distribution, given n events happened, 
becomes: 


v(t n ,n) =D(t = t n )n E 

(5.17) 

ro o 

(A t(n)) = / texp(tu(t n ,n))dt 

J 0 

=A~ l e tn/r n E 

(5.18) 


where t n states thè time of thè nth event. The constraint introduced above is 
respected for: 


r > (A(5.19) 

ln(r) > — + Eìn(n) + ln(A) ~ — (5.20) 

r r 

(5.21) 


This distribution was implemented in NetGrowth, thè results are confronted with 
thè BEST model in Fig. 5.2 After thè interval time is extracted thè exact branch 
where thè split happens is selected by a weighted choice basic algorithm. The weight 
depends dynamically by thè branch order in thè tree and by thè symmetry parameter 
S introduced above. 


5.1.2 Interstitial branching 

For thè lateral branching a similar approach is proposed following thè results 
in |Szebenyi et al., 1998 . Each segment has a certain probability of generate a new 
growth cone, which depends only by a ocal excess of actin. The new branch appears 
laterally and starts a new growth cone. The GC is considered equivalent as it’s 
appeared after a GC split event. For a correct characterization and simulation of 
thè process we investigate thè explicit time dependence fo lateral branching events, 
and thè spatial distribution of new GCs over thè branch length. The cited paper 
describes thè process accurately and reports numerical results in terms of branching 
rate and distance frorn thè leading growth cone, some relevant data are represented 
in Fig. ??. By thè success reported in previous section thè sanie rnethod is applied: 
thè rate is assumed to decay with a certain function, D(t) = t°\ and thè expected 
interval is drawn after thè last branching event. The spot where thè branch rises is 
chosen by an exponential decaying distribution, as reported in |Sz ebenyi et ah, 1998 
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Figure 5.2. Branching times for BEST model Pictures 1 and 2 describe thè statistic 
properties of topological tree generated by thè BEST model, both thè histograms repro¬ 
duce a simulatimi whose number of branches is dynamically computed, thè distributions 
are averaged with their ensembles. This test suite was created to verify whether thè 
equation approximated above reproduce thè results obtained by thè BEST model. 

The branching events are generated in two ways: by their rate, drawing a stochastic 


number for each time step (Eq. 5.6 1 as prescribed from thè model, or extracted from 


thè interval distribution (Eq. 5.18l which is illustrated above. The parameters are equal 
for both thè process: r = 1000, A = 1. The symmetry parameter is not considered 
because thè clioice of thè branching cone is marginai in respect to thè time interval. The 
matching of these two methods is tested varying thè competition parameter: E = 0.3 in 
1 and E = 0.7 in 2 

The fìxed time approximation proposed above fìts adequately with thè BEST model 
distribution, pictures A. The Inter Branching Interval distribution (IBI) in B and 
average time of branching occurrences are well reproduced too. The approximation is 
weak for large interval, as explained in Eq. |5.21[ and this results appear clearly in thè 
third histogram. 
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5.1.3 Conclusions 

The branching models for Growth Cone splitting and lateral branching have been 
adapted for minimize thè computation costs. The exact estimation of thè branching 
rates as functions of parameters and nurnber of cones were hard to resolve exactly 
for thè explicit time dependence, conversely a computational method was introduce 
with great success in reproducing thè originai model. The whole neurite will be 
compared with reai neurites in thè end of thè chapter, here we limit to comment 
some trees in Fig. ?? 
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Figure 5.3. Dendritic trees produced by branching models in NetGrowth The 

pictures represent some of possible dendritic trees realized with thè models described 
above. A and B are arborization produce with Van Pelt model varying thè symmetry 
parameter and thè competition between cones. B, with less GC has E = 0.7 while 
E = 0.3 for A. 

In C thè arborization is produced with lateral branching model. 
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5.2 Competition for criticai resource 


In chapter |3.4.1| we have introduced a simple model for criticai resources competition 
among growth cones. This model is based on thè hypothesis that GC might require 
a specific molecule to grow, which would be available through a given amount 
throughout thè neurite. In a complex neurite thè GCs compete for this resource 
and their dynamics (elongation, retraction, stalling or splitting) depends on thè 
amount of resource that is available to them. In thè following chapter some results 
are obtained for a simplified model. Our model avoids thè integration of thè diffusive 
equation over thè whole neurite; conversely, it follows thè solution of a System of 
SDE, one for each GC. I have to underline this model is thè last of a vast amount of 
competition rnodels, defined and studied during rny period in Paris. I will describe 
thè previous attempts in thè Appendix ?? to avoid an useless complication of thè 
manuscript. We will proceed introducing thè ideas behind and then relating them 
to thè equation. The stochastic differential equations are managed by thè formalism 
used in |Boffetta Guido, 2012 |e Parisi, 2010 


5.2.1 Strive to grow: from ideas to formai models 


The ideas shaping this model are borrowed by (Ooyen, 2004 . The fundamental 
principle is thè amount ot thè criticai resource received is necessary to synthesise 
microtubules or to keep them together |Hely, 2001]. Hence thè amount of criticai 
resource sets thè elongation rate and, possibly, thè branching rate too. 

The first models tested were equivalent to thè stochastic distribution of a quantity 
over thè leaves of a tree. In this frame we had introduced thè geometrical factor 
Q, it reduces thè received amount for high order branches. The received resource 
cannot be completely consumed because some other molecules are necessary: these 
elements can be produced in thè GC itself with a certain rate u. Same how thè 
resource corrupts with a rate rp The consumption rate should not be fìxed but 
depend by externally triggered or stochastic variations, x- Thus thè distribution 
follows thè consumption; this picture of thè process motivates thè second term in 
r.h.s of Eq. 5.22| Eventually thè criticai resource produced in thè soma can fluctuate 
with its own timescale, imposing thè second equation in Eq ]5.22| Starting from this 
sketch, for a growth cone i, we get: 


/ 


CLi — di 


U + 


1 

n 


\ 


(1 + Xi) + 


A (iKdi 

T~d X) Cj Kdj 


(5.22) 


À = Am A - - + C 

ta r d 


1 


— ~(Am — A) + C 


where £ G Af(0,crg,p) and Xì £ A/"(0, ai, pi) are potentially correlated Gaussian 
random variables. In these equations: 


di is thè quantity of resource available, 
u is thè consumption rate, 
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r i is thè leak timescale, 

Ci is thè weight factor (geometrical) of thè GC, 

A is thè amount of resource available inside thè neurite and which can be delivered 
to thè GCs, 

r is such that r _1 = r^ 1 + r^ 1 

A m = t Am 

Once thè quantity of criticai resource is obtained, we can define thè dynamics of thè 
growth cone’s speed v based on thè thresholds: 


<0 if diti < 0 rs 

0 if 0 rs < ciiU < 0 se (5.23) 

>0 if Ose < UiU 

where: 

v r is thè maximum retraction speed when a GC gets 0 resource 

v e is thè maximum elongation speed when a GC gets all thè resource generated by 
thè soma 

Eventually we also define thè threshold 0s, above which thè GC can split. 


0 rs below which thè GC retracts 
0 se above which thè GC elongates 
such that: 

{ CLjU @rs n i 
9rs Vr 

CL%U Ose ny 
Ali-Ose Ue 


5.2.2 Solutions for A 


Before implementing thè models we wanted to forecast its average behaviour and 
thè dynamics that can arise by it. Some of thè used mathematical tools are reported 
in thè Appendix [C] The stationary distribution for thè amount of criticai resource 
produced is: 


f(A) = C exp 


(A-A m ) 2 

2a ì 


(5.24) 


And thè solution of thè average variable (A) is given by: 


(A)(t) — Am + (Aq — A«)e t ^ T 


(5.25) 


Which has been represented in Fig. 
threshold. 


5. I 


in respect to elongation and retraction 
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Density ( A=10.0) s=1.0) u=5.0) tau=2.0) 



A 


Density evolution ( s=1.0; u=5.0; tau=2.0) 



A 

Figure 5.4. Single Growth Cone with criticai resource picture describes thè station- 
ary probability distribution for a Growth Cone in thè criticai resource model, thè vertical 
lines indicate thè thresholds for retracting or elongate. The models parameter can set 
thè mean above or below thè threshold, intuitively we expect for an increasing number 
of Growth Cones thè mean rate of received resource falls below thè elongation threshold, 
further studies should be performed in this direction, studying thè distribution for a 
fixed number of growth cones, and eventually solve thè model for an increasing number 
of Growth Cones. 
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5.2.3 Solution for one growth cone (white noise) 

Let’s move our attention to thè dynamics of thè growth cones, assumed thè equa- 
tions 5.22 hold for a System with N growth cones thè case for N = 1 follows: 


à = —an (1 + x) H- 

Td 


À — — {Am — A) + £ 

T 


(5.26) 


This lead to thè average and squared average: 


A 


(A a) = —n{a)At H-A t 

1~d 

((A a) 2 ) = a 2 n 2 cr 2 At + o(At) 


(5.27) 

(5.28) 


which differs from a Ornstein-Uhlenbeck process for thè presence of a as a factor of 
thè stochastic component. Considering thè average equation for A and assuming (if 
k / t _ 1 ) we can expand 


5.26 


(<*>W = — + Woe"* 1 + , T 

KT d T d (TK ~ 1) 


(Aq - A m )e-^ 


or to 


<o)(t) = — + (a) 0 e~ Kt + A ° Au te- Kt 

KT d Tfj 


(5.29) 


(5.30) 


if they are (we do not consider this case since we will consider k t a 1 , r d 1 ). 


Conditional Fokker-Planck 

Let g A (a,t) = f(a,t\A) be thè conditional distribution for a at fixed A. 

We get: 

dtGA = (cr 2 - n)f + (^2aa 2 + y - naj d a gA + d 2 a g A (5.31) 

(aKa x ) 2 


= d a 


= d a 


— - na) gA + d a 
T d 


+ an(Kcr 2 - 1)^ g A + 


-9 A 


( aKa x ) 2 o 7 ( 

-x- o a g A 


with, since a > 0, 


Vi, d a gA(0,t) = 0. 


(5.32) 

(5.33) 

(5.34) 


5.2.4 Competition between two GCs 

Let there be two GCs, denoted be thè subscripts 1 and 2, thè evolution of their 
available resource in time is given by: 

A ai 


ài = -ftai(l + Xi) + 


Ò2 = -rea 2 (l + X 2 ) + 


Td ai + a 2 

A a 2 
T d ai + a 2 


(5.35) 
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where Xì £ AA(1 ,ai,pi). This System can also be describe by thè symmetric and 
antisymmetric components: 
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- K ai M 

i 

) /3 - K(aiXi ~ C 12 X 2 ) 


(5.36) 


with a = ai + a 2 and (3 = ai — 02 - The equation is integrated numerically and thè 
results are presented in Fig. ??. The method chosen for thè integration is thè Heun 
method, also known as Improved Euler Method, thè method is recali in Appendix |IV| 
The algorithm has been tested in a test suite, here we have verified thè results were 
consistent, i.e. thè growth cones where advancing or retracting by thè amount of 
received criticai resource, and there was actually a competition. Once thè model 
was considere effective we have implemented it to thè simulator, thè weigth Q was 
expanded into two weigths, one topological and one for thè diameter, in thè idea 
thè larger thè neurite, thè most thè received resource. In Fig. |5.5| thè case of two or 
three competing GCs is discussed. 


5.2.5 More than two Growth Cones 


In case thè neurite starts to brandi repeatedly, which is normally thè case, thè result 
becomes soon unpredictable, thè space of parameters of thè model is broad, and 
whether we hope they can be biologically related with a research project focused on 
it, for thè moment we have to guess thè right set of parameters to produce relevant 
and plausible neurons. In thè following pages three case will be evaluated: firstly thè 
effect of an unbalance between thè outgrowing cones, this inhomogeneity is driven 
by thè neurite diameter in this case, Fig. 5.6 Second thè relevance on thè final 


morphology by thè amplitude of stalling regime, among thè retracting and elongating 


thresholds Fig. 5.7 Eventually thè morphology obtained changing thè correlatimi of 
thè stochastic variable £, in case of a large variance, is reported in Fig ??. 


5.2.6 Conclusion 

The elongation rate of Growth Cone has been approached with a competition model, 
gathering some previously discussed methods |van Pelt et al., 2003b| , ||Hely, 2 001 . 
After an intense effort of modelling, discussed in Appendix ??, a System of two 
SDE has been devised to solve thè problem. The differential equations governing 
thè model have been analytically resolved, where possible, and integrated with thè 
Heun method (Appendix |IV| ) . With this model thè neurite is simulated as a complex 
structure, where thè relative outgrowth of some components affect negatively thè 
elongation of other growth cones. Some experiments can be devised to test whether 
this approach is meaningful and relevant, these experiment can be arranged with 
neurons deprived of tubulin or whose MAP2 protein is inhibited. In this case thè 
rate of resource produced and diffused can be related to thè parameters of thè model 
and test which conditions can be reproduced with NetGrowth and which not. Some 
relevant comparation are presented in previous iniages. 
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Figure 5.5. Base study for competition model Here it is reported thè separated study 
performed for thè competition model. The equations in |5.22| are implemented with thè 
Heun method in for two and three growth cones. This helped to study thè overall effect 
of competition on outgrowing cones without thè distracting presence of neurites. The 
picture A reports thè evolution of criticai resource access during one simulation, between 
two growth cones, thè colour proxies for thè time. Pictures B indicate thè distance from 
thè soma for thè growth cones as a function of time, thè position changes during thè 
time and gcs elongate or rtract depending by thè criticai resource received. The pictures 
C shows thè amount of CR received, obviously it is symmetrical for two growth cones 
and less trivial for three. 
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Figure 5.6. Competition enhanced by neurite diameter In thè picture two neuron 
growth with thè criticai resource model are presented. The neurons are produced with 
same parameters except for thè weight Q which triggers thè competition by diameter. 
In thè neuron on thè left all thè brazes receive, approximately, thè same amount of 
criticai resource, which is consumed independently by thè neurite size. In thè neuron 
on thè left, instead, thè fìux of substance received depends strongly by thè diameter, 
this dependence prevents shorter branches to elongate and promote larger ones. This 
mechanism can be enhanced more with a topological dependent weight, that forces thè 
tree to grow symmetrically or not. The weight becomes more relevant when thè pruning 
process is taken into account, in this case thè weaker branches not only stay behind but 
get pruned on thè long run. 




Figure 5.7. Competition with low or high elongation threshold This picture 
presents thè dendrograms of two neurons growth with different parameters in respect to 
thè elongation and retraction threslrolds 9 e i, 6 re , botlr thè case are diameter independent. 
The dendrogram on thè left A present some pausing growth cones, i.e. thè piate branches, 
while thè one on thè right B has only outgrowing branches. The latter dendrogram 
presents, indeed, a greater value for thè elongation threshold, so thè stalling regime is 
reduced and thè inactive growth cones get pruned, leaving more substances to longer 
growth cones. This confìguration can be particularly relevant for thè axon, which can 
choose thè direction following an evanescent gradient of chemiotaxic field and set all thè 
neurite effort to reach thè target. 
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In thè next chapter an overall estimatimi of thè NetGrowth simulateci neurons is 
offered, highlighting strengths and flaws of thè project. 



5.2 Competition for criticai resource 


91 



Figure 5.8. External stimuli with different timescales This example shows thè effect 
of randomness when related to thè consumed criticai resource. The idea is: there are 
fìuctuations in thè employment, or in thè leakage of CR. these fluctuations are due to 
internai or external factors and have a certain intensity and time scale. The intensity can 
he set by thè variance of thè stochastic variable £ while thè time scale can be regulated 
with thè correlation factor. In thè example two neurons -growth with same parameters 
and different correlation times- present characteristic behaviours: in thè first case A thè 
neuron is quite symmetric on all its branches, indeed thè stochastic variable is quickly 
averaged about its mean value (i.e. zero) and neither retraction or elongation is enhanced. 
In thè second case B some branches are stretched way more than their siblings, we 
expect also more pruned branches in this case. 
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simulator 
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5.3 Validating Netgrowth morphologies 

The validation of NetGrowth simulator is a fundamental part of thè project: thè 
morphologies and networks generated should be a dose reproduction of those observed 
in neurobiology laboratories. Some test cases on which NetGrowth could be used 
are thè following: 

• How does thè inter population connectivity change when thè aver- 
age branching rate is increased? To aid neurobiologist in relating thè 
macro-scale properties of single neurons and neuronal networks on Petri dish, 
i.e. connectivity, graph closeness measures, arbor extension, etc..., to thè 
micro-scale processes and parameters: The modularity of thè simulator allow 
for thè elaboration and implementation of new models in a parallelized simula¬ 
tor, and we hope once thè simulator will be ready to use some labs will pay 
attention to its possibilities. 

• What is thè resulting network of a certain culture structure? How 
do topologies affect synaptic plasticity? To furnish plausible networks 
and morphologies to theoretical researchers fuels thè investigation of more 
plausible network topologies, which, integrated into activity simulator, e.g. 
NEST, allow for thè simulation of cultured neurons dynamics. The verisimili- 
tude is a necessary requirements to keep thè investigations relevant and prone 
to applications. 

Having stressed thè importance of obtaining realistic neuron morphologies and 
effective environment interactions, then cornes thè hardest and trickiest part of thè 
whole project: testing thè verisimilitude of thè generated neurons and networks. 
As generally done in Science thè dividi et impera approach should make thè case, 
indeed thè properties to test are many and verifying thè branching rate together 
with thè elongation model could not unveil thè insufficiencies of one or thè other 
formai descriptions. The simulator should be conhrmed by appositely conceived 
experiments, as those used to benchmark thè environment interactions. These 
experiments require time and competence, and were not performed during thè short 
tirne of rny master thesis. 

On thè other hand one of severe lessons I learnt during these project is what 
happens inside a celi is never obvious: The chemo-electrical processes going on 
beyond cellular membrane are tightly connected, and to study them separately is not 
an easy achievement. It is possible, for this purpose, to inhibit some mechanisms, 
as thè branching, with drugs or genetic manipulation, but it requires an entire 
research project and often offers only qualitative results. Furthermore thè behaviour 
of thè single neuron depends thoroughly by thè environment, i.e. thè rest of thè 
culture. Neurons are not solitary, either their lifespan is strongly correlated to their 
environment and to thè presence of other cells. Glia, for instance, are essential 
to feed thè neurons and form synaptic connections, same how thè neuron will not 
survive if they don’t establish synaptic connections, and when they do they prune 
large part of their body which is not functional to its scope: receive and transmit 
electrochemical signals. With these premises to test thè NetGrowth simulator would 
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require a long session of experiments, thè parameters should be fine tuned with 
observation and accurately related to live cells behaviour; different environment 
should be testes and see whether they affect or not thè growth and how. This was 
not possible in thè context of my thesis, therefore only a short hint of validation will 
be presented in this manuscript. 

The validation will be divided in two sections, in thè first I will review two 
articles which offer a consistent blueprint of measures, towards an integrated and 
effective statistic of neuron morphology. In thè second I will produce these measures 
for some morphologies generated with NetGrowth. 


5.3.1 Measuring neuron’s morphologies 


The articles we are referring for morphologies validation are |Costa et al., 2010 


Snider et al., 2010 . The authors of these review concentrated on generai properties 


of neurons’ morphology, stressing thè relevant aspects of neurites, Fig |5.9| and thè 
consequent measures can be performed to classify these neurons and to validate in 
silico generated neurons. These approach unveils thè most of neurons in a reduced 
area of thè 2D space, which is very reduced in respect to thè neuron that can 
be generated in silico both with generative and statistical procedures. Anyway a 
relevant properties of neuronal morphology is thè spatial density function, a function 
that specifies thè probability of finding a brandi of a particular arbor at each point 
of thè reai space. From thè analysis of over a thousand arbors from many neuron 
types in various species, jSnider et al., 2010 has discovered an unexpected simplicity 
in arbor structure: all of thè arbors examined, both axonal and dendritic, can be 
described by a Gaussian density function truncated at about two standard deviations. 
These result, very surprising indeed, is valid for neuron growth with a completely 
isotropie space and can be very useful to verify thè quality of morphologies produced 
with NetGrowth. In respect to thè previous analysis thè classification is reduced 
to four parameters in 3D and three in 2D: thè total length of dendritic tree and 
thè 3(2) parameters for thè standard deviation in each dimension, an example is 


reported in Fig. 5.10 Another standard measure can is thè Sholl Analysis: a method 
of quantitative analysis commonly used in neuronal studies to characterize thè 
morphological characteristics of an imaged neuron Sholl, 1953 . It was first used to 
describe thè differences in thè visual and motor cortices of cats. While methods for 
estimating number of cells have vastly improved since 1953, thè method Sholl used 
to record thè number of intersections at various distances from thè celi body is stili 
in use and is actually named after Sholl. In order to study thè way in which thè 
number of branches varies with thè distance from soma, it is convenient to use a 
series of concentric spherical shells as co-ordinates of reference these shells have their 
common centre in thè soma. Sholl also realized his method is useful to determine 
where and how big is thè region where synapses are possible, an example of Sholl 


analysis on a reai neuron is offered in Fig. 5.10 


In thè next section thè morphologies produce with NetGrowth in isotropie space 
will be validated with thè method of arbor density. 
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Figure 5.9. The variables of neuromorphology spaceThe neurites are complex struc- 
tures and their geometry and topology can be classified on thè base of 20 parameters, 
picture A. The variables can be global, as thè nurnber of total branches, or locai, (i.e. 
thè distance between two branching points, thè angle between two sibling branches, 
etc.... The author proposes a PCA description, which maps thè many variables in a 2D 
space, picture B which is more easy to visualize and compare, in thè picture thè grey 
spots are thè morphologies produced in in silico while thè black spots are those obtained 
by thè NeuroMorpho database. |Costa et al., 2010 



Figure 5.10. Density of dendritic arbor and Sholl analysis A broad study, on 
tlrousands of reconstructed morphology, shows thè neurite has a density function which 
is Gaussian isotropie space, picture A, and proposes that four parameters are enough to 
describe thè dendritic arbor, thè dendritic arbor length and thè standard deviation in 
each direction. [Snider et al., 2010] A similar analysis, picture B, is used since 1953 to 
classify neuron morphologies, its inventar, which gave thè name to thè method, asserted 
Sholl analysis is valid to estimate thè probability of creating a synaptic bottoni at a 
certain distance from thè soma. [Sholl, 1953] 
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5.3.2 Simulateci and cultured neurons 

This section is indeed very relevant for thè scientifìc validity of thè work I have 
done for thè NetGrowth simulator. Sirice thè simulator will produce morphology 
for confined culture a sound validation of thè morphologies should be done in that 
direction. To studies these morphologies in vitro it is, sadly, very hard. The arising 
complications are many, one over all thè death of thè selected neuron. Furthermore 
thè acquisition of thè irnage into a descriptive format is not trivial, many software 
exist to do it and thè most are proprietary and sold with thè microscope itself. 
Eventually we have to consider thè evaluation of thè morphology cannot be done 
on generai basis, indeed thè morphology depend deeply by thè parameters of thè 
experiment: thè area of thè brain whose cells belong, thè species of thè animai, thè 
conditions of growth and thè time of thè neuron. It is evident that a sound work of 
validation should be done with thè help of a neurobiologist and with thè facility of a 
neurobiology laboratory. These accomplishments are well beyond thè possibilities, 
in terrns of time and knowledge, of this manuscript and are left open. Anyway thè 
measures themself will be performed on thè produced neurons and attached, both to 
prove thè results are similar to those we can find on publications and in thè web in 
generai, and to figure out how thè conclusion of thè NetGrowth project will appear. 
Even to find neurons in 2D in thè web is a hard task, indeed thè NeuroMorpho 
database collects only neurons obtained from in vivo, and it lacks of neurons growth 
on Petri dishes. Since thè following analysis is completely arbitrary I opted to pay 
tribute to Santiago Ramón y Cajal, thè first neuroanatomist.He was fascinated by thè 
complexity of neurite arbor and draw a huge quantity of hand-made reproductions of 
neurons from in vivo, his work was thè basis for dendritic arbor studies for half of a 
century. A set of neurons drawn from in vivo observation are presented in Fig. |5.11[ 
thè morphologies reproduced are thos of Chandelier neuron (Fig. 
celi Fig. ?? and Multipolar celi Fig |5.13 

Chandelier celi like 

5.4 NetGrowth simulateci cultures 

The scope of NetGrwoth is not only to produce valid isolated neurons but to create 
complex networks simulating thè growth of neurons in complex confinements. In 
this section two example of such ability will be offered, thè first is a very classic 
confìguration, with hundreds of neurons contained into a Petri dishes. The second 
example is derived instead from thè set of experiments realized by R. Renaud at 
IPGG, in this case thè neurons will be confined in a two-chambers culture. The 
number of neurons is usually around hundreds of thousands in these cultures, but thè 
simulations were performed on a personal computer, whose limited RAM and CPU 
does not allow for such a huge simulation. The analysis of thè networks obtained is 
beyond thè scope of these thesis: thè scenarios which are opened with this simulation 
are numerous and wide. For instance thè degree distribution can be investigated 
and compared with in vitro measured degree distribution. Hence it is possible to 
compute thè centrality of neurons, and see whether their inhibition slows down thè 
activity of thè whole culture. 


5.12), Purkinje 
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Figure 5.11. Different Types of Neurons drawn by S. Ramón y Cajal A. Purkinje 
celi B. Granule celi C. Motor neuron D. Tripolar neuron E. Pyramidal Celi F. Chandelier 
celi G. Spindle neuron H. Stellate celi, based on reconstructions and drawings by Cajal 
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Figure 5.12. Chandelier neurons or chandelier cells are a subset of GABA-ergic cortical 
interneurons. Chandelier neurons synapse exclusively to thè axon initial segment of 
pyramidal neurons, near thè site where action potential is generated. [6] It is believed 
that they provide inhibitory input to thè pyramidal neurons, but there is data showing 
that in some circumstances thè GABA from chandelier neurons could be excitatory. 
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Figure 5.13. Multipolar neurons A multipolar neuron is a type of neuron that possesses 
a single axon and many dendrites (and dendritic branches), allowing for thè integration 
of a great deal of information from other neurons. These processes are projections from 
thè nerve celi body. Multipolar neurons constitute thè majority of neurons in thè centrai 
nervous System. They include motor neurons and interneurons and are found mostly in 
thè cortex of thè brain, thè spinai cord, and also in thè autonomie ganglia. Total neurite 
length=23195.95 firn Number of terminals=35 
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The very important step is to connect thè neurons each other throughout thè 
synapses; these chemoelectrical buttons are among thè most amazing and complicated 
elements of neuronal dynamics studies: they can be excitatory or inhibitory and 
thè understanding of thè most complex functions of nervous System pass by them. 
The generai mechanism is thè following: thè pre-synaptic neuron releases Chemical 
elements into thè synapse, these are absorbed by thè post-synaptic and activate or 
inhibit thè ion channels, favouring or preventing thè rise of dendritic spikes. Synapses 
are not considered now, in next sections thè neurons will be considered connected 
whether thè intersection among thè dendritic tree and thè axon’s branches is not 
ernpty. 



Postsynaptic neuron ID 


5.4 NetGrowth simulateci cultures 


101 



Presynaptic neuron ID 

Density of neurites Density of synaptic buttons 




Figure 5.14. NetGrowth generateci neurons in Petri dish culture The figure 
presents results from a Petri disli like confìnement with 100 neurons.The neurons are 
growth within thè environment with BEST model branching model and run and tumble 
navigation model. The data produced are saved in SWC format and post processed with 
tools I developed within NetGrowth project. Picturc A reproduces thè adjacency matrix 
obtained from neurite intersection method: whether an axon intersects thè dendrite of 
another neuron a synapse is generated. Pictures A.l and A.2 describes thè density of 
neurites, both axon and dendrites and thè occurred synapses. The correlation is evident, 
this is due to thè triviality of intersection method. The last picture shows thè position 
of thè soma of generated neurons. The further develop of NetGrowth will consider thè 
soma migration. 
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Figure 5.15. NetGrowth generateci neurons in diode-like culture The culture 
implemented by R. Renault in jRenault, 2015| is reproduced with NetGrowth and thè 
relative network is outlined. While thè density of neurite is homogeneous in thè two 
chambers thè effect of axons invading from left to right is showed by thè higlrer density 
of synapses in thè right chamber: thè axons from thè left forms several synaptic buttons 
within thè chamber. This confìguration is clear in picture C: left axons in red invade 
thè other chamber with more intensity of right blue axons, this inlromogeneity is due to 
thè shape of thè communication channels. 


Q 

fi 

o 

Sh 


G 

<D 

G 


u 


& 

cd 

G 

>> 


C/D 

in 

O 

CL. 



Presynaptic neuron ID 


Figure 5.16. NetGrowth generateci neurons in diode-like culture, adjacency 
matrix The adjacency matrix for thè diode like culture shows clearly thè global properties 
of thè network: neurons from thè left, in red with IDs less then 100, have much more 
postsynaptic connection in right (ID greater than 100) chambers then opposite. This 
confìguration reproduces a perceptron and is one of basic step towards a neuronal based 
computing device. 
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Part IV 
Conclusions 
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The NetGrowth simulator has been discusseci in thè manuscript, focusing on 
thè biophysical models implemented and their relation compared to state of thè 
art models. NetGrowth is a morphology generator for 2D neurons in cultures. It 
has been designed to produce realistic neuronal morphologies leading to relevant 
network topologies. It can also be a valid instrument in thè quest for relating thè 
microscopie properties of neurons to thè macroscopic evidences of cultures and their 
computational abilities. The introductive work dealt with thè positioning of this 
research project in thè vast and quickly evolving held of neuroscience. A huge 
bibliographic research has been done to evaluate thè publications and thè models 
that could be useful for thè development of thè simulator. 

A relevant research project The NetGrowth project is thè result of thè collab- 
oration among thè Neurophysics group at MSC Lab, thè group of Catherine Villard 
at IPGG, Paris, thè group of J.Soriano at UVB, Barcelona and E.Moses, Weizmann 
Institute, Rehovot. These research groups are involved at different levels of neuronal 
culture studies and they have showed great interest in thè realization of thè simulator. 
The existence of a simple IT tool to generate realistic networks, combined to thè 
impressive quality of NEST activity simulator, opens interesting scenarios in thè 
field of neuronal cultures. The costs in terms of time and biological material required 
to grow a culture are very huge. Hence thè opportunity to simulate thè experiments 
reduce considerably thè risk of mistakes, NetGrowth can help researchers to choose 
thè relevant shape for thè neuronal culture. Although simulations cannot substitute 
thè experimental results and evidences, it can help into predict thè outeomes of an 
experimental setup and offer a comparison tool. The project has revealed interesting 
also for studies in culture with different substrates. The opportunity offered from 
thè environment interactions are indeed very broad and thè realization of mechanical 
diode, as those proposed in [?], is just one of thè possible implementations. An 
experiment that will use thè NetGrowth simulator as benchmark will be realized in 
next months by Soriano and collaborators. The culture will be devised in two or 
more topographical areas, these areas have different heights and allow upper neurons 
to descend, while offer conhnement to neurons in lower levels. Simulations with 
NetGrowth will help to devise these areas forecasting thè resulting networks. 

Devising and implementing models Most of thè work I have done on thè 
project is thè realization of thè simulator infrastructure and thè devising and 
evaluation of thè models. The models were presented by increasing complexity: at 
fìrst thè navigation is studied then thè mechanical interaction. The second part 
deals with models of branching and competition. These last models were thè more 
challenging from thè analytical point of view: thè simplicity they are presented may 
cover thè intense work of modelling in favor of readability of thè whole project. Part 
of thè analytical computations are reported in thè Appendixes. The navigation 
part describes thè study of random walks and thè related generating algorithms: 
reproducing thè navigation has been challenging since a simple random walk was 
not adequate to reproduce thè neurite outgrowth. Usually in physics we deal with 
random walk that self-cross repeatedly and we look only at thè average distance 
from thè origin. In this case thè GC navigation must both reach thè target and 
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leave a plausible dendritic tree behind itself. Up to now thè neurites intersect 
each other, we are working on fasciculation and self avoiding model to prevent this 
inconvenient. The models have been evaluated in this case also for thè usability 
they offer, and for thè needs of a clear parameter space for thè user. The study 
on mechanical interactions instead had presented an intense work of algorithmic 
optimization, switching from a biophysical simple idea to its implementation has 
revealed tricky with respect to effìciency. The model implemented in thè end was 
efficient and reproduced thè experimental evidences. Regarding thè branching I had 
to study and solve non trivial stochastic problems, thè first was to derive thè interval 
distribution from a time dependent rate. I studied thè problem and approached 
with classical stochastic instruments without effective success, then I accepted to 
relax some requirements and obtained an interval distribution that is an effective 
solution for thè needs of thè simulator, reproducing correctly thè results of thè BEST 
model. Eventually thè most of analytical effort regards thè ideatimi, and resolution, 
of a SDE systems. It was devised to reproduce competition for criticai resources. 
The model was thè last of three models which were not analytically solvable. These 
model are sketched in thè thesis but are not deepen since, for thè state of thè art 
of NetGrowth, they are not relevant. I have enhanced my ability in approaching 
stochastic problems and, most of all, in evaluation of thè quality of thè proposed 
Solutions. 

Research perspective The NetGrowth simulator will be soon released for thè 
scientifìc community. In next future a serious study on thè performance will be 
accomplishes and potential flaws will be solved. Hence thè next step is thè systematic 
implementation of methods for synapses generation. Up to now thè intersection 
method is being used, which estimates thè numbers of synapses generated by thè 
number of intersections between axon and dendrites, this algorithm can be improved 
with density field methods. 

The implementation of fasciculation among neurites is a required improvement too. 
On longer terms we expect to develop a simulation environment (SENEC initiative)]^] 
towards thè realization of a growth-activity simulator. We aim to join NetGrowth 
with thè NEST activity simulator through thè NNGT library developed by T.Fardet. 
Either thè project misses a sound validation with experiment realised on purpose. 
For this scope, we pian to ideate and implement experiments with thè researchers of 
IPGG. 

Me, student in neuroscience As a theoretical physicist my knowledge in bio- 
physics and neurobiology was very limited. This research project was an opportunity 
to discover thè complexity behind neuronal outgrowth. During thè thesis I had thè 
opportunity to meet researchers with very different backgrounds, from people I met 
in Sissa Neuron Technology Summer School to thè researchers of IPGG. Discussing 
thè NetGrowth project with them opened my mind with respect to thè multi-faced 
problems of neuroscience. 

I do not think I will keep my way on thè biophysics of neuronal growth, but now I 

J https://github.com/SENeC-Initiative 




106 


know what is possible and what not in a in vitro culture experiment. This experience 
will help me in collaborating with neurobiologists and experimentalists. 



107 


Appendix A 

Random Walk characterization 


A.0.1 Random Walk, a short hint 

Random Walk approach is applied in many fields and thè repertoire of problems is 
pretty huge, we will follow some works on Random Walk in biology [?], trying to 
characterize results from in vitro experiments and from NetGrowth simulator. 

The study of long chain of molecules, as polymers, is a lucky application of Random 
Walk approach. It’s relevant to study their properties depending on thè micro physics 
of thè anchorage. Another relevant approach descends from statistical mechanic and 
shades thè light on thè macro scale behaviour: a sub-chain of strongly connected 
elements, can be considered under a coarse grain operation, as an unit, hence 
thè constraints between different units are weaker and thè long range behaviour 
presents other properties than short range. An important objective, arnong others, 
is to describe thè random walk in a statistical field theory formalismi, studying thè 
properties of thè polymer chain with thè renormalization group. In such formalism 
thè Mean Square Displacement, i.e. thè absolute distance between then start and 
thè end of thè path, is wjv, and N is thè nunrber of segments. 


V 


lim RT 

AT->oo 


(A.l) 


is fixed finite and thè correct exponent 7 is investigated [Procacci et al., 2 008 . The 
algorithms we are implementing are short-range models, since thè partial correlation 
function is an exponential decaying one, and can be demonstrated they ultimately 
belong to diffusive random walkers as defined in ??. The literature we are referring 
to is generally related to polymers. [Cioletti et al., 2014]. 


A.0.2 Characterization of Random walks: tortuosity, contraction 
and MSD 

To evaluate thè morphology produced by thè different algorithms we need to search 
for some statistical properties: Measures offering relevant information of thè studied 
System too, i.e. thè distance reached in a certain time by thè extruding growth 
cone. All thè statistical measures are computed over thè reai distance from thè soma, 
measured in micrometers. The assumption rnade throughout this Appendix is thè 
process is not stationary: thè dynamics is not governed only by relative distances 







108 


and thè origin time has any special meaning. This approximation is relaxed since 
we expect a strong dependence by time for Growth Cone dynamics. 

On thè other side we assume thè System will relax towards a diffusive Random 
Walk and thè persistence length obtained for thè first interval of thè neurite is thè 
longest possible during neurite evolution. That is, we expect thè GC dynamics will 
be influenced less by thè soma position as far it is from it. 


Square distance from thè origin 

The first element we are evaluating is thè distance reached by thè Growth Cone. 
The neuron extrudes its active units, thè GC, and they will tread a certain amount 
of distance before attach to a target dendrite. Understand this is a diffusive or a 
ballistic behaviour is a cruciai point. In this paragraph thè length of thè segments is 
considered unitarian and homogeneous throughout thè process. Roughly, we can 
divide thè dynamics into two intervals, a super-diffusive (anomalous) and a diffusive 
one. The former ends when thè distance travelled by thè growth cone is far larger 
than thè persistence length of thè System. In case of reai neurons, if presents, thè 
diffusive lirnit is barely reached, as can be seen in Fig. |4.7| To simulate thè growth 
pattern of reai neurons, it’s required to know thè persistence length of such neurons 
and reproduce it with thè random walker. 

In thè 2D RW thè absolute distance from thè initial point is: 


(AX 2 ) = ((x - x 0 ) 2 + (y - yof 


(A.2) 


thè ratio over thè curvilinear abscissa can be linear(o: = 1 ), diffusive (a = 2 ), or 
present anomalous diffusion with higher powers: 


(AX 2 ) oc t° 


(A-3) 


This simple measure is enough to characterize thè persistence length of thè System, 
recalling thè work of|Tchen, 1989 , assuming stationary state, thè MSD is: 


(x 2 ) = Y J {hb j ) = b 2 Y J Y, c ^-j) 

0 0 0 

Neglecting thè ends thè previous equation can be transformed in 


(A.4) 


! ^^C(s) + ò 2 ^C(0) 

s=0 j=s j =0 

n 

(A.5) 

S = Y,C(s) 

s=1 

(A.6) 

n 

M = J2 sC(s) 

s —J 

(A.7) 


(A-8) 


where S and M can be interpreted as thè volume below thè correlation curve and 
its first moment. 
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Thus, defìned = nb 2 thè diffusive MSD. It’s possible thè obtain thè mean square 
displacement: 


r 2 M 

(-j) = 1 + 2(5 -— 

?’5 n 

for large n: (A.9) 

r 2 

(“ 2 ) = 1+25 

r o 

(A.10) 


We have verified this relation for thè MCRW in Fig. ??, The volume S and thè 
moment M are computed numerically. 

For an exponentially decaying correlation function states: 


5 = l p , M = l 2 p 


The relation between thè memory parameter a and thè persisterne length should be 
valuated frorn thè partial correlation coefhcient = (bibi + k), indeed it’s possible to 
demonstrate that for a finite nurnber of partial correlation: 


(r 2 /r’o 


nj = i(l + Cj) 

n"=i(i-c,) 


1+2 ip 


(A.ll) 
(A. 12) 


But nevertheless we know thè partial correlation is an exponentially decaying is stili 
hard to compute it analytically for thè MCRW. 

Some attempts follow in next paragraph. In next section we will show thè relation 
between thè MSD here presented and thè statistica! evolution of angle variable 9. 


Mean Square Displacement of angle variable: 

The second relevant measure, linked with previous is thè Mean Square Displacement 
of 9 over thè distance s. 

A 9 2 = (9 S - 9 0 ) 2 = 9 2 s 


In this case it’s easy to derive a consistent relation with thè persistence length for a 
short interval. Asserted thè global correlation is an exponential dumping function: 
C(s) ~ exp (-l/lp). 

For l <C l p 


C\s)^l-^+0{l- 2 ) 

Lp 

(A-13) 

C(.) -1 - <A f > 

(A. 14) 

« l~A 

CN 

II 

CM 

< 

(A-15) 
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Figure A.l. Preliminary study for exponential decaying memory algorithm 

a. The MCRW algorithm, described by Eq. 4.1.3| is executed varying a, then thè 
persistence length is measured as described in |4.4| 

b. The persistence length as a function of a and a, thè graphic is fitteci by A * exp(a * £) 


but it sliows a super-exponential behaviour. c. Eq. The convergence of |A.0.2|is verifìed: 
red dots correspond to wliile blue dots are computed through thè integrai of thè 
correlation curve. 

d. The System converges to linin^oo R/R0 faster with lower a. The equilibrami state 
can be reached in a very long time when a approach to one, and in generai reveal thè 
presence of a large transient. In thè transient thè memory constrained CRW presents a 
non diffusive behaviour. 
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This measure naturally defines thè persistence length as thè inverse of thè slope: 
l p = m^e 2 ) • This approach is supported by thè article [Martin, 2013| which uses 
this approximation to estimate thè persistence length of microtubules. 

To be consistent about persistence length and correlation time let’s reformulate 
it taking into account thè variable length of each segment, thè previous equation 
becomes 


< A92 <"» = 2 p 

where b is thè length of thè segment. 


(A.16) 


Contraction 


The contraction is a non standard statistical measure, often encountered in neurons 
and generai biological morphology evaluations. It’s dehned as thè average ratio 
between thè euclidean distance of thè growth cone from thè soma and thè length of 
thè neurite itself; 7] is always less than 1. It’s showed to be an exponential decay for 
thè CRW we are considering and thè r is obtained. The Contraction curve inforni 
us whether thè diffusion limit was reached or not, indeed thè ratio will change, in 
exponential manner, until thè neurite length is much longer then thè persistence 


length, then will establish on a fixed value which correspond to Eq. A.0.2 


A.0.3 Analytic approach to MCRW 


Our interest remains to obtain thè persistence length from thè parameters in Eq. 4.1.3 
Let’s make some consideration on this function: 


• It loses of physical sense when thè angle get larger than n, 

• Taylor expansion of Eq. |4.4| requires n remains little. 

• To simplify thè computation thè limit of a 11 —> 0 would be useful. 

To restore limit a n ->0a simple consideration can be offered: for large a we don’t 
need to be very dose to thè origin, since thè angular variation will require larger 
distance to affect thè Taylor expansion, while for little a thè correlation will reduce 
to zero in few steps and thè limit is valid. Let’s look at thè average values of thè 
stochastic process, keeping a n for exactness. 


@n +1 — @n +1 ~b & r n+l 


We note that (A.0.2) can be rewritten: 

1 — a 


@n +1 — 


1 - a n+l 


^Ta n k 9 k + ar n+ i 


k=0 


(A.17) 


(A.18) 
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Recursive relations For 9: 



(0 n - ar n ) + ar n+ì 


1 + a( 1 — a) — a 



1 - a n+l 
let 


1 + a( 1 — a) — a n 
1 - a n+1 
1 - a n 


~ a (r n+ 1 - r n ) 


(T y n+1 i — a n+i Tn 


For 9 2 : 


0n+ 1 = £ 2 # 2 + 2az (r n+ i + r n ) 6 n + <r 2 (r n +1 + r„) 2 


Average values If we average, recalling r is a Gaussian with mean zero 



[1 + a(l — a)] ( 9 n } 


this is a sirnple geometrie suite that leads to: 


(9 n ) = [1 + a(l - a)] n 9 q 


(A.19) 


For thè squared value: 


{9 2 n+ i) ~ 1 + z 2 (9 2 ) + 2az n {(r n+ i - r n ) 6 n ) + cr 2 ((r n+1 - r n ) 2 ) 
~ z l(°l) + 2 ^n((l - P)<?rl) + <x 2 ((r 2 +1 - 2r n+ i + r 2 )) 
~ z n(9n) + 2ct"(1 — j3)z n + 2cr^(l — /3) 

~ z n(9n ) + 2cr(l — (3)(z n + 1) 

~ fln(^n) 9-n 


This equation is exact and can be rewritten starting from 9 q: 


(9i,+ 1 )— n o,ì{9q) +^ n ak 


(A.20) 


i i k=1 


we should linearize thè previous expression for little n and relate to thè persistence 
length through |A.16 since thè variable n is discrete thè simplest way is to evaluate 
thè difference between thè MSD in n = 1 and n = 0, and setting thè initial angle to 
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zero: 


d_ 

dn 


(^n+l) l n =l — 


1 



(oi(ao(0“) + bo ) + òi) — (ao(0 2 ) + &i) 

= 61 + &o«o - b 0 = bi + bo(ao - 1) = &(zi + 1) + 6(z 0 + l)(a 0 - 1) 
= b(zi - z 0 + Zq) 


2a\l-f3) 


1 + a — (1 — a) 2 / 1 — a(l — a) 

1 — a 2 1 — a 


31 


which is wrong. Indeed thè a dependent behaves like an exponent and thè persistence 
length goes to zero for high values of a, which is heuristically and numerically verihed 
incorrect. 

The precious computation remains useful to study thè long range behaviour of S 2 . 
It’s calculated that for large n there is a arithmetico-geometrie suite which can be 
rewritten as: 

<0-7 = [1 + a(l — a )] 2 (<0 2 ) — 7 ) (A.21) 


with 


b 2 er 2 (l — f3) [2 + a( 1 — a)] 

1 — a a( 1 — a) (2 — a( 1 — a )) 


(A.22) 


hence: 

<^n+ 1 ) = [1 + a(l - a)] 2n (<# 2 )o - 7 ) + 7 (A.23) 


This result will be useful in later section when considering thè angle distribution on 
a longer tirne scale. 

Another approach that can be attempted is compute it numerically: sagemath [j] a 
powerful python tool for symbolic computation has been used for it. 6 2 has been 
computed up to thè third terni, and then thè linearization performed, subtracting 
thè second terni. Hoping thè function maintains thè same rate in thè transient. 


(Ad 2 ) = (d 2 ) - (9 2 ) 

= (5a 6 + 10a 5 + 2 (3a 6 + 9a 5 + 13a 4 + 12a 3 + 7a 2 + 3a + l)) /3 3 
—2a 3 - 2 (2 a 6 + 4 a 5 + a 4 - 4 a 3 - 3 a 2 ) j3 2 
-2 a 2 + 2 (3 a 5 + 5 a 4 - 3 a 3 - 5 a 2 - 3 a - l) /? 

a 2 

a 6 + 4 a 5 + 8 a 4 + 10 a 3 + 8 a 2 + 4 a + 1 

Neither in this case thè function l p (cx, f3, a) presents thè right behaviour and further 
attempts to perform an analytical estimatimi of thè persistence length won’t be done. 


Angle Mean Square Displacement and Persistence Length In this section 
we have tried to measure thè persistence length with numerical simulations. To 
reconstruct thè proper function, thè one to convert thè requested persistence length 


1 www. sagemath. org 
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into thè right adimensionai parameters, let’s tackle thè problem one parameter at 
time. First l p is obtained as function of a. in a purely memory random walk. This 
function is illustrated in Fig. ?? to be super exponential, but silice we are interested 
in moderate values of l p we can accept thè error of considering an exponential 
function. Then another fit changing a was attempted, and thè following is thè 
formula obtained. 


l p = 7 exp(—5cr + 12a)o: = (ln(Z p /7) + 5<r)/12 (A.24) 

Despite thè expectations thè implementation of such formula didn’t work, indeed 
thè fit for a it’s to rude and a little change in its value causes a great variatimi in l p . 
The valid alternative was to fit over a 3D manifold, but such computation requires 
sophisticated approaches, especially if thè functional form is unknown. To track 
values usable in thè project have we searched thè subset of parameters corresponding 
to a plausible persistence length, a range between 50 um and 300 um. The table is 
resumed in a 3D plot in Fig. |A.2| 
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(b it bj) for max memory coeff. 
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Figure A.2. Plausible values for a, fi and a 

Here we present thè set of parameters offering a plausible persistence length as result, 
in thè range l p £ (50,300)/uto The grid search was performed in this manner: fìrst thè 
correlated Gaussian parameter fi and thè angle of view a were fìxed, then thè memory 
parameter a was varied from 0, up to obtain thè maximum persistence length. < bibj > 
was measured through thè cosine oi Ad. 

(A-B) show thè persistence length reached for each simulation thè yellow stripes indicate 
that was impossible to set a proper a for that set of fi, a. In (A) thè maximum values of 
correlation length reached during thè grid search. In B only a stripe show values about 
50 firn: for some parameters even very short memory, a ~ 0, presents a persistence 
length longer than thè required interval. (B-C) show thè max and min values of a 
for thè given persistence length. The graph show thè importance of a in setting thè 
persistence length. The valid values of a are between zero and one, but wlren thè a is to 
high (0.6), any values is enough to keep thè persistence length in thè required range. 
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Appendix B 

Correlateci Random Walk: 
Coloured Noise 


The author in |Deserno, 2011] offer a simple algorithm to generate correlateci Gaus- 
sian, with a correlation coefficient / = exp(—1/r): 


r n = f ■ r n -1 + \Jl ~ P ■ in (B.l) 

f o = Co (B.2) 

C(m - n) = (r n r m ) = f m ~ n (B.3) 


In thè same paper such System is discusseci: Once defìned thè correlated and 
uncorrelated random walks as: 


N 

Rn ■■= Y r n 

71=0 

N 

Gn := Y, tn 

n =0 


it’s easy to show that: 


Rn = 


Rn = 


' 1 + / 

1 -/ 


G n ~ 


f 


-r N + l- 


coth — Gn 
2 T 


x/W 2 

i — / JV 1 ‘ i-/ 

Gn : t <s 1 

\J~2tGn :t>1 


(B.4) 

(B.5) 


(B.6) 

(B.7) 


While thè second and third term in Eq |B.6| will be of Constant magnitude with time 
thè fìrst term will reflect thè behaviour of a classic Random Walk; hence it can be 
seen that thè correlated Gaussian random walk is “faster” than thè uncorrelated 
random walk, even though in both cases thè increments are Gaussian deviates with 
zero mean and unit variance! 

More precisely, thè mean square displacement of thè correlated walk will grow 
stronger - in thè case r > 1 - by a factor of 2 t: Correlation gives distance! This 
System can be described as Markovian if thè space of phase is enlarged to thè 
parameter r n , and this procedure is thoroughly generai. 
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Appendix C 

Fokker Plank 


Fokker-Planck 

Let’s thè variable A evolve by thè process: 

À = — ( Am — A) + £ 

T 

Let f(A, t ) be thè probability density of A at time i; we have: 

f(A, t + At)= j f(A - A', t)P At (A% A - A')dA' 

= /(A, t) - fd A (AA) - (A A)d A f + l -fd\{{AA) 2 ) 

+d A fd A ((AA) 2 )+ l -((AA) 2 )d 2 A f 


Hence: 


$/ = --/ + 


A 


r A 


a: 


- — ) d A f + ^d 2 A f = d A 


1 


( Am — A)f + -^-d A f 


with thè boundary condition, since A > 0, 

Vi, d A f(0,t) = 0. 


(C.l) 

(C.2) 

(C.3) 

(C.4) 

(C.5) 

(C.6) 


Stationary regime 


Once thè permanent regime is reached, thè Fokker-Plank equation (C.5) becomes: 


d A 


{Am — A)f + ~^~d A f 


T 


= 0 


1 07 

Hence —(Am — A)f -\ — -d A f = 0 (regularization at A 
t 2 

— (A — Am) 2 + est, which leads to: 


f(A) = C exp 


(A — AmY 


(C.7) 

od). Thus ln(/) = 

(C.8) 
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where ctj = ^. C is thè normalization coefficient given by thè integrai of thè 

probability distribution to 1 between 0 and oo: 


C~ x = 



(C.9) 

(C.10) 


Here a reference to Heun http://mathfaculty.fullerton.edu/mathews//n2003/heunsmethod/Heun’sMetl 
* Neuron, neurite and growth cone properties 
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Appendix D 

Variables of thè simulator 


//! number of neurites for createci neuron [int] 
extern const std::string num_neurites; 

#defin SOMA_RADIUS // for neurite starts 10 [micrometers] 
extern const std::string soma_radius; 

#define AX0N_DIAMETER 0.6 //[micrometers] 
extern const std::string axon_diameter; 

#define DENDRITE_DIAMETER 0.6 //[micrometers] 
extern const std::string dendrite_diameter; 

#define INIT_BRANCHING 0 //initial branching length 0 [micrometers] 
extern const std::string initial_branch_length; 


/* 

* SPACE SENSING MODEL 
*/ 


extern 

const 

std : 

:string 

extern 

const 

std : 

:string 

extern 

const 

std : 

: string 

extern 

const 

std : 

: string 

extern 

const 

std : 

: string 

extern 

const 

std : 

: string 

extern 

const 

std : 

: string 

extern 

const 

std : 

:string 

extern 

const 

std : 

:string 

extern 

const 

std : 

:string 

extern 

const 

std : 

:string 


duration_retraction; 
filopodia_min_number; 
filopodia_finger_length; 
filopodia_wall_affinity; 
max_sensing_angle; 
proba_down_move; 
proba_retraction; 
scale_up_move; 
sensing_angle; 
speed_ratio_retraction; 
substrate_affinity; 


#define DURATI0N_RETRACTI0N 200. 

#define FIL0P0DIA_MIN_NUM 24 
#define FIL0P0DIA_FINGER_LENGTH 50. 

#define FIL0P0DIA_SUBSTRATE_AFINITY 0.1 
#define FIL0P0DIA_WALL_AFFINITY 2. 

#define MAX_SENSING_ANGLE 1.5707963267948966 // 100 degrees max for 1 s 
resol 

#define ONE DEGREE 0.017453292519943295 
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#define PROBA_RETRACTION 0.001 
#define PR0BA_D0WN_M0VE 0.008 
#define RW_DELTA_C0RR 100. 

#define RW_MEM0RY_TAU 100. 

#define RW_PERSISTENCE_LENGTH 10. 

#define SCALE_UP_M0VE 20. 

#define SENSING_ANGLE 0.1433 
#define SPEED_RATIQ_RETRACTIQN 0.2 
#define SPEED_GROWTH_CQNE 1. 

#define WALL_AFNTY_DECAY_CST 19.098593171027442 // inverse of 3 deg in 
radians 


/* 

* CRITICAI. MODELS 

*/ 

//! Oparam tub_topo_coefficient 0.1 [naturai] 
extern const std::string use_critical_resource; 
#define USE_CRITICAL false 
extern const std::string CR_use_ratio; 

#define CRITICAL_USE_RATI0 1 
extern const std::string CR_leakage; 

#define CRITICAL_LEAKAGE 6 

extern const std::string CR_correlation; 

#define CRITICAL_CQRRELATI0N 0.1 
extern const std::string CR_variance; 

#define CRITICAL_VARIANCE 0.1 // 

extern const std::string CR_weight_diameter; 

#define CRITICAL_WEIGHT_DIAMETER 1. 

extern const std::string CR_weight_centrifugal; 

#define CRITICAL_WEIGHT_CENTRIFUGAL 1. 

extern const std::string CR_elongation_factor; 
#define CRITICAL_ELQNGATION_FACTQR 0.5 
extern const std::string CR_elongation_th; 

#define CRITICAL_ELQNGATI0N_TH 0.35 
extern const std::string CR_retraction_factor; 
#define CRITICAL_RETRACTION_FACTQR 0.1 
extern const std::string CR_retraction_th; 

#define CRITICAL RETRACTION TH 0.15 


extern const std::string CR_neurite_split_th; 

#define CRITICAL_SPLIT_TH 250. 

extern const std::string CR_neurite_available; 

#define CRITICAL_AYAILABLE 100. 

extern const std::string CR_neurite_variance; 

#define CRITICAL_GEN_VAR 5. 

extern const std::string CR_neurite_generated; 

#define CRITICAL_GENERATED 150. 

extern const std::string CR_neurite_generated_tau; 
#define CRITICAL GEN TAU 100. 
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extern const std::string CR_neurite_delivery_tau; 
#define CRITICAL_DEL_TAU 50. 

#define CRITICAL GEN CORR 0. 


/* 

* RANDOM WALK MODEL 
*/ 

extern const std::string random_walk_submodel; 

//! Spararli speed_growth_cone 10 [micormeter/second] 
extern const std::string speed_growth_cone; 
extern const std::string speed_variance; 

//! Spararli persistenc_length 2000 [micrometer] 
extern const std::string rw_persistence_length; 
extern const std::string rw_memory_tau; 
extern const std::string rw_delta_corr; 

//Spararli sensing angle is choosen from experimental 
// data and it 1 s 8.2 degrees 

//! RUN AND TUMBLE 

extern const std::string rt_persistence_length; 

#define RT_PERSISTENCE_LENGTH 100. 

//SELF REFERENTIAL MODEL 

// 

#define SFR_AV0IDANCE_F0RCE 1 

extern const std::string srf_avoidance_force; 

#define SFR_AV0IDANCE_DECAY 2 

extern const std::string srf_avoidance_decay; 

#define SFR_INERTIAL_F0RCE 1 

extern const std::string srf_inertial_force ; 

#define SFR_INERTIAL_DECAY 2 

extern const std::string srf_inertial_decay ; 

#define SFR_S0MATR0PIC_F0RCE 1 

extern const std::string srf_somatropic_force; 

#define SFR_S0MATR0PIC_DECAY 2 

extern const std::string srf_somatropic_decay; 


/* 

* GR0WTH C0NE SPLITTING PARAMETERS 
*/ 

extern const std::string gc_split_angle_mean; 
extern const std::string gc_split_angle_std; 

//! Spararli van_pelt model for branching probability and direction default: 
True 

extern const std::string use_van_pelt; 

//! Van_Pelt BEST model parameters 
extern const std::string B; 



extern const std::string E; 
extern const std::string S; 
extern const std::string T; 

#define GC_SPLIT_ANGLE_MEAN 98.0 / 180 * 3.14 
#define GC_SPLIT_ANGLE_STD 10. / 180 * 3.14 
#define USE_VAN_PELT true 
#define VP_B 5. 

#define VP_E 0.05 
#define VP_S 1. 

#define VP T 0.01 


/* 

* LATERAL BRANCHING PARAMETERS 
*/ 


extern const 
extern const 
extern const 
extern const 
extern const 
extern const 
extern const 
extern const 


std: :string 
std: :string 
std: :string 
std: :string 
std: :string 
std: :string 
std: :string 
std: :string 


use_flpl_branching; 
flpl_branching_rate; 
use_uniform_branching; 
uniform_branching_rate; 
lateral_branching_angle_mean 
lateral_branching_angle_std; 
diameter_variance; 
diameter_eta_exp; 


#define ANGLE_IN_DEGREES true 
#define DIAMETER_ETA_EXP 1.5 
#define DIAMETER_VARIANCE 0.1 

#define LATERAL_BRANCHING_ANGLE_MEAN 90 * 3.14 / 180 
#define LATERAL_BRANCHING_ANGLE_STD 1. / 180 * 3.14 
#define UNIF0RM BRANCHING RATE 0.001 


/* 

* ACTIN WAVE MODEL 
*/ 

//! actin wave trigger or not False [bool] 
extern const std::string use_actin_waves; 
#define USE_ACTIN_WAVES false 
//! Actin Waves model parameters 
extern const std::string actin_content; 

#define ACTIN_C0NTENT 0. 

extern const std::string actin_content_tau; 
#define ACTIN_C0NTENT_TAU -1. 
extern const std::string actin_wave_speed; 
#define ACTIN_WAVE_SPEED 150. 
extern const std::string actin_freq; 

#define AW GENERATION STEP -1. 
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* RECORDERS 
*/ 


extern 

const 

std : 

: string 

event_type; 

extern 

const 

std : 

: string 

interval; 

extern 

const 

std : 

: string 

level; 

extern 

const 

std : 

:string 

observable; 

extern 

const 

std : 

:string 

observables; 

extern 

const 

std : 

:string 

record_to; 

extern 

const 

std : 

:string 

restrict_to; 

extern 

const 

std : 

:string 

targets; 

extern 

const 

signed char 

lateral_branching; 

extern 

const 

signed char 

gc_splitting; 

extern 

const 

signed char 

gc_deletion; 

extern 

const 

std : 

:string 

num_growth_cones; 
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Appendix E 

Neuron morphologies standard 
format 

E.l Cultured neuron image analysis 

I will discuss briefìy about image acquisition from experiments and relative analysis. 
This appendix results by a document I provided to experimental neurobiologist in 
Septembers, to state our requirements for a effective image analysis and NeetGrowth 
validation. 

E. 1.1 Required measures 

To test and validate *Net growth* models morphological information from experimen¬ 
tal observation are necessary.. Silice each model concentrate on a different aspect of 
neuronal morphology development then there are varied possible measures, we can 
roughly divide into: 

Growth cone direction, elongation and guidance To study thè mechanical 
or chemical guidance of thè growth cone a detailed information on its path is 
required, in thè growth phase thè dynamic of thè GC is not expected to be straightly 
continuous and homogeneous, it’s expected to retract, pause and elongate with 
different interval. A sound comparison between thè simulator and thè experiment 
requires information on thè dynamics of thè GC. Generally these measure can be 
achieved with fixed camera in bright field, but thè quality won’t be satisfying; also 
can be obtained with timelapse of a fhiorescent-dye based illumination. NetGrowth 
can furnish this information with thè Recorders, they are stili not fully working but 
we expect we will be able to investigate and compare thè dynamics of thè Growth 
process. Generally thè dynamics can be stored in TIFF image format. 

Neurite branching and topology The shape of neurites requires to be matched 
with experiment through a graph-like structure. Some sirnple topological measures 
(number of internai nodes and leaves) can be performed with a merely topological 
information, as in Van Pelt and VanOoyen works. Other measures, like Sholl Analysis, 
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require geometrical information. The TREESToolbox |Cunt z~t al., 20l2] can help 
in this. 

E.1.2 Cultured neurons images 

Formats and Software In order to relate thè data generated from a simulator 
to those acquired through an optical instruments it’s required to develop a common 
language. Such a common language is, computationally speaking, a format. A 
format is a way to stock data which is standard and uniformly applied by all thè 
contributors. Whichever standard will provide for a finite set of information, this 
set has to be as large as possible to avoid any loss of details and information. 

In neuroscience thè process of digitalization has a 20 years backward history. 
One of thè rnost active researcher in this field is A. Ascoli 0 He played a main role 
in thè developing of Neuromorpho.ortj^the Online facility to store and share neuron 
morphologic data, developed at thè George Mason University. The NeuroMorpho 
(NM) staff works to maintain an updated database of 3D neurons, a research group 
can upload thè result of its research following thè application line of thè project. The 
file will be acquired and evaluated and then added to thè database. The standard 
format which NM utilizes is thè SWC, an example of stored neurons is Fig. ?? 

The standard format for experimentalist is thè TIFlQ (Tagged Image File 
Format) Tiff file can contain time-lapses, they are stored as a series of images on 
thè temporal dimension, like a video. Furthermore it can store different tagged 
channels, this is particularly relevant when thè sample is observed with more than 1 
fluorescent dye, e.g. one for tubulin, one for actin, one for calcium, and so on. Tiff 
format is commonly used by other communities (graphics, physician, etc... ) but 
it’s a proprietary format. The most common tool to visualize and obtain data from 
Tiff file is IrnageJ, ImageJ allows for channels and timelapses sequences. 

Drawbacks 

• Swc format has no information on time and dynamics, thè file is a shoot of 
thè neuron morphology at a precise instant during thè growth. 

• Tiff format has no usable information on neuron topology or geometry, all 
thè information is encoded in pixel. 

The NetGrowth team has decided to adhere NeuroMorpho standard format: 
SWC (Cannon et al., 1998).It is documetned at thè end of thè appendix. Since thè 
importance this format has in thè neuroscience concerning with morphology rnany 
tools where developed to translate proprietary format to SWC one: one of them is 
NeuronLanc^ To perforai analyses on thè neuron we are using btmorpì¥ ]tool. It was 
developed in thè framework of another simulator by Torben Nielser^j At present 

http://dx.doi.org/10.1016/j.neuron.2013.03.008 
2 http://neuromorpho.org/index.j sp 
'https://en.wikipedia.org/wiki/TIFF 
J http://www.neuronland.org/NL.html 

’http://b-torbennielsen.home.oist.jp/btmorph/index.html 
'http://b-torbennielsen.home.oist.jp/Software.html 
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Figure E.l. Dynamics or Morphology Top: Neuron image obtained from Swc data 
troni Netgrowth simulato, bottoni Example of Tifi single timelapse image, with curtesy 
of F. Ulloa Severino from Sissa Vincent Torre group 
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moment thè Btmorph github repository is in a messy state since many group were 
indipendently contributing to it. 

We developed a compatibility tool for Netgrowth and we hope keep working on 
it. 

Tiff to Swc converter In order to work together, us(NG team) and whoever 
want to use our NetGrowth to evaluate is experiments, it is required to have a 
software to convert TIFF files into SWC. Detect and recognize neurites branching is 
not an easy task and many tools are available on thè Internet, each one with prò and 
cons. A full reviews of this tool would require a long time, so we concentrated on 
some tools which are most promising and whose developers are available to internet. 
Maybe TREESToolbox is thè best. 

Here it is a list of tools: 

http://www.northeastern.edu/neurogeometry/resources/tools/ 
http://imagej.net/Simple_Neurite_Tracer#Installation 

https://wiki.colby.edu/display/BI474/ImageJ+Software+-+Neurite+Measurement 
https://github.com/sgbasel/neuritetracker 

https://imagescience.org/meij ering/publications/download/cyto2004.pdf 
http://onlinelibrary.wiley.com/doi/10.1002/dneu.20866/full 
http://fournierlab.mcgill.ca/styled-6/NeuriteTracer.html 


Probabilistic hypotesis density filtering 

We had a contact with thè developer of this ImageJ plugirj^] Until now thè swc file 
generated frorn thè program is unusable, mainly for thè presence of noise in thè file. 
The following improvement are necessary to use it: 

• The program should at least distinguali between different objects. Two neurons 
should not be in thè samefìle, or at least they cannot be attached each other 
in thè swc chain. Much more if thè attached segments are just noise. 

• The program should attach thè neurite segments each other and in thè proper 
order. I couldn’ check this is done since I cannot plot he SWC with btmorph. 

• Should distinguish between neurite and soma, and set thè appropriate identifier 

• The program should measure thè diameter of thè neurite and associate each 
point a dimension, even relative. 

• This is thè information obtained with phd from sanie image above. 


E.2 SWC format 

The three dimensionai structure of a neuron can be represented in a SWC format. 
SWC is a simple Standardized format. A swc neuron is a sequence of cylinders which 
diameter reflects thè neurite diameter in that point. Each elemnt has a parent which 

'https://bitbucket.org/miroslavradojevic/phd/overview 
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its connected. Every compartment has only one parent and thè parent compartment 
for thè fìrst point in each file is always ‘-1‘ (if thè file does not include thè soma 
information then thè originating point of thè tree will be connected to a parent of 
-1). The index for parent compartments are always less than child compartments. 
Loops and unconnected branches are excluded. All trees should originate from thè 
soma and have parent type 1 if thè file includes soma information. 

Soma can be a single point or more than one point. When thè soma is encoded 
as one line in thè SWC, it is interpreted as a "sphere". When it is encoded by more 
than 1 line, it could be a set of tapering cylinders (as in some pyramidal cells) or 
even a 2D projected contour ("circumference"). Each line has 7 fields encoding data 
for a single neuronal compartment: 

http://www.neuronland.org/NL.html 

1. an integer number as compartment identifierW 

2. type of neuronal compartmentW 

* ‘0‘ - undefined 

* ‘1‘ - soma 

* ‘2‘ - axon 

* ‘3‘ - basai dendrite 

* ‘ 4 ‘ - apical dendrite 

* ‘5‘ - fork point 

* ‘ 6 ‘- unspecified neurites 

* ‘ 7 ‘ - glia processes 

3. x coordinate of thè compartment 

4. y coordinate of thè compartment 

5. z coordinate of thè compartment 

6. radius of thè compartment 

7. parent compartment 
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